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ABSTRACT 


This  study  gives  a  means  of  classifying  synoptic  pressure  maps  for  Alberta 
using  an  objective  correlation  method  initially  developed  by  Lund.  It  also 
illustrates  the  development  of  a  climatological  or  forecast  aid,  namely 
probabilities  of  precipitation.  Classification  of  500  mb.  flows  for  the  summer 
period.  May  to  September,  1963-1971,  has  delineated  35  distinct  flow  patterns 
(Appendices  3  and  4).  A  year-round  Alberta  region  surface  map  catalogue  (33 
distinct  types.  Appendices  1  and  2)  was  also  developed  using  26  years  of  data 
(January  to  December,  1946-1971).  Both  catalogues  contain  statistical 
information  on  map  patterns  for  each  month,  so  they  are  valuable  climatological 
aids  in  themselves. 

Probabilities  of  24  hour  summer  precipitation  (06Z-06Z,  £  0.01  inches) 
for  500  mb.,  surface,  and  combinations  of  these  map  types  (May- September, 
1963-1971)  are  calculated.  They  illustrate  the  application  of  the  map  type 
catalogues  to  precipitation  forecasting.  The  results  indicate  that  significant 
variations  in  precipitation  probabilities  are  obtained  for  different  map  types  or 
combinations  of  types. 

The  precipitation  probabilities  derived  from  the  map  catalogue  data  are 
verified  for  an  independent  data  sample  (06Z-06Z,  May-September,  1972-1973) 
using  the  Brier  score.  The  results  indicate  that  the  forecast  probability  nearly 
equals  the  observed  frequency  of  measureable  precipitation.  However,  little 
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improvement  over  climatology  is  obtained  for  many  stations  because  there  are  too 
many  forecasts  near  the  climatological  frequency  of  precipitation. 

The  derived  probabilities  are  also  used  for  verifications  of  the  same  period 
(May  -  September,  1972-1973)  using  18Z-18Z  precipitation  amounts  rather  than 
06Z-06Z.  Comparable  results  are  obtained. 

The  results  using  the  precipitation  probabilities  are  good 
considering  that  only  two  predictors,  namely  the  500  mb.  and  surface  flows,are 
used. 

Two  final  studies  are  made.  The  first  consists  of  comparing  maps  for  a 
winter  period  with  the  summer  500  mb.  map  catalogue.  Indications  are  that  the 
500  mb.  summer  catalogue  is  useful  in  other  seasons.  The  second  study  tests  the 
ability  of  numerical  weather  prognoses  to  forecast  the  occurrence  of  the  correct 
map  type.  Numerical  prognoses  forecast  the  occurrence  of  the  actual  map  type 
about  fifty  percent  of  the  time.  With  improved  prognoses,  the  forecast 
verifications  will  undoubtedly  increase. 

The  500  mb.  and  surface  map  classifications  should  enhance  our 
knowledge  of  synoptic  pressure  systems  in  Western  Canada.  The  development  of 
climatological  and  forecast  aids  for  which  the  synoptic  flow  pattern  is  the 
dominant  factor  is  greatly  aided  by  map  catalogues.  The  use  of  these  map 
catalogues  in  a  forecast  office  may  prove  to  be  very  beneficial. 
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CHAPTER  I 


INTRODUCTION 


1 . 1  History 


Out  of  the  south  cometh  the  whirlwind 
And  cold  out  of  the  north  (Job  37:9). 

Fair  weather  cometh  out  of  the  north 
(Job  37:22) . 

When  the  east  wind  toucheth  it, 
it  shall  wither  (Ezekiel  17:10). 

Over  the  years  man  has  tried  to  correlate  weather  variables.  With 

the  coming  of  surface  maps  and  their  patterns  of  highs  and  lows,  other 

significant  predictors  of  the  weather  evolved. 

Winds  to  your  back,  storms  (lows)  on  your 
left  (Buys  Ballot's  Law). 


1.2  Synoptic  Climatology 
1.2.1  Background 

During  the  Nineteenth  Century,  the  relationship  of  wind  to 
pressure  systems  greatly  advanced  synoptic-scale  climatology  and  a  new 
foundation  was  provided. 

The  classic  approach  to  the  study  of  surface  patterns  has  been 
to  investigate  the  effects  of  the  surface  airflow  using  six  classes  of  isobars: 
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straight,  anticyclonic  or  cyclonic  curvature,  anticyclone  and  cyclone  centers 
and  indeterminate  situations.  The  emphasis  of  early  classifications  was  on 
patterns  that  remain  more  or  less  recognizable  for  several  days.  When  upper  air 
data  became  more  available,  classifications  were  carried  out  for  upper 
tropospheric  patterns  and  the  concept  of  large-scale  "steering"  of  surface 
features  was  introduced. 

Synoptic  climatology  developed  considerably  during 
and  after  the  Second  World  War .  Jacobs  (1946)  described  the  spatial 
distribution  of  certain  meteorological  elements  for  various  air  flows  over  Japan. 
He  suggested  that  at  that  time  the  Weather  Division  of  the  Army  Air  Forces  had 
better  knowledge  of  weather  processes  over  the  Japanese  Islands  than  did  the 
Japanese  themselves. 

The  large-scale  approach  through  similarity  of  patterns  for  a 
period  of  days  was  developed  about  the  same  time,and  is  discussed  by  Elliot 
(1951). 

Synoptic  climatological  studies  were  carried  out  for  many 
parts  of  the  world  in  the  1950's,  and  in  the  Soviet  Union  classifications  for  the 
entire  hemispheric  circulation  were  developed.  In  the  late  1950's  and  1960's 
objective  means  of  classifying  weather  maps  were  devised.  The  methods  and 
applications  of  synoptic  climatology  are  well  documented  by  Barry  and  Perry 
(1973)  for  this  period. 

Recently,  considerable  interest  has  been  regenerated  in 
weather  classification  schemes  because  of  the  availability  of  computers. 
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1.2.2  Classification  Approaches 

Various  ways  have  been  developed  in  which  pressure  maps  or 
the  weather  elements  themselves  may  be  classified.  Consider  first  an  individual 
map  of  the  pressure  field  (or  contour  field).  The  "static"  pattern,  may  be 
classified  into  categories  from  two  main  standpoints:  (  i)  the  identification  of 
individual  features  of  circulation  such  as  high  and  low  pressure  cells,  ridges  and 
troughs,  in  specified  locations;  (  ii)  the  description  of  the  complete  pressure 
pattern,  or  its  most  significant  features,  by  subjective  or  objective  methods.  The 
second  approach  covers  an  immense  range  of  methods  and  is  usually  referred  to  as 
"map  typing." 

The  pressure  (or  streamline)  map  may  also  be  regarded  in  a 
"kinematic"  framework,  and  two  broad  lines  of  work  of  comparable  importance 
have  been  developed  here.  They  are:  (  i)  classification  of  cyclone  and 
anticyclone  paths,  and  (  ii)  the  classification  of  circulation  or  air  flow  patterns. 
Examples  of  the  kinematic  approach  are  illustrated  in  the  works  of  Elliot  (1949, 
1951)  and  the  California  Institute  of  Technology  (1943).  Inevitably,there  is  some 
overlap  between  these  four  methods  in  practice  since  the  type  of  flow  pattern  is  a 
function  of  the  pressure  distribution  and  its  changes.  The  general  aims  of  all 
approaches  overlap  and  differences  between  them  are  essentially  a  matter  of 
emphasis. 

1.2.3  Pressure  Pattern  Classification  (Map  Typing) 

Surface  pressure  pcttern  classification  has  been  attempted 
since  the  turn  of  the  century  and  is  well  documented  by  Barry  and  Perry  (1973). 
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The  subjective  typing  of  pressure  patterns  presents  several  problems.  Since 
atmospheric  patterns  are  continuous,  the  delimitation  between  classes  is  somewhat 
arbitrary  and  therefore  rather  unsatisfactory.  Another  problem  arises  with  limiting 
the  number  of  map  types.  The  decisions  of  any  two  meteorologists  may  well  differ, 
so  that  a  unique  classification  cannot  be  obtained. 

Two  further  problems  are  often  excluded  from  a  map  typing 
scheme.  One  is  the  variable  intensity  of  pressure  systems  and  the  second  is  the 
seasonal  variation  in  type  characteristics. 

Much  thought  has  been  given  in  recent  years  to  numerical 
classification  techniques.  Three  major  approaches  have  been  used  -  correlation 
methods,  specification  techniques  and  empirical  orthogonal  functions  or 
eigenfunctions  (Barry  and  Perry,  1973). 

In  this  study,  the  correlation  method  is  applied  to  the  Alberta 
region  based  mainly  on  similar  studies  by  Lund  (1963)  and  Hartranft  et  al  (1970a, 
1970b).  A  complete  description  of  this  method  is  presented  before  Alberta  results 
for  500  mb.  and  surface  are  given.  Once  the  map  types  are  obtained,  they  are 
used  to  illustrate  the  development  of  a  climatological  and  forecast  aid,  i.e., 
probabilities  of  precipitation.  Forecast  verifications  are  included  for  independent 
data  samples.  A  summary  and  conclusions  are  offered  in  the  final  chapter. 


CHAPTER  II 


COMPUTERIZED  MAP  TYPING  AND  APPLICATIONS 

2.1  Introduction 

Recently,  considerable  interest  has  been  regenerated  in  weather  map 
typing  schemes  through  the  availability  of  computers.  The  map  pattern 
classification  technique  developed  by  Lund  (1963)  is  especially  noteworthy  since 
it  is  a  practical  and  efficient  means  of  objectively  identifying  recurring  synoptic 
map  configurations.  The  Lund  technique  demonstrates  that  by  utilizing  linear 
correlation  methods,  gridded  or  station  data  of  certain  meteorological  parameters 
(e.g.,  sea  level  pressure,  constant  pressure  heights,  temperature,  etc.)  can  be 
correlated  to  arrive  at  sets  of  similar  configurations.  This  technique  can  be  used 
either  to  classify  a  historical  period  of  data  into  various  surface  or  upper  air 
weather  types  for  subsequent  climatological  or  forecast  studies,  or  to  determine 
specific  map  types  associated  with  a  particular  weather  phenomena.  Because 
of  the  apparent  potential  of  the  Lund  technique  (Lund,  1963;  Hartranft  et  al, 
1970a),  the  Map  Typing  Project  Office  (12th  Weather  Squadron,  Colorado 
Springs)  has  been  applying  the  procedure,  with  extensive  refinements  and 
modifications,  to  develop  a  fully  computerized  system  of  classifying  day  to  day 
map  patterns. 
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2.2  Map  Pat-fern  Classification  by  Correlation  Methods 
2.2.1  Correlation  Coefficients 

The  formula  for  the  product  moment  linear  correlation 
coefficient  (r)  for  n  pairs  of  values  (in  this  context,  n  points  on  a  map)  is: 


E  lx, -  xHy,-  y} 

i :  1 


{  Z  [X-X]2 

i :  1 


where  the  X.  are  the  n  values  for  a  given  day  and  the  Y.  are  the  n 
corresponding  values  for  another  day.  The  results  of  correlating  each  element  of 
a  map  with  all  the  other  elements  of  a  set  of  maps  is  a  matrix  of  correlation 
coefficients  with  the  diagonal  values  equal  to  1 .00  (each  map  correlates  with  itself 
1 .00).  The  lowest  possible  correlation  coefficient  is  -1 .00  (when  two  maps  are 
completely  opposite,  i.e.,  highs  are  replaced  by  lows  and  vice  versa). 

The  effect  the  above  formula  has  on  r  when  comparing 
various  map  patterns  (isobars  or  contours)  can  be  illustrated  mathematically  or  by 
using  specific  examples.  The  theoretical  derivation  becomes  tedious  even  for  the 
simplest  patterns,  therefore  specific  examples  computed  by  Scholefield  (1973)  will 
be  discussed.  In  this  paper,  nine  hypothetical  weather  maps  were  constructed, 
each  being  composed  of  isobars  that  were  concentric  about  the  center.  Isobars 
were  drawn  at  four  millibar  intervals  for  all  maps  and  pressure  values  at  twelve 
grid  points  were  abstracted  and  computer  correlations  between  all  maps  were 
made.  The  results  obtained  are  summarized  on  the  following  page. 


7 


(  i)  Maps  correlate  perfectly  with  themselves  at  a  linear  correlation  coefficient 
of  1.00. 

(  ii)  Maps  that  have  identical  isobaric  configurations,  directions  of  flow  and 
horizontal  pressure  gradients,  will  correlate  at  a  value  of  1.00  regardless 
of  the  pressure  magnitudes,  provided  it  is  constant  at  any  respective  point. 
In  other  words,  there  would  be  no  difference  in  the  correlation  coefficient 
between  a  low  of  central  pressure  9 75  mb.  or  1009  mb.  as  long  as  the 
pressure  gradients  are  the  same. 

(iii)  Maps  that  have  identical  isobaric  configurations,  directions  of  flow  and 

different  but  constant  horizontal  pressure  gradients.,  correlate  identically. 

(  iv)  Maps  that  have  identical  isobaric  configurations,  direction  of  flow  and 
different  but  variable  horizontal  pressure  gradients  do  not  correlate 
identically  but  are  still  very  similar.  The  degree  of  similarity  between  such 
maps  decreases  with  the  degree  of  gradient  variability. 

(  v)  Maps  that  have  identical  isobaric  configurations  and  horizontal  pressure 
gradients  but  completely  opposite  flow  patterns  correlate  with  each  other 
at  -  1  .00  and  are  therefore  completely  opposite,  i  .e.  a  high  in  the  center 
compared  to  a  low. 

Considerable  insight  into  how  the  linear  correlation 
coefficient  equation  works,  when  used  to  correlate  surface  weather  map  patterns, 
is  given  from  the  above  study.  The  examples  used  shr  v  that  the  degree  of 
correlation  between  maps  is  almost  entirely  dependent  upon  the  isobaric 
configurations  and  the  directions  or  flow  of  the  maps  being  correlated.  Various 
horizontal  pressure  gradients  have  no  effect  on  the  correlation  as  long  as  the 
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gradients  on  each  map  are  constant.  If  they  are  variable,  that  is,  stronger  on 
some  portion  of  the  map  than  on  others,  the  correlation  between  them  will  be 
reduced . 

This  means  that  if  a  large  number  of  maps  is  being 
correlated  to  arrive  at  a  set  of  similar  map  types,  each  map  type  subset  will  be 
composed  of  maps  with  similar  isobaric  configurations  and  directions  of  flow.  It 
will,  however,  be  quite  possible  and  very  likely  that  maps  of  a  similar  type  will 
have  widely  different  horizontal  pressure  gradients.  Because  the  gradient  is 
closely  related  to  the  intensity  of  a  pressure  system,  it  becomes  significant  when 
one  tries  to  associate  certain  weather  conditions  with  map  types.  Therefore,  as 
a  second  approximation,  further  consideration  should  be  given  to  typing  maps 
according  to  the  pressure  gradient  strength  in  an  effort  to  increase  the  degree  of 
association  between  map  types  and  various  weather  conditions. 

2.2.2  Threshold  Correlation  Coefficient 

Once  all  maps  are  correlated  with  each  other,  there  are  a 
number  of  ways  to  analyze  the  data  and  select  types.  Type  selection  depends  upon 
the  choice  of  a  threshold  value  of  the  correlation  coefficient  which  is  used  to  divid 
the  sample  of  maps  into  groups  having  a  degree  of  similarity  in  pattern.  The 
question  arises:  At  what  level  of  correlation  does  the  similarity  betv/een  pressure 
patterns  cease?  It  has  been  found  that  if  the  value  is  too  low,  the  data  sample 
will  be  divided  into  fewer  map  types,  containing  a  large  number  of  maps  in  each 
type.  The  map  types  would  be  rather  general  and  there  would  be  wide 
variations  in  the  detailed  circulation  patterns  within  each  type.  With  each 
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increase  in  the  threshold  value  of  the  correlation  coefficient,  there  is  an 
accompanying  increase  in  the  number  of  types  generated.  This  is  because  each 
map  type  would  contain  maps  with  a  higher  degree  of  similarity  and  therefore  it 
requires  more  types  to  depict  the  range  and  orientations  of  map  patterns  which 
occur  over  the  area  of  interest.  More  significantly,  an  increase  in  the  threshold 
correlation  coefficient  yields  a  large  increase  in  the  number  of  uncorrelated  cases, 
i.e.,  the  number  of  maps  which  do  not  correlate  with  any  type. 

Lund  (1963)  and  Bjorem  (1966)  used  0.70  as  a  threshold  in 
their  surface  map  typing  procedure;  for  upper  air  patterns,  Augulis  (1969,  1970a, 
1970b,  1970c)  used  0.80.  The  Map  Typing  Project  Office  experimented  with 
several  different  threshold  values  (0.70,  0.80,  0.90).  The  practical  value  for 
typing  upper  air  maps  was  found  to  be  0.90,  while  a  value  of  0.70  was  the 
optimum  for  typing  surface  maps  (Hartranft,  1970a).  Later  tests  and  experience 
in  Surface  map  typing,  indicated  that  0.80  was  the  practical  threshold  value 
(Nieman  and  Sabin,  1972).  The  results  were  based  on  the  geographical  size  of 
the  area,  similarity  of  weather  for  stations  within  the  area  and  similarity  of  maps 
within  each  map  type. 

There  is  still  no  overall  consensus  on  the  choice  of  a 
threshold  correlation  coefficient.  Proponents  of  the  high  correlation  coefficient 
(e.g.  0.90)  want  better  discrimination  for  the  relatively  rare  events.  Conversely, 
the  main  argument  for  a  lower  correlation  coefficienf  (e.g.  0.70)  is  to  limit  the 
total  number  of  map  types  to  include  more  cases  in  ecch  type  and  increase  the 
reliability  of  the  climatological  statistics  generated  for  each  type. 


Some  authors  have  suggested  the  use  of  two  thresholds  in 
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determining  types.  One  is  to  use  a  high  value  to  derive  the  types  initially  and 
then  use  a  lower  threshold  to  assign  cases  to  the  types.  Another  method  is  to  use 
a  high  threshold  initially,  then  those  cases  which  remain  uncorrelated  would  be 
correlated  at  a  lower  threshold  to  obtain  additional  types.  This  would  result  in 
two  types  of  classification  -  strong  and  weak. 

The  selection  of  a  threshold  value  of  the  correlation 
coefficient,  therefore,  requires  careful  consideration  since  this  is  the  control  on 
how  similar  the  map  patterns  included  in  a  given  type  must  be.  The  choice  will 
depend  upon  the  map  size  being  used,  upon  the  type  of  data  field  being 
correlated  and  upon  the  kind  of  climatological  guide  being  developed. 

2.2.3  Map  Typing  Area  Windows) 

Lund  (1963)  selected  a  window  which  included  surface 
pressure  data  for  twenty-two  irregularly  spaced  reporting  stations  in  the 
northeastern  United  States,  map  size  area  approximately  ten  degrees  latitude  by 
thirteen  degrees  longitude.  Bjorem  (1966)  used  different  windows  for  his 
surface  and  500  mb.  map  types  in  the  northwestern  United  States.  Bjorem's 
surface  map  type  area  included  data  from  nineteen  irregularly  spaced  reporting 
stations,  approximately  ten  degrees  of  latitude  by  thirteen  degrees  longitude  , 
while  his  500  mb.  map  types  used  seventeen  grid  points  from  a  data  grid, 
approximately  fifteen  degrees  latitude  by  fifteen  degrees  longitude. 

Augulis  (1969,  1970a,  1970b,  1970c),  used  gridded  data 
for  500  mb.  from  the  National  Meteorological  Center  (NMC)  octagonal  grid  as 
the  data  base  and  included  fifty-two  data  grid  points,  approximately  thirty 
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degrees  latitude  by  fifty  degrees  longitude  .  This  large  map  size  was  selected  to 
emphasize  large-scale  flow  features,  i.e.,  long-wave  and  jet  stream  positions. 
Minor  features  such  as  migratory  short-wave  troughs  and  ridges  were  weighed 
lightly  in  their  precipitation  climatology  studies.  It  should  be  noted,  however, 
that  excessive  map  area  sizes  introduce  some  interpretation  problems  when 
attempting  to  apply  the  map  types.  Small-scale  features  of  the  circulation  will 
affect  only  a  relatively  small  number  of  data  points  within  the  network  and 
consequently  will  have  a  negligible  effect  on  the  overall  correlation.  Therefore, 
one  would  expect  to  find  two  maps  which  are  highly  correlated  but  which  are 
highly  dissimilar  in  one  or  more  portions  of  the  network. 

Initial  work  at  the  Map  Typing  Project  Office  in  Colorado 
used  gridded  data  of  surface  and  7C0  mb.  comprising  six  NMC  diamond  shaped 
grid  points,  with  an  additional  six  grid  points  obtained  by  interpolation  (tv/elve 
points  in  all).  The  map  size  covered  an  area  with  dimensions  of  ten  degrees  of 
latitude  and  fifteen  degrees  of  longitude.  In  more  recent  work,  map  catalogues 
were  produced  for  selected  areas  or  "windows"  in  the  U.S.A.  These  windows 
included  from  nineteen  to  twenty-four  NMC  grid  points  spaced  at  381  kilometers. 

From  the  above  discussion,  it  can  be  seen  that  the  selection 
of  the  map  size  or  "window"  for  a  given  data  network  is  a  subjective 
consideration  which  is  probably  best  resolved  by  experimentation.  There  are, 
however,  important  factors  to  be  considered  when  determining  the  optimum  map 
typing  area,  and  these  are  summarized  below: 

(  i)  Grid  or  data  spacing  -  Area  must  have  a  sufficient  number  of  data  points 
to  adequately  distinguish  the  flow. 
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(  ii)  Problem  to  be  solved  -  Size  of  area  depends  on  whether  the  problem  is 
global,  synoptic  or  mesoscale. 

(iii)  Size  of  correlation  coefficient  -  If  a  large  correlation  coefficient  plus  a 
small  number  of  types  is  desired,  a  small  area  must  be  used.  These  three 
factors  are  inter-related. 

(  iv)  Persistence  -  Grid  size  must  be  small  enough  so  that  a  map  type  does  not 
persist  for  too  long  a  period.  Changes  in  pattern,  not 
persistence,  are  preferred,  at  ieast  in  zonal  westerlies  studies. 

(  v)  Computer  limitations  -  The  computer  programs  and  computer  may  be  limited 
as  to  the  number  of  data  points  and  the  area  used. 

(vi)  Period  of  Forecast  -  Area  should  be  sufficiently  large  to  describe  the  flow 
pattern  that  will  govern  the  weather  during  the  forecast  period. 

2.2.4  Basic  Analysis  Technique 

Each  map  describing  the  synoptic  pattern  is  correlated  with 
all  other  maps  in  the  data  sample  using  equation  1 .  A  count  is  made  among  all 
the  maps  to  determine  how  many  maps  correlate  at  the  threshold  correlation 
coefficient  or  greater  with  each  map  in  the  sample. 

After  counting  maps,  a  choice  has  to  be  made  in  selecting 
"typical"  weather  map  types.  The  method  described  by  Lund  (1963)  and 
adopted  by  the  Map  Typing  Project  Office  is  to  determine  the  map  which  has  the 
largest  number  of  other  maps  correlating  with  it  at  or  above  the  threshold  value 
and  designating  it  Type  1 .  Following  the  Lund  method,  the  next  step  is  to 
remove  from  the  sample  all  of  the  maps  initially  assigned  to  and  including  Type  1 . 
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The  counting  procedure  is  then  applied  to  the  remaining  maps  and  the  procedure 
repeated.  This  counting-selection-elimination  procedure  is  repeated  creating 
Types  2,  3,  4,  etc.  until  only  a  few  maps  remain.  It  is  apparent  that  the  object 
of  this  procedure  is  to  identify  map  types  which  are  mutually  exclusive,  i,e., 
where  one  type  does  not  correlate  highly  with  another  type. 

Generally  there  is  a  cut-off  value  preset  below  which  no 
further  types  will  be  chosen.  This  cut-off  value  is  essentially  arbitrary  and  varies 
from  author  to  author.  Lund  (1963)  terminated  selection  of  types  when  he  had 
selected  ten  types.  Hartranft  et  al  (1970a)  terminated  the  selection  and 
assignment  when  the  number  of  maps  in  the  sample  fell  below  two  percent  of  the 
original  sample  or  when  ten  map  types  had  been  selected  -  whichever  occurred 
first.  Later  work  by  the  Map  Typing  Project  Office  in  producing  surface  map 
catalogues  terminated  selection  of  types  when  the  number  of  maps  remaining  to  be 
typed  fell  below  tv/o  percent  of  the  original  sample  or  v/hen  thirty  map  types  had 
been  selected  -  whichever  occurred  first. 

Maps  are  assigned  to  a  selected  weather  type  on  the  basis  of 
the  correlation  coefficient  being  equal  to  or  above  the  threshold  value.  In  the 
initial  assignment,  no  distinction  is  made  concerning  how  much  greater  than  the 
threshold  value  a  map  correlated  with  a  given  type.  It  is  often  Found  that  one 
map  may  correlate  above  the  threshold  value  with  two  or  more  types.  However, 
during  the  initial  analysis,  the  map  is  assigned  to  that  weather  type  v/ith  which 
it  first  correlates  with  a  coefficient  equal  to  or  greater  than  the  threshold  value. 

It  is  then  eliminated  from  eligibility  into  later  created  types. 

Once  the  weather  types  are  chosen  and  all  maps  are  assigned 
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to  a  type  (or  to  an  uncorrelated  "residual"  group),  the  maps  are  recorrelated 
against  the  selected  types  to  determine  which  type  each  map  correlates  with  most 
highly.  This  reassignment  step  is  significant  in  that  it  not  only  assigns  a  map  to 
the  type  it  most  resembles,  but  also  identifies  map  types  which  initially  had  a  low 
frequency  of  occurrence,  and  after  the  reassignment  phase,  become  major  types. 

There  are  minor  variations  to  this  procedure.  If  one  map 
correlated  above  the  threshold  with  two  or  more  map  types,  some  investigators 
would  assign  it  to  both  types  in  an  attempt  to  increase  the  data  base  for 
climatological  studies.  However,  a  problem  arises  in  applying  the  results  of  the 
classification  scheme  to  an  objective  forecast  study.  A  map  may  correlate  above 
the  threshold  value  with  two  weather  types  and  the  pressure  pattern  differences 
may  be  small.  Yet  the  objective  forecast  method  applied  to  the  map  in 
consideration  to  determine  the  probability  of  occurrence  of  some  weather 
phenomena  may  be  entirely  different. 

If  the  number  of  duplicated  classifications  with  more  than  one 
type  are  numerous  it  may  be  advantageous  to  define  "hybrid"  groups.  These 
hybrid  groups  would  increase  the  number  of  types  but  may  be  associated  with 
significantly  different  weather  phenomenon  than  the  two  types  considered 
separately.  The  map  type  assignment  will  depend  upon  the  type  of  forecast  or 
climatological  aid  being  developed. 
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2.3  Applications  of  Objective  Pattern  Classifications 

2.3.1  Case  Studies  of  Special  Weather  Phenomena 

The  primary  purpose  of  this  application  is  to  identify 
recurring  synoptic  patterns  or  map  types  which  occur  before  or  during 
troublesome  and  hazardous  weather  phenomena.  Special  map  type  files  can  be 
prepared  of  surface  and  upper  air  maps  associated  with  these  critical  weather 
situations  and  can  be  used  as  a  valuable  guide  in  diagnosing  the  current  situation, 
interpreting  prognostic  charts  and  alerting  forecasters  to  potentially  dangerous 
synoptic  situations.  This  particular  concept  is  not  new.  The  value  of  preparing 
case  studies  on  critical  weather  phenomena  has  been  publicized  in  meteorological 
literature  for  many  years.  Examples  which  have  been  prepared  for  Alberta  are 
provided  by  Thompson  (1950),  Burns  (1970),  Storr  (1953)  and  Hage  (1957) .  This  approach 
is  commonly  referred  to  as  the  case  study  method  of  investigating  forecast  problems. 

In  developing  a  case  study,  weather  maps  are  collected  in 
one  of  two  ways.  One  method  involves  the  accumulation  of  surface  and  upper 
air  maps  on  an'tos-you-go"  basis,  in  which  maps  and  data  are  collected  and  filed 
each  time  a  particular  weather  event  occurs.  The  second  method  involves  a 
search  of  historical  weather  records  for  days  and  time  periods  in  which  a 
particular  weather  event  occurred.  The  associated  weather  maps  for  these 
occurrences  are  then  obtained  from  historical  weather  map  files. 

Attempts  are  then  made  to  establish  common  relationships 
between  certain  circulation  parameters  and  patterns  and  the  weather  event.  It  is 
difficult  to  identify  objectively  certain  recurring  patterns  or  map  types.  This  is 
an  area  where  computer  map  typing  procedures  play  a  significant  role,  not  only 
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from  the  standpoint  of  saving  a  considerable  amount  of  time,  but  primarily 
because  they  yield  more  objective  and  reliable  results. 

Even  if  certain  map  types  are  identified  subjectively,  the 
problem  that  still  arises  is  that  the  selected  types  suggest  only  the  necessary 
conditions.  Further  testing,  referred  to  as  the  "negative  check,"  is  required  to 
screen  the  entire  weather  map  file  and  select  all  similar  weather  types  which 
conceivably  could  produce  the  weather  phenomenon.  The  negative  check  must  be 
applied  impartially  to  all  situations.  If  the  check  is  valid,  only  a  small  number  of 
similar  map  types  will  be  found  where  the  phenomenon  did  not  occur.  The 
negative  check  will  therefore  reveal  whether  some  additional  screening  must  be 
devised.  Since  this  phase  of  case  studies  is  extremely  tedious  and  time 
consuming,  computerized  map  typing  procedures  can  play  a  significant  role. 

2.3.2  Objective  Map  Typing  for  Case  Studies 

The  general  approach  as  used  by  the  Twelfth  Weather 
Squadron  for  developing  case  studies  is  illustrated  in  Figure  1  .  The  main 
weakness  of  the  method  is  that  a  separate  data  base  must  be  compiled,  and  map 
types  generated,  for  each  different  weather  phenomenon  to  be  investigated. 

The  first  step  is  to  obtain  a  listing  of  the  occurrence  data 
for  the  selected  weather  event.  This  includes  the  data  and  the  period  of  the  day 
for  each  occurrence,  and  any  other  pertinent  information.  The  occurrence  data 
base  should  be  as  large  as  possible  to  ensure  reliable  results. 
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FIGURE  1 

General  Case  Study  Approach  for  a  Special  Weather  Phenomenon 
(from  Hartranft  et  al,  1970a,  p.  18) 
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The  map  typing  procedure  described  in  Section  2.2  is 
applied  to  the  event  occurrence  data  using  the  appropriate  synoptic 
data.  The  choice  of  map  data  depends  on  the  type  of  case  study  being 
developed.  The  pertinent  map  data  are  correlated  for  the  occurrence  dates  in 
order  to  identify  the  map  types  (viz.,  A  ,  B  ,  C  ,  .  .  .).  Those  maps  which 

XXX 

are  not  assigned  to  any  type,  i.e.,  correlated  less  than  or  equal  to  the  threshold 
correlation  coefficient  ,  are  grouped  together  in  'an  uncorrelated'  group  (U  ). 

X 

The  degree  to  which  the  selected  map  types  represent  necessary  conditions  is 
inversely  proportional  to  the  size  of  the  U  group.  If  the  size  of  the  U  group  is 

X  X 

too  large  it  means  that  there  is  too  large  a  variety  of  dissimilar  map  patterns, 
representing  the  weather  event  in  question.  The  circulation  pattern  is  probably 
an  insignificant  causative  factor.  Based  upon  experience  gained  by  the  Twelfth 
We ather  Squadron,  the  size  of  the  U  group  should  not  exceed  ten  percent  of  the 

X 

total  sample. 

Whether  or  not  a  reliable  positive  forecast  can  be  issued 
during  positive  situations  is  determined  only  by  the  results  of  the  negative  check. 
The  negative  check  involves  correlating  all  the  remaining  data  with  each  of  the 
selected  map  types  A  ,  B  ,  C  ,  etc.  Each  non-occurrence  is  assigned  to  a 
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selected  map  type  if  the  correlation  coefficient  exceeds  or  equals  the  threshold 

value.  If  the  non-occurrence  correlates  v/ith  more  than  one  map  type,  it  is 

assigned  to  that  map  type  with  which  it  correlates  me  t  highly.  This  procedure 

will  yield  qroups  of  non-occurrences,  A  ,  B  ,  C  ,.  .  .  U  .  The  A  group, 
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for  example,  would  contain  the  sample  of  non-occurrences  in  which  the  map 
pattern  was  similar  to  the  occurrence  map  type,  A  .  The  relative  frequency  of 
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occurrences  to  non-occurrences  of  an  event  for  a  given  map  type  will  determine 
whether  additional  screening  procedures  are  required.  If  the  relative  frequency 
is  low,  this  is  the  case  and  other  parameters  such  as  moisture,  stability, 
temperature,  or  history  should  be  investigated. 

A  study  of  snowfall  using  this  approach  is  seen  in  Hartranft  et 

al  (1970a,  1970b). 

2.3.3  Area  Map  Type  Catalogues 

A  second  approach  toward  applying  the  weather  map 
classification  scheme  is  to  generate  map  types  using  a  ful.l  data  sample.  A  major 
weakness  of  this  method  when  used  in  the  initial  step  of  a  forecast  study  is  the 
lack  of  control  in  choosing  map  types.  There  will  be  in  all  probability  no 
specific  map  type  which  will  enable  a  reliable  occurrence  or  non-occurrence 
forecast  to  be  made.  As  a  result,  greater  stratification  of  data  is  required.  The 
primary  advantage  is  that  a  catalogue  of  map  types  can  be  generated  for  a 
particular  geographical  area,  and  this  catalogue  used  to  study  a  variety  of 
weather  phenomena.  This  would  eliminate  the  need  to  apply  map  typing 
procedures  to  separate  sets  of  data  when  investigating  different  weather 
phenomena.  Since  the  entire  data  sample  is  used  no  negative  check  is  required. 

The  series  of  steps  which  must  be  followed  to  obtain 
representative  map  types  for  a  particular  area  Is  illustrated  in  Figure  2  and  has 
been  explained  in  Section  2.2. 
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FIGURE  2 

Procedure  for  Determining  Map  Types  for  Area  Map  Type  Catalogue 
(from  Nieman  and  Sabin,  1972,  p.  3) 
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The  typing  procedure  is  straightforward  but  a  problem  arises 
when  the  number  of  maps  becomes  large,  i.e.,  the  excessive  number  of 
correlations  required  initially.  The  correlation  matrix  is  a  symmetric  matrix 
about  the  principal  diagonal  (diagonal  elements  are  all  1  .00).  As  a  result, 
n(n-l)  /  2  correlations  are  required  for  a  data  sample  of  n  maps.  For  one  year 
of  data  (one  map  per  day),  this  entails  approximately  66,430  correlations 
whereas  for  ten  years  of  data,  this  requires  6,666,726  correlations!  Techniques 
which  correlate  the  minimum  amount  of  map  data  are  required,  two  of  which  are 
described  in  the  following  sections. 

2.3.4  Representative  Month  (REMPO)  Technique 

The  REMPO  technique  for  developing  seasonal  catalogues  of 
map  types  is  illustrated  in  Figure  3. 

First  it  is  necessary  to  decide  upon  the  "homogeneous 
seasons"  for  the  particular  region  in  question.  The  seasons  will  not  necessarily 
correspond  to  the  standard  "calendar"  seasons.  The  next  step  is  to  determine  a 
representative  month  (or  similar  period)  within  each  season.  Both  of  the  above 
should  be  determined  from  climatological  data. 
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FIGURE  3 

Representative  Month  Technique  (REMPO)  to  Develop  Seasonal  Map  Type  Catalogue 
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The  actual  procedure  can  be  followed  by  referring  to  Figure 
3.  Data  for  each  day  in  the  representative  period  (several  years  of  data)  are 
correlated  against  all  other  days,  and  map  types  are  developed  as  discussed 
previously.  Complete  data  for  each  month  of  the  season  are  then  correlated 
against  the  types  and  a  type  assigned  to  each  map  time.  If  after  this  step  the 
number  of  uncorrelated  maps  is  high,  it  is  possible  that  the  initial  types  do  not 
adequately  describe  the  "homogeneous"  season.  Therefore,  all  the  uncorrelated 
data  are  assembled  and  correlated  to  develop  additional  types.  These  additional 
types  added  to  the  original  types  provides  map  types  for  the  "homogeneous" 
season . 

2.3.5  Batch  Processing  (BATCH)  Technique 

The  BATCH  technique  for  developing  map  types  is 
illustrated  in  Figure  4.  This  method  may  be  applied  to  develop  types  for  the 
entire  year  or  for  "homogeneous"  seasons  as  with  the  previous  method.  In  the 
latter  case,  since  all  the  data  are  used  for  the  "homogeneous"  seasons,  it  is  not 
necessary  to  determine  a  climatological  representative  month  or  months. 
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FIGURE  4 

Batch  Processing  Method  Used  to  Develop  a  Map  Type  Catalogue 


' 
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Map  types  are  developed  independently  for  each  month. 
Many  of  the  types  from  one  month  to  another  will  be  similar  to  each  other  since 
they  were  developed  independently.  Therefore,  in  order  to  decrease  the  total 
number  of  types,  it  is  necessary  to  combine  types.  This  can  be  done  by  grouping 
the  types  developed  into  homogeneous  seasons  or,considering  all  the  data,  a 
yearly  map  catalogue  can  be  obtained.  Either  procedure  is  done  by  grouping  the 
types  together  and  correlating  against  each  other,  in  developing  the  final  types, 
a  threshold  correlation  value  must  be  applied  which  is  higher  than  or  equal  to 
the  threshold  value  used  to  determine  the  monthly  types.  By  using  a  threshold 
value  of  0.80  or  greater,  this  ensures  that  only  highly  similar  types  are  grouped 
together.  All  uncorrelated  cases  are  considered  to  be  unique  types  and  are 
included  in  the  final  types. 

The  final  phase  is  to  recorrelate  each  map  in  the  data 
sample  and  assign  it  to  a  map  type  on  the  basis  of  the  value  of  the  correlation 
coefficient . 

2.3.6  Comparison  of  REMPO  and  BATCH  Techniques 

Sabin  and  Nieman  (1972)  have  compared  the  two  methods 
used  to  produce  synoptic  map  type  catalogues.  They  found  that  both 
techniques  produce  approximately  the  same  number  of  types.  However,  there  is 
no  exact  one-to-one  correlation  between  a  REMPO  c  d  BATCH  map  type. 
Because  of  the  typing  procedure,  it  would  be  pure  coincidence  if  a  particular 
map  date  was  found  to  represent  a  map  type  in  both  catalogues.  Nevertheless, 
there  are  similarities  between  the  two  sets  which  can  be  seen  in  the  frequency 


. 
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of  occurrences,  the  pressure  regimes  depicted  and  the  actual  correlation  of 
types  in  the  two  sets. 

Since  the  two  methods  produce  comparable  map  type  sets,  the 
choice  of  which  method  to  use  is  primarily  a  matter  of  computer  time  and  storage. 
The  computer  running  time  is  nearly  equal  for  REMPO  and  BATCH.  Since  1972, 
all  Twelfth  Weather  Squadron  map  typing  projects  use  the  BATCH  processing 
method.  The  decision  was  based  on  the  fact  that  BATCH  is  much  more 
straightforward,  requires  less  manual  manipulation  of  data,  and  considers  data  for 
the  entire  period  rather  than  data  for  only  a  representative  month. 

2.4  Disadvantages  of  Correlation  Methods 

Objections  to  map  pattern  classification  using  correlation  methods 
are  sometimes  raised.  Clearly  the  procedure  to  produce  types  are  not  entirely 
free  of  subjectivity.  Factors  which  affect  the  computational  problems  such  as  the 
threshold  correlation  coefficient,  area  size  and  the  number  of  data  points  have 
been  covered  in  previous  sections. 

In  the  assignment  of  days  to  the  catalogue  types,  some 
misclassification  may  occur  because  a  map  may  correlate  above  the  threshold 
value  with  more  than  one  type.  The  usual  procedure,  is  to  place  those  situations 
in  the  group  for  which  the  correlation  coefficient  with  the  'pattern  day1  is  a 
maximum  but  it  may  be  preferable  to  define  'hybrid'  groups,  if  duplicated 
classifications  are  numerous.  The  decision  of  defining  'hybrid'  types  or  not  will 
depend  on  the  weather  phenomenon  to  be  investigated. 

Most  important  of  all,  it  must  be  recognized  that  the  correlation 
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method  relates  each  day  of  a  given  group  to  the  'pattern  day'  and  does  not 
consider  the  pressure  magnitudes  or  correlations  within  members  of  the  group. 

This  clustering  method  may  not  produce  the  optimum  grouping  of  the  data. 

It  may  be  instructive  to  determine  the  mean  pressure  field  for  each 
group,  although  the  features  on  this  are  weaker  than  those  on  the  corresponding 
pattern  day.  The  choice  depends  on  the  time  scale  involved  and  the  type  of  study 
being  undertaken.  It  should  be  stressed  that  in  the  case  study  approach  a  single 
'key'  map  should  be  used  rather  than  a  mean  map  to  represent  the  map  type.  If  a 
mean  map  was  used,  many  days  would  not  correlate  with  the  mean  map.  A  mean 
map  may  be  preferable  if  one  is  looking  for  large-scale  circulation  features  or 
trying  to  minimize  variations  between  maps  within  a  type  to  facilitate  visual 
typing. 

From  a  statistical  point  of  view,  two  major  questions  concern  the 

effects  of  spatial  and  temporal  autocorrelation  and  of  non-linearity.  The  product- 

moment  correlation  coefficient  (equation  1)  theoretically  assumes  that  the  data  used 

in  its  calculation  are  not  serially  correlated,  i.e.,  there  is  no  correlation  between 

X.,  ,  and  X.  .  ,  ,  (k,  m  ~  1 ,  2, . . .).  T  he  fact  that  one  is  spatially  and 

i  t  i  +  k,  t  +  m 

temporally  sampling  a  continuous  field  variable  (pressure  or  height)  puts  the  validity 
of  this  assumption  in  jeopardy.  Caution  must,  therefore,  be  used  in  interpreting  such 
correlations,  certainly  in  terms  of  their  significance  level.  The  second  difficulty 
is  that  fields  of  pressure  or  geopotential  height  are  generally  non-linear  and 
considerable  information  may  therefore  be  lost  in  correlating  fields  on  the  basis  of 
a  linear  technique. 
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On  the  whole,  the  advantages  of  using  linear  correlation  methods  far 
outweigh  the  disadvantages.  The  development  of  synoptic  climatological  aids  based 
upon  map  types  allows  an  integration  of  local  area  climatology  into  each  forecast. 


CHAPTER  III 


SURFACE  MAP  CATALOGUE  FOR  ALBERTA 

3.1  Introduction 

The  Map  Typing  Project  Office  (Twelfth  Weather  Squadron,  United 
States  Air  Force,  Colorado  Springs)  has,  since  1968,  been  conducting  applied 
research  into  the  field  of  computer  derived  weather  map  types.  In  1972,  this 
office  began  the  production  of  surface  pressure  map  catalogues  for  fourteen 
regions  (windows)  in  the  United  States.  A  surface  map  catalogue  has  been 
produced  by  the  above  office  for  Alberta  which  v/ill  be  described  in  detail  in  the 
following  sections. 

3.2  Data  and  Area  Used  for  the  Alberta  Window 

The  data  used  in  the  preparation  of  the  map  type  catalogues  were 
obtained  from  the  Fleet  Numerical  Weather  Central  (FNWC),  Monterey, 
California  and  consists  of  grid  point  mean  sea  level  surface  pressure  data 
(standard  grid  spacing  of  381  kilometers)  for  an  eighteen  by  seventeen  point  grid 
covering  the  United  States  and  most  of  Canada.  The  total  period  of  record  was 
January  to  December  1946  -  1971 .  The  pressure  data  at  the  12Z  (Greenwich 
Mean  Time)  reporting  time  are  available  for  each  grid  point  for  the  entire  period 
while  values  for  the  00Z  times  are  available  only  from  1955-1959  and  1962-1971 . 
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The  area  covered  by  the  Alberta  surface  map  catalogue  is 
illustrated  in  Figure  6,  Appendix  1  .  Ideally,  the  area  in  which  one  is 
interested  should  be  centered  in  the  window.  This  was  not  possible  for  the 
Alberta  window  because  of  data  availability.  Since  most  weather  systems  which 
affect  Alberta  originate  in  the  Pacific,  it  was  decided  to  include  more  data 
grid  points  in  the  west  rather  than  the  east.  The  uppermost  boundary  of  available 
grid  point  data  coincides  with  the  upper  boundary  of  the  window  (grid  points  one, 
two,  three).  The  computer  program  used  by  the  Map  Typing  Project  Office 
allows  a  maximum  of  twenty-four  grid  points  in  the  production  of  a  catalogue. 
Once  derived,  the  map  types  apply  to  all  locations  which  fall  within  the  window 
but  it  has  been  found  that  the  similarity  of  maps  within  a  type  generally  extends 
far  beyond  the  area  used  to  determine  the  type  for  windows  of  the  size  used. 

3.3  Analysis  Procedure 

The  analysis  procedure  follows  the  BATCH  processing  method 
(section  2.3.5)  with  some  modifications  because  of  the  large  amount  of  data. 

The  initial  phase  in  the  selection  of  map  types  involves  the  computation  of 
correlation  coefficients  of  each  map  with  all  other  maps  for  a  given  month.  The 
initial  correlations  are  carried  out  using  12Z  map  data  for  the  period  1954-1971  . 
As  the  coefficients  are  determined  the  number  of  correlations  that  equal  or 
exceed  0.70,  0.80  and  0.90  are  tabulated. 

The  next  phase  Involves  the  determination  of  map  types  for  each 
month  in  the  following  manner  (for  a  more  detailed  description,  see  Chapter  I  I): 
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(  i)  The  map  that  correlates  with  the  greatest  number  of  other  maps  for  a  given 
month,  with  a  correlation  coefficient  of  0.80  or  higher,  i  .e. ,  threshold 
value  0.80  ,  is  selected  as  a  map  type.  If  two  or  more  maps  have  the  same 
number  of  correlations  equal  to  or  above  the  threshold  value,  the  most 
recent  map  is  selected  as  the  type. 

(  ii)  All  maps  that  correlate  with  the  selected  key  map  with  a  coefficient  of 
0.80  or  higher  are  eliminated  from  further  consideration. 

(iii)  Steps  ( i)  and  (ii)  are  repeated  until  either  the  number  of  types  reaches  a 
predetermined  maximum  (usually  30),  or  the  number  assigned  to  a  type  is 
less  than  a  predetermined  minimum  (usually  six).  . 

(  iv)  Once  the  map  types  are  selected,  each  map  for  that  month  is  assigned  to 
the  map  type  with  which  its  correlation  is  the  highest.  However,  if  its 
highest  correlation  does  not  equal  or  exceed  0.80,  no  assignment  is  made. 

(  v)  Any  map  type  having  fewer  than  two  percent  of  the  total  number  of  maps 
for  that  month  assigned,  is  eliminated  as  a  map  type.  The  dates  and  grid 
point  pressure  values  for  each  type  (key  map)  are  stored  for  further 
processing . 

The  data  are  processed  one  month  at  a  time,  i.e.,  all  maps  occurring  in  January 
for  the  period  of  1954-1971  are  processed  first,  then  February,  March,  etc. 

Two  reasons  for  processing  the  data  in  this  manner  are  computer  core  availability 
and  data  storage  in  monthly  increments. 

Once  the  map  types  are  determined  for  each  of  the  twelve  months, 
the  correlation  of  each  monthly  map  type  with  all  other  monthly  types  is  done. 
This  determines  the  final  map  types  for  the  area  by  consolidating  similar  types. 
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The  final  types  are  determined  in  the  following  manner: 

(  i)  The  monthly  map  type  that  correlates  with  the  greatest  number  of  other 
map  types  (  >  0.80)  is  selected  as  a  final  map  type.  Again,  the  most 
current  monthly  type  is  selected  if  two  or  more  monthly  types  have  the  same 
number  of  correlations. 

(  ii)  All  monthly  map  types  that  correlate  with  the  selected  type  with  a 

coefficient  of  0.80  or  higher  are  eliminated  from  further  consideration. 

(iii)  Steps  (i)  and  (ii)  ore  repeated  unless  the  number  of  final  types  reaches  a 
predetermined  maximum  (usually  40). 

(  iv)  All  monthly  map  types  that  have  not  been  assigned  to  a  final  type  are 

considered  unique  and  are  designated  as  final  map  types  unless  the  number 
of  types  exceeds  forty. 

(  v)  The  above  procedure  yields  a  map  catalogue  of  recurring  synoptic  patterns 
for  the  entire  year.  Once  the  final  types  are  determined,  each  map  in  the 
data  period  (1946-1971)  is  assigned  to  the  map  type  with  which  it 
correlates  the  highest  (primary  type),  provided  the  correlation  coefficient 
is  greater  than  or  equal  to  0.80.  If  this  is  not  the  case,  the  map  is 
assigned  to  the  uncorrelated  (U)  group. 

(  vi)  The  second  highest  correlation  coefficient  is  also  computed  for  each  map, 
and  if  it  exceeds  the  threshold  value,  the  map  type  associated  with  that 
coefficient  is  recorded  (secondary  map  type),  n  addition,  if  the  second 
highest  correlation  coefficient  is  within  0.05  of  the  first,  as  well  as  being 
equal  to  or  greater  than  0.80  ,  the  occurrence  of  a  secondary  highly 
similar  type  is  noted . 
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3.4  Further  Processing  by  Colorado  Springs 

Once  the  assignment  is  complete,  a  tabulation  of  map  type 
occurrences  and  frequency,  information  on  types  most  frequently  following  a  given 
type,  and  information  on  persistence  of  each  type  is  performed.  The 
characteristics  of  the  types  are  stratified  for  each  month,  providing  very  useful 
information  on  the  behaviour  of  synoptic  systems. 

3.5  Results  for  the  Alberta  Surface  Map  Catalogue 

In  the  production  of  the  map  types  for  Alberta,  a  total  of  6,579 
mapswas  used  in  the  initial  correlations  and  14,554  maps  were  assigned  to  the 
final  map  types.  The  catalogue,  run  on  a  mini-computer  system,  required 
approximately  ten  hours  of  computer  processing  time,  which  Is  about  equally 
divided  betv/een  the  central  processing  procedure  and  the  input/output  section. 

3.5.1  Monthly  Surface  Types 

Table  1  illustrates  the  results  obtained  for  each  month  when 
the  12Z  maps  were  processed  for  the  period  1954-1971  .  Due  to  an  error  in  a 
data  date,  1953  data  for  July  were  also  processed  giving  585  days  (maximum  was 
to  be  558).  This  data  error  has  no  effect  on  the  final  generation  of  the  map 
types  since  no  key  map  was  selected  from  the  1953  data.  It  is  interesting  to  note 
that  the  average  number  of  map  types  for  each  month  is  eighteen  and  that  this 
remains  relatively  constant  throughout  the  year.  The  number  of  uncorrelated 
maps  averages  twenty-four  percent  but  a  noticeable  difference  can  be  observed. 
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TABLE  1 

Surface  Map  Types  for  Alberta  for  Each  Month  and  the  Final  Yearly  Totals 


MONTH 

NUMBER  OF 
DAYS 

NUMBER  OF 

MAP 

TYPES 

NUMBER  OF 

UNCORRELATED 

MAPS 

(MONTHLY 

PERCENTAGE) 

January 

558 

19 

126  (23) 

February 

508 

17 

110  (22) 

March 

554 

15 

143  (27) 

April 

538 

17 

147  (27) 

May 

557 

21 

137  (25) 

June 

538 

19 

157  (29) 

July 

585 

18 

161  (28) 

August 

558 

20 

142  (25) 

September 

540 

18 

149  (28) 

October 

548 

16 

95  (17) 

November 

539 

19 

103  (19) 

December 

556 

18 

122  (22) 

TOTALS 

1954-1971 

6,579 

217 

1 , 597  (24) 

FINAL 

CATALOGUE 

1946-1971 

14,554 

33 

3,836  (26) 

H 

. 
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The  maximum  number  of  uncorrelated  da/s  occurs  during  the  summer  period,  with 
a  minimum  during  the  October  -  December  period.  This  is  probably  the  result  of 
stronger  thermal  gradients  in  the  late  fall  and  early  winter,  compared  with  the 
summer.  The  same  results  are  evident  for  the  uncorrelated  group.  Appendix  2. 

3.5.2  Yearly  Surface  Map  Types 

The  map  types  for  each  month  (total  217)  have  been 
consolidated  into  a  catalogue  for  the  year  comprising  thirty-three  map  types. 
These  types  account  for  seventy-four  percent  of  the  days  in  the  development 
sample  when  each  map  time  is  correlated  with  the  key  pattern  of  the  thirty-three 
types. 

In  the  development  of  the  surface  types.  Map  Type  1 
was  the  most  predominant  recurring  pattern  of  the  monthly  types;  Map  Type 
2,  the  second  and  so  on. 

However,  when  the  reassignment  of  each  date  in  the  period 
of  record  was  completed.  Map  Type  4  became  representative  of  the  most  days, 
i.e.,  the  most  predominant  excluding  the  uncorrelated  group  ,  then  Map  Type 
2,  etc. 

3.5.3  Display  of  Surface  Types  and  Characteristics 

Appendix  2  contains  the  characte-istics  of  each  map  type  and 
the  associated  key  map  pattern.  The  map  type  most  frequently  following  and  the 
most  frequent  duration  statistics  are  computed  for  the  types  when  data  dates  are 
separated  by  a  twelve-hour  period  (recall  that  all  00Z  maps  were  not  available). 
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In  all  cases,  the  statistics  are  stratified  by  month  and  show  the  seasonal 
characteristics  in  the  map  types.  For  example.  Map  Type  1  is  predominantly  a 
winter  phenomenon,  occurring  infrequently  in  the  summer.  This  type  is 
characteristic  of  a  cold  outbreak  of  Arctic  air.  The  opposite  is  true  for 
(predominantly  a  summer  and  fall  map).  This  type  results  because  the 
strongest  thermal  gradient  is  shifted  northward  at  this  time  of  the  year. 

Other  characteristics  can  best  be  explained  by  referring  to 
Map  Type  1,  Appendix  2.  The  statistics  on  the  map  type  most  frequently 
following  in  the  next  twelve  hours  is  tabulated  for  all  map  types  but  only  the  five 
most  frequent  are  listed.  For  Map  Type  1  in  January,  Map  Type  1  will  follow 
itself  in  the  next  twelve  hours,  forty-seven  percent  of  the  time.  Other  map 
types  which  have  been  found  to  follow  are:  Map  Types  27  (16%),  16(8%), 
uncorrelated  U  (8%)  and  30  (5%). 

In  the  persistence  statistics,  all  consecutive  twelve-hour 
occurrences  of  the  same  map  type  are  tabulated  and  frequencies  computed.  The 
four  most  frequent  durations  are  listed  under  this  heading.  Referring  to  Map  Type 
1  in  January,  fifty-three  percent  of  the  maps  processed  occurred  for  a  single 
twelve-hour  period,  twenty-three  percent  for  two  twelve-hour  periods,  thirteen 
percent  for  four  and  nine  percent  for  three.  The  map  type  statistics  for 
following  types  and  persistence  show  the  seasonal  migration  of  pressure  systems. 

3.6  H igh I y  Similar  Surface  Map  Types 

The  selection  of  map  types  was  carried  out  such  that  no  two  types 
correlate  above  or  equal  to  the  threshold  correlation  value  with  one  another. 
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Intercorrelation  of  the  surface  map  types  have  confirmed  this  hypothesis.  This 
procedure  also  provides  a  means  of  checking  the  selection  of  the  map  types. 

Map  T ype  2  and  17  illustrate  the  problems  in  visual 
typing.  The  correlation  between  these  two  types  is  below  the  threshold  value  of 
0.80  but  when  the  correlation  is  output  by  the  computer,  round-off  produces  a 
value  of  0.80.  Inspection  of  the  key  patterns  show  that  the  maps  are  indeed 
similar,  the  difference  being  related  to  the  intensity  of  the  system  and  the 
variation  in  one  area  of  the  window.  The  flow  pattern  over  Alberta,  however, 
is  definitely  similar. 

In  order  to  obtain  more  information  on  these  two  map  types,  a 
comparison  of  the  primary  and  secondary  correlations  for  each  date  in  the  period 
May  to  September  1963-1971  was  made.  This  time  period  was  chosen  since  the 
500  mb.  map  catalogue  and  the  precipitation  probabilities  described  in  other 
chapters  used  the  same  time  period.  The  comparison  indicated  that  the 
pattern  occurred  more  frequently  than  Map  Type  17.  It  was  also  found 

that  when  the  primary  map  type  for  a  date  was  Type  2,  the  secondary  was 
frequently  Type  17.  The  reverse,  however,  never  occurred,  i.e.,  primary 

Type  17,  secondary  Type  2. 

3.7  Highly  Uncorrelated  Maps 

Another  aspect  of  the  intercorrelations '  'hich  were  computed,  was 

the  highly  uncorrelated  maps,  i.e.,  more  specifically  the  correlations  less  than 

or  equal  to  -0.80.  This  implies  that  the  maps  are  similar  but  have  reversed  flow 

directions.  Table  2  indicates  the  map  types  which  correlate  with  each  other  less 
than  or  equal  to  the  negative  of  the  threshold  value. 
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17 
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TABLE  2 

lntercorrelation  of  Surface  Map  Types. 
Coefficients  Less  Than  or  Equal  -0.800. 


CORRELATION 
MAP  TYPE  COEFFICIENT 


2 

-  0.956 

9 

-  0.846 

17 

-  0.831 

10 

-  0.827 

6 

-  0.822 

30 

-  0.847 

5 

-  0.844 

27 

-  0.863 

12 

-  0.871 

15 

-  0.828 

10 

-  0.895 

14 

-  0.837 

23 

-  0.897 

33 

-  0.814 

39 


Visual  inspection  of  the  key  pattern  maps  show  that  the  flow 
directions  are  completely  reversed  for  those  types  indicated  in  Table  2.  The  best 
example  is  Map  Types  1  and  2,  which  show  that  the  locations  of  the  high  and 
low  pressure  areas  are  essentially  interchanged. 
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CHAPTER  IV 


SUMMER  500  MB.  MAP  CATALOGUE 


4.1  Introduction 

A  summer  500  mb.  map  catalogue  of  recurring  map  patterns  was 
produced.  This  catalogue,  when  used  in  conjunction  with  the  surface  map 
catalogue,  will  reveal  information  concerning  the  three-dimensional  nature  of 
synoptic  pressure  systems.  The  following  sections  describe  the  data,  areal  size  and 
the  analysis  procedure  used  for  the  generation  of  this  catalogue. 

4.2  500  mb.  Summer  Map  Catalogue 
4.2.1  Data  and  Areal  Size 

Height  data  for  500  mb.  v/as  abstracted  for  twenty 
rawinsonde  stations  in  Western  Canada  and  the  Northwestern  United  States  using 
Monthly  Upper  Air  Summaries  and  analyzed  maps  from  the  Edmonton  Weather 
Office.  The  data  period  covered  was  May  to  September  1963-1971  (2,754  maps, 
00Z  and  12Z).  Figure  7,  Appendix  3,  illustrates  the  stations  used  in  the  final 
catalogue. 

Initial  data  abstraction  included  ship  4YP  but  after  some 
preliminary  computations  were  rejected  since  great  variation  could  occur 
in  maps  between  4YP  and  the  coast  of  British  Columbia.  A  large  area  was 
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purposely  chosen  since  the  location  of  certain  weather  systems  are  known  to 
change  dramatically  in  a  twelve-hour  period  (Hage,  1957). 

In  order  to  obtain  a  complete  data  base  for  the  period 
chosen,  values  were  interpolated  from  analyzed  charts  if  they  were  missing  for 
any  reason.  This  was  also  the  case  for  stations  which  had  been  relocated  at  some 
time  during  the  period  of  record,  i.e.,  values  were  interpolated  at  the  present 
(1973)  locations  (UIL,  XD).  The  Vernon  rawinsonde  (VK)  was  not  operational 
until  the  fall  of  1971  so  values  were  interpolated  from  analyzed  charts  for  the 
entire  period  1963-1971 .  This  makes  the  catalogue  compatible  with  the  present 
day  reporting  locations. 

The  areal  extent  of  the  500  mb.  catalogue  is 
approximately  twenty  degrees  latitude  by  thirty  degrees  longitude.  This  is 
double  the  recommended  size  of  windows  (Hartranft,  1970a)  and  is  slightly  larger 
than  the  area  used  by  most  authors  except  Augulis.  This  large  area  was  deemed 
necessary  due  to  the  scarcity  of  station  data  and  the  unavailability  of  grid  point 
data.  A  larger  area  also  seems  necessary  to  adequately  describe  the  flow 
pattern  over  Alberta. 

4.2.2  Analysis  Procedures  and  Results 

The  map  types  were  generated  using  the  REMPO 
technique,  and  the  final  results  are  contained  in  Appendix  4. 

June  and  July  were  chosen  as  the  representative  period  to  be 
used  initially.  This  period  was  chosen  because  it  coincides  with  the  maximum 
mean  monthly  precipitation  for  most  Alberta  stations.  Map  types  obtained  during 
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this  period  would  hopefully  be  useful  in  studying  the  occurrence  of 
precipitation.  Only  the  12Z  maps  were  analyzed  in  order  to  decrease  the 
computations. 

The  Map  Typing  Project  Office  has  found  that  the  optimum 
threshold  correlation  coefficient  for  upper  air  maps  is  0.90.  This  threshold  was 
applied  to  the  data  sample  but  it  became  evident  that  this  value  was  too  large, 
probably  a  result  of  the  large  area  used.  The  threshold  value  was  decreased  to 
0.85  for  the  final  map  types  for  the  catalogue.  Table  3  shows  the  comparison  of 
the  results  for  each  threshold  value.  In  the  selection  of  map  types,  a  minimum 
cut-off  value  of  one-half  percent  of  the  data  sample  was  used  below  which  a  map 
would  not  be  considered  a  type.  With  a  threshold  value  of  0.90,  twenty-seven 
types  were  initially  selected  and  after  assignment  of  all  data  in  every  month, 
thirty-five  percent  remained  uncorrelated .  A  threshold  value  of  0.85  produced 
thirteen  types  and  twenty-four  percent  remained  uncorrelated.  The  uncorrelated 
groups  for  threshold  correlation  coefficients  of  0.90  and  0.85  were  considered  too 
large.  The  uncorrelated  groups  were  separately  correlated  and  additional  types 
selected  (0.90  and  0.85  threshold  values  with  selection  terminated  at  one-half 
percent).  The  higher  threshold  (0.90  produced  forty-six  additional  types  for  a  total 
of  seventy-three.  The  0.85  threshold  produced  an  additional  twenty-two  types  for  a 
total  of  thirty-five.  The  uncorrelated  group  had  dropped  to  eighteen  and  twelve 
percent  respectively.  Although  the  use  of  0.90  produces  more  distinct  types,  the 
number  of  them  appears  to  be  excessive.  The  results  obtcined  for  a  threshold 
value  of  0.85  were  used  as  the  map  types  for  the  final  summer  500  mb.  catalogue. 
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TABLE  3 

Map  Typing  Results  for  500  mb.  using  0.90  and 
0.85  Threshold  Correlation  Coefficients 


Data:  500  mb.  heights  at  nineteen  stations.  Period  of  record.  May  - 

September  1963-1971  (2,754  maps). 

RepresentaKve  Period: 

June  and  July  1963-1971  .  12Z  maps  only. 


Threshold  Coefficient  0.90  0.85 

Number  of  Initial  Types 

(cut-off  less  than  3  maps)  27  13 

Uncorrelated  Maps  After 

Assignment  956  (35%)  649  (24%) 

Additional  Types  From 

Uncorrelated  Maps  46  22 

Total  Number  of  Map 

Types  73  35 

Total  Number  of 
Uncorrelated  Maps  After 

Final  Reassignment  482  (18%)  322  (12%) 
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It  is  interesting  to  note  that  the  date  chosen  for  the  key  map 
for  Map  Type  1  in  the  final  catalogue  would  have  been  the  same  if  the 
threshold  value  were  0.90.  The  number  of  maps  occurring  for  each  of  the  final 
types  in  the  initial  selection  and  the  final  assignment  is  shown  in  Table  4.  The 
first  thirteen  types  were  developed  from  the  representative  period  (June  and  July, 
549  12Z  maps)  while  the  remainder  were  selected  from  the  uncorrelated  group  of 
649  maps.  The  importance  of  the  final  reassignment  phase  of  the  typing 
procedure,  i.e.,  map  assigned  to  the  type  with  which  it  correlates  most  highly, 
is  shown  in  Table  4.  Map  Type  8  was  the  eighth  pattern  to  be  selected 
initially  but  after  the  final  reassignment  of  the  data  base  for  all  months,  this  type 
becomes  the  fourth  most  predominant  (excluding  the  uncorrelated  group). 

In  the  production  of  the  surface  map  catalogue,  a  map  type 
which  had  less  than  two  percent  of  the  development  sample  assigned  to  it  was  not 
regarded  as  a  map  type.  This  procedure  was  disregarded  for  the  500  mb. 
catalogue  and  all  types  selected  were  retained  no  matter  what  the  outcome  of  the 
reassignment  phase.  If  a  one  percent  cut-off  value  had  been  used,  one-half  the 
types  would  have  been  rejected,  whereas  a  one-half  percent  cut-off  would  have 
resulted  in  the  loss  of  three  types,  namely  29,  31,  and  34. 
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TABLE  4 

Number  of  Maps  Associated  With  Each  500  mb.  Map  Type 


Map  Type 

Number  of  Maps 
Initially  Assigned 

Number  of  Maps 

Final  Assignment 

Percent  of 
Total  Data 

1 

199 

525 

19 

2 

58 

385 

14 

3 

54 

249 

9 

4 

32 

126 

5 

5 

21 

58 

2 

6 

16 

65 

2 

7 

11 

40 

1 

8 

9 

187 

7 

9 

8 

19 

1 

10 

8 

27 

1 

11 

7 

75 

3 

12 

6 

28 

1 

13 

6 

15 

1 

14 

38 

59 

2 

15 

32 

66 

2 

16 

28 

51 

2 

17 

23 

24 

1 

18 

22 

35 

1 

19 

19 

55 

2 

20 

18 

18 

1 

21 

16 

27 

1 

22 

16 

26 

1 

23 

13 

17 

1 

24 

13 

15 

1 

25 

9 

64 

2 

26 

9 

20 

1 

27 

9 

17 

1 

28 

9 

36 

1 

29 

9 

13 

0 

30 

8 

18 

1 

31 

8 

13 

0 

32 

7 

15 

1 

33 

7 

19 

1 

34 

7 

7 

0 

35 

7 

18 

1 

U 

322 

12 

. 
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Characteristics  of  the  map  types  were  computed  for  each  map 
type  as  for  the  surface  catalogue.  All  statistics  were  tabulated  by  month.  The 
statistics  illustrate  the  various  motions  which  can  occur  for  a  variety  of  map 
patterns.  No  extensive  description  can  be  attempted  because  of  the  large  amount 
of  data  obtained.  The  next  chapter  when  precipitation  and  the  combination  of 
surface  and  500  mb.  map  types  are  considered,  will  cover  some  map  type 
characteristics.  It  is  worthwhile  to  note  that  visual  inspection  of  the  key  map 
patterns  and  characteristics  indicate  that  the  area  chosen  and  modifications  to  the 
analysis  procedure  has  not  produced  undesirable  results.  The  area  used  for  the 
analysis  was  purposely  made  relatively  large  based  on  the  assumption  that  weather 
patterns  may  change  drastically  over  a  relatively  short  period  of  time.  This  is 
evident  from  the  statistics  on  the  following  types.  Map  Type  3,  which  has  a 
closed-low  circulation  centered  on  the  coast  of  British  Columbia  occasionally 
changes  in  twelve  hours  to  a  Map  Type  15,  which  has  a  closed  circulation  in 
Northern  Alberta. 

Retention  of  all  map  types  obtained,  regardless  of  the 
outcome  of  the  reassignment  phase,  may  enable  the  catalogue  to  be  used  for 
seasons  other  than  the  one  for  which  it  was  developed.  The  weather  patterns 
which  occur  in  other  seasons  should  bear  a  resemblance  to  the  map  types  already 
obtained.  The  differences  will  probably  be  evident  only  in  the  map  type 


characteristics. 
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CHAPTER  V 


DEVELOPMENT  OF  SYNOPTIC  CLIMATOLOGICAL  AIDS 

5.1  Climatological  Aids  From  Map  Types 

Given  a  catalogue  of  map  types,  a  variety  of  synoptic 
climatological  aids  can  be  developed  for  each  map  type  using  weather  records 
from  individual  stations.  Any  weather  phenomena  for  which  the  synoptic  flow 
pattern  is  the  dominant  factor  sjch  as  precipitation,  visibility  and  ceilings 
could  be  investigated.  Since  the  catalogue  contains  statistical  information  on 
the  type  occurrence,  the  following  map  type  and  type  duration,  it  is,  in  itself, 
a  climatological  aid  providing  useful  information  on  synoptic  patterns  occurring 
within  the  area. 

Conditional  climatological  tables  of  ceiling  and  visibility 
(commonly  referred  to  as  "persistence  probability"  or  PP  tables)  have  been  in  use 
for  many  years  by  most  weather  detachments  in  the  United  States  Air  Force  and 
Army.  A  weakness  of  these  PP  tables  is  their  inclusion  of  ceiling  and  visibility 
conditions  which  are  due  to  many  different  synoptic  and  smaller  scale 
situations.  PP  tables  of  weather  variables  based  upon  a  factor  such  as  wind 
direction  or  surface  pressure  could  likely  be  improved  if  they  were  based  on  map 
types. 

Hartranft  et  al  (1970a,  1970b)  have  used  map-typing  in  predicting 
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heavy  snowfall  (greater  than  or  equal  to  two  inches)  for  Colorado 
Springs.  By  using  a  combination  of  surface  and  700  mb.  map  types,  reliable 
"no"  forecasts  were  obtained, but  in  order  to  separate  the  snow  producing  storms 
v/ithin  each  map  type,  additional  screening  and  parameters  such  as  moisture 
(surface  and  aloft)  and  pressure  gradients  were  needed.  Lund  (1963)  found  that 
useful  information  could  be  obtained  about  probable  weather  by  knowing  the  map 
type,  i.e.,  precipitation  amount,  snowfall  amount  and  percent  of  possible 
sunshine  at  Boston,  Massachusetts. 

The  major  advantage  of  synoptic  climatological  aids  based  on 
computer-derived  map  types  is  that  the  synoptic  flow  pattern  is  taken  into 
consideration  and  prognostic  chart  information  can  be  incorporated  for 
operational  use. 

5.2  Comparison  of  500  mb.  and  Surface  Map  Types 

Although  the  500  mb.  and  surface  map  types  have  no  one-to-one 
relationship,  an  examination  of  surface  patterns  associated  with  each  500  mb. 
type  will  reveal  information  concerning  the  intensity  of  systems  and  mean 
temperature  patterns.  No  quantitative  mean  thermal  patterns  can  be  obtained 
since  the  data  values  for  each  catalogue  were  not  identical  (number  of  data 
points  and  spacing);  however,  qualitative  thermal  patterns  can  be  inferred  for 
well  developed  situations. 

After  map  types  for  500  mb.  and  surface  were  developed 
separately,  contingency  tables  of  their  mutual  occurrence  were  prepared  for  the 
period  May  1  to  September  30,  1963-1971  .  In  all,  1,224  combinations  are 
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possible  with  the  number  of  types  obtained.  Considering  only  primary  (i.e., 
highest  correlation)  500  mb.  and  surface  types,  it  was  found  that  814  (66.5%) 
combinations  did  not  occur,  whereas  when  both  primary  and  secondary  (second 
highest  correlation  coefficient)  types  were  considered,  628  (51 .3%)  did  not  occur. 
One  hundred  and  seventy-seven  combinations  of  primary  500  mb.  and  surface  map 
types  (14.5%)  accounted  for  seventy-seven  percent  of  the  days  in  the  nine-year 
period.  Every  surface  type  was  possible  with  only  one  upper  air  type,  namely 
number  one.  This  could  be  expected  since  this  type  occurs  most  frequently  and 
numerous  thermal  patterns  can  be  envisioned  depending  on  the  intensity  of  the  upper 
air  flow. 

The  five  most  frequent  primary  surface  types  associated  with  each  500 
mb.  type  are  shown  in  Table  5,  The  most  frequent  surface  type  associated  with  an 
upper  air  pattern  is  the  uncorrelated  group,  a  result  not  surprising  since  this  group 
occurs  most  frequently  in  the  summer  and  is  a  result  of  poorly  organized  weather 
systems.  The  strong  similarity  between  surface  Types  2  and  17  has  previously  been 
mentioned  in  section  3.4.  Map  Type  17  occurs  more  frequently  than  Type  2  with 
the  500  mb.  Types  6,  18,  20,  and  29,  each  characterized  by  a  ridge  west  of 
Alberta.  An  inspection  of  upper  air  patterns  which  are  associated  with  Map  Type 
2,  indicates  that  they  all  are  characterized  by  a  trough  west  of  Alberta. 
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TABLE  5 

Predominant  Surface  Map  Types  and  Frequency  Associated  With  Each 
500  mb.  Map  Type  (May  -  September,  1963-1971) 


Primary  500  Number  of 

mb. Map  Type  Maps  Primary  Surface  Map  Types  (Frequency  %) 


1 

525 

U  (32) 

4(16) 

2(15) 

9(  7) 

13(4) 

2 

385 

U(38) 

2(15) 

4(12) 

9(  6) 

7(4) 

3 

249 

U(34) 

4(11) 

2(10) 

9(  6) 

14(5) 

4 

126 

U(45) 

4(17) 

6(  9) 

7(  4) 

13(4) 

5 

58 

U(38) 

9(16) 

2(12) 

25  (  9) 

14(7) 

6 

65 

U(52) 

17(  9) 

27  (  9) 

1(  6) 

2(5) 

7 

40 

4(35) 

U(28) 

14(10) 

10(  8) 

25(8) 

8 

187 

U(34) 

9(14) 

12(14) 

3(11) 

2(6) 

9 

19 

U(68) 

1(11) 

19(11) 

9(  5) 

10(5) 

10 

27 

U(63) 

27(15) 

23(.  7) 

K  4) 

2(4) 

11 

75 

U(28) 

4(25) 

2(12) 

3(  5) 

9(4) 

12 

28 

1(32) 

U(18) 

4(14) 

19(11) 

6(7) 

13 

15 

U(57) 

25(27) 

2(  7) 

19(  7) 

26(7) 

14 

59 

U(34) 

4(17) 

6(  8) 

7(  8) 

9(8) 

15 

66 

4(53) 

U(12) 

14(  8) 

2(  5) 

6(5) 

16 

51 

U(22) 

9(18) 

2(12) 

28(10) 

3(8) 

17 

24 

U(46) 

1(13) 

6(  8) 

19(  8) 

2(4) 

18 

35 

U(40) 

10(11) 

27(11) 

K  9) 

17(6) 

19 

55 

U(49) 

2(11) 

10(  9) 

5(  7) 

26(7) 

20 

18 

U(39) 

9(17) 

17(17) 

2(11) 

12(6) 

21 

27 

U(44) 

4(26) 

25(11) 

H  7) 

9(4) 

22 

26 

U(54) 

12(12) 

5(  8) 

19(  8) 

1(4) 

23 

17 

U(41) 

9(18) 

7(12) 

27(12) 

4(6) 

24 

15 

U(60) 

19(13) 

K  7) 

6(  7) 

12(7) 

25 

64 

2(25) 

4(20) 

U(19) 

13(17) 

6(5) 

26 

20 

U(55) 

23(15) 

1(  5) 

2(  5) 

11(5) 

27 

17 

4(41) 

U(18) 

6(12) 

14(12) 

9(6) 

28 

36 

U(47) 

27(19) 

10(  8) 

12(  8) 

2(3) 

29 

13 

U(54) 

4(  8) 

10(  8) 

12(  8) 

17(8) 

30 

18 

4(56) 

9(17) 

’(  6) 

14(  6) 

25(6) 

31 

13 

U(62) 

4(  8) 

7(  8) 

9(  8) 

23(8) 

32 

15 

U(40) 

30(27) 

1(13) 

19(  7) 

26(7) 

33 

19 

4(58) 

U(21) 

13(11) 

19(  5) 

25(5) 

34 

7 

U  (71 ) 

20(14) 

30(14) 

35 

18 

4(28) 

U(28) 

26(11) 

6(  6) 

7(6) 

U 

322 

U(41) 

4(  9) 

6(  6) 

!4(  6) 

19(6) 
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By  considering  the  upper  air  and  surface  patterns,  the  intensity  and 
mean  temperature  patterns  of  systems  could  be  inferred.  The  large  amount  of 
information  precludes  any  complete  discussion  of  results.  The  evolutionary  aspect 
of  each  paired  type  for  the  day  of  occurrence  and  succeeding  time  periods  could 
be  important  in  connection  with  forecasting  applications. 

5.3  Probabilities  of  Precipitation  From  Map  Types 

5.3.1  Precipitation  Data  and  Analysis  Procedures 

Precipitation  data  were  obtained  from  Monthly  Records 
published  by  the  Atmospheric  Environment  Service  (AES)  for  twenty-three  first 
order  observing  stations  in  Northeastern  British  Columbia  and  Alberta  for  the 
period  May  to  September  1963-1971  (1,377  days).  The  station  locations  and 
identifiers  used  are  illustrated  in  Figure  5.  Since  the  map  type  catalogues  were 
obtained  by  using  data  at  00Z  and  12Z,  precipitation  amounts  should 
preferrably  be  centered  at  these  times.  This  would  allow  both  map  type  times  to 
be  used  in  developing  precipitation  frequency  (hereafter  referred  to  as 
precipitation  probabilities  or  PoPs).  The  only  easily  available  data  are  a  twenty- 
four  hour  precipitation  amount  based  on  the  standard  climatological  day  (06Z- 
06Z).  Precipitation  amounts  were  not  used  for  climate  stations  because  the 
times  of  observation  are  so  variable  and  difficulty  would  be  experienced  in 
associating  the  precipitation  to  the  map  type  times.  No  attempt  was  made  to 
distinguish  convective  and  continuous  precipitation .  Snows  of  early  May  and 
late  September  were  also  not  differentiated. 
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FIGURE  5 

Stations  for  Precipitation  Study  and  Climatic  Probabilities 
of  Measureable  Precipitation 
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The  use  of  twenty-four  hour  precipitation  amounts  required 
some  decision  as  to  the  map  type  time  to  be  used  in  computing  PoPs.  The  12Z 
(6  AM  MDT)  map  type  was  chosen  to  compute  the  probabilities  as  this  was  only 
six  hours  after  the  beginning  of  the  precipitation  day  and  was  usually  prior  to  the 
development  of  convective  precipitation. 

Various  correlations  for  obtaining  PoPs  from  the  data  and 

map  types  are: 

(  i)  12Z  500  mb.  type 
(  ii)  12Z  surface  type 

(iii)  1 2Z  500  mb.  and  surface  types  (referred  to  as  500  mb.  -  surface) 

(  iv)  12Z  500  mb.  and  following  00Z  500  mb.  types  (referred  to  as  500  mb.  - 
500  mb.) 

(  v)  12Z  surface  and  following  00Z  surface  types 

(  vi)  12Z  500  mb.  and  surface  and  following  00Z  500  mb.  and  surface  types. 

PoPs  were  calculated  using  the  first  four  methods  (i  -  iv).  The  first  two  were 
not  expected  to  give  the  best  results  but  hopefully  would  provide  probabilities 
which  were  significantly  different  from  climatological  probabilities.  The 
combination  of  the  500  mb.  -  surface  types  would  be  expected  to  produce  the 
best  results  since  the  thermal  pattern  is  being  taken  into  account.  Hage  (1957) 
has  found  that  abrupt  changes  in  the  500  mb.  flow  occur  during  intense  lee 
cyclogenesis  in  Alberta.  This  was  the  principal  reason  for  computing  PoPs  using 
the  12Z  500  mb.  and  00Z  500  mb.  map  types  and  should  produce  better  results 
than  considering  the  12Z  500  mb.  type  only. 
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The  occurrence  of  precipitation  was  defined  using  a 
threshold  value  of  0.01  inches  or  more,  and  was  based  on  work  by  other  authors. 

It  may  be  more  useful  to  increase  this  value  or  compute  PoPs  for  different  ranges 
of  precipitation  values.  The  present  objective  of  this  work,  however,  was  to 
show  that  significant  improvement  over  climatological  probabilities  is  obtained  by 
taking  the  synoptic  flow  pattern  into  consideration. 

A  map  for  any  given  day  can  correlate  above  the  threshold 
correlation  coefficient  with  one  or  more  of  the  catalogue  types.  To  take  this 
into  account,  the  highest  and  second  highest  correlations  were  computed  and 
retained.  This  would  indicate  that  the  pattern  for  thisparticular  day  is  a  combination 
of  two  map  types  and  could  be  classed  as  a  "hybrid"  type.  In  obtaining  the 
PoPs,  two  approaches  can  be  taken.  First,  the  pattern  can  be  considered  unique 
and  PoPs  can  be  computed  for  this  "hybrid"  type,  or,  the  precipitation 
occurrence  or  non-occurrence  can  be  associated  with  each  type  separately.  The 
latter  method  would  have  the  effect  of  increasing  the  data  base  and  stabilizing 
the  probabilities.  This  method  would  prove  to  be  beneficial  provided  there  was 
no  drastic  increase  or  decrease  in  the  probabilities.  A  test  was  carried  out  when 
the  PoPs  were  computed  for  the  500  mb.  map  types.  Initially,  only  the  primary 
map  type  was  used.  This  provided  a  data  base  of  1,377  observations.  The 
probabilities  were  then  recomputed  using  both  primary  and  secondary  type,  if 
any0  This  procedure  increased  the  data  base  to  2,24  ' .  An  inspection  of  both 
results  indicated  an  average  mean  change  in  the  probabilities  of  approximately 
five  percent,  with  changes  rarely  in  excess  of  ten  percent.  The  largest  changes 
occurred  for  those  map  types  which  had  relativelyfew  occurrences.  Because  the 
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results  of  the  above  showed  a  large  increase  in  the  data  base  and  little  change  in 
the  probabilities,  it  was  decided  that  the  use  of  both  the  primary,  and  secondary  map 
type,  if  any,  would  be  most  beneficial. 

Table  6  indicates  the  increase  in  the  data  base  when  both 
primary  and  secondary  map  types  were  used  for  each  of  the  four  methods. 


TABLE  6 

Increase  in  Data  Base  for  Probabilities  of  Precipitation 


Classification  Method 

Primary 

Types 

Primary  & 
Secondary 

Percent 

Increase 

12Z  500  mb. 

CO 

p-— • 

2246 

63 

12Z  Surface 

1377 

1791 

30 

12Z  500  mb.  &  12Z  Surface 

1377 

2948 

109 

12Z  500  mb.  &  00Z  500  mb. 

1377 

3862 

180 

The  climatological  probabilities  of  measureable 
precipitation  for  each  station  were  calculated  for  the  entire  summer  season  (May 
to  September)  using  the  latest  thirty  year  precipitation  normals  (1941-1970) 
published  by  the  AES.  The  probabilities  for  each  station  for  any  one  day  are 
illustrated  in  Figure  5  (page  52). 
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5.3.2  Results 

The  probabilities  of  precipitation  for  the  500  mb.  map  types 
are  shown  in  Appendix  5.  The  results  indicate  that  significant  departures  from 
climatological  probabilities  are  obtained.  An  investigation  of  the  PoPs  for  the 
most  frequently  occurring  types  indicate  interesting  results.  The  most 
predominant  pattern  is  Map  Type  1 .  This  type  is  characterized  by  a  trough  in 
central  B.C.  with  a  west  to  southwest  flow  over  Alberta.  The  PoPs  show  an 
increased  probability  from  climatology  in  the  far  northern  stations  (YE,  OJ,  MM), 
a  decrease  in  southern  Alberta  and  little  or  no  change  elsewhere.  The  decrease 
in  the  south  is  attributable  to  the  subsidence  flow  over  the  Rocky  Mountains, 
while  the  increase  in  the  north  is  due  to  the  cyclonic  curvature  of  the  flow.  Map 
Type  2  results  show  a  decrease  over  the  southern  half  of  Alberta  as  a  result  of  a 
flat  ridge  predominating.  Map  Type  3  produces  an  increase  over  climatology 
because  of  increasing  cyclonic  flow,  especially  in  the  northern  areas ,while 
decreases  with  Map  Type  8  are  produced  by  ridging. 

The  uncorrelated  group,  which  contains  various  map 
patterns,  also  provides  encouraging  results.  This  type  yields  lower  PoPs  than 
climatology  over  northeastern  Alberta  and  higher  values  in  the  south. 

An  examination  of  all  results  indicate  that  areas  of  high 
and  low  PoPs  are  obtained  by  considering  different  flow  patterns,  with  ranges  of 
probabilities  from  zero  to  one  hundred  percent.  The  values  of  one  hundred 
percent  however,  cannot  be  considered  significant  since  this  type  occurred  only 
five  times.  The  results  appear  reasonable  and  can  be  explained  qualitatively  in 
terms  of  our  knowledge  of  map  patterns  and  precipitation.  It  is  also  obvious  that 
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for  various  map  types,  the  precipitation  probabilities  are  near  the  climatological 
probabilities.  This  indicates  that  more  parameters  need  to  be  investigated. 

The  surface  map  PoPs  are  contained  in  Appendix  6  and  the 
results  are  similar  to  those  for  the  500  mb.  types.  The  ranges  of  values  are  again 
from  zero  to  one  hundred  percent.  The  values  of  one  hundred  percent  should  not 
be  considered  significant  due  to  the  small  number  of  cases,  and  also  the 
probabilities  for  Map  Types  16,  18,  21,  22  and  29  which  had  five  or  fewer 
occurrences  in  nine  years. 

The  PoPs  for  surface  types  2  and  17  are  especially 
interesting.  It  has  previously  been  mentioned  (sections  3.4  and  5.2)  that  these 
two  patterns  are  highly  similar  with  Type2  being  more  developed.  This 
supposition  about  development  is  confirmed  by  the  PoPs.  The  probabilities  for 
Map  Type  2  indicate  lower  than  climatology  for  southern  Alberta  and  along  the 
Rocky  Mountains  as  far  as  XJ,  and  higher  probabilities  in  far  northern  Alberta. 
The  results  for  Map  Type  17  show  low  probabilities  for  all  stations,  not  only 
significantly  lower  than  climatology,  but  generally  one-half  the  values  of 

Type  2. 

The  results  for  the  uncorrelated  group  are  much  less 
encouraging  than  the  500  mb.  The  surface  uncorrelated  group  is  the  most 
predominant  summer  type,  but  little  or  no  improvement  over  climatology  is 
indicated . 

The  probabilities  for  the  combinations  of  12Z  500  mb.  and 
12Z  surface,  or  12Z  500  mb.  and  following  00Z  500  mb.  ere  contained  in 
Appendices  7  and  8,  respectively.  The  final  results  retain  only  those 
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combinations  which  were  considered  important.  Initially,  the  500  mb.  -  surface 
and  the  500  mb.  -  500  mb.  PoPs  resulted  in  490  and  415  distinct  classifications 
respectively.  Many  had  few  occurrences  so  it  was  decided  to  disregard  any 
classification  which  had  fewer  than  ten  occurrences.  These  cases  v/ere  not 
totally  disregarded,  but  were  used  to  compute  PoPs  when  no  significant  surface  or 
00Z  500  mb.  type  occurred.  An  example  of  the  use  of  these  PoPs  can  best  be 
explained  by  referring  to  Map  Type  1  (Appendix  7).  If  the  500  mb.  is  Map  Type 
1,  the  important  surface  map  classifications  are  Types  2,  3,  4,  5,  6,  7,  9,  12, 
13,  14,  17,  28  and  U.  The  PoPs  for  ail  other  surface  types  had  fewer  than  ten 
occurrences, so  they  were  ail  combined  and  are  listed  in  the  last  classification. 

If  the  surface  type  is  not  one  of  the  important  classification,  the  last  is  the  one  to 
use.  This  procedure  decreased  the  number  of  combinations  to  107  and  108  for  the 
500  mb.  -  surface  and  500  mb.  -  500  mb.,  respectively. 

Classification  of  the  500  mb.  and  surface  pattern 
simultaneously  generally  yields  improved  results  over  considering  the  types 
separately.  This  result  was  expected  since  the  three-dimensional,  i.e.,  thermal 
aspect  of  the  atmosphere  is  being  taken  into  account.  There  are,  however, 
certain  combinations  which  do  not  change  the  probabilities  from  those  of  the  500 
mb.  or  surface.  No  attempt  was  made  to  test  the  significance  of  a  particular 
combination  other  than  by  the  number  of  occurrences.  It  is  again  worthwhile  to 
note  that  differences  in  PoPs  for  surface  May  Types  2  and  17  with  an  upper  air 
pattern  are  still  apparent  as  was  the  case  when  they  were  considered  separately. 

The  combination  of  12Z  500  mb.  and  the  following  00Z  500 
mb.  PoPs  (Appendix  8)  indicate  some  improvement  over  using  the  I2Z  500  mb. 
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type  alone.  When  there  is  no  change  in  type  from  12Z  to  00Z,  the  probabilities 
are  virtually  identical  with  those  of  Appendix  5.  When  there  is  a  change  in 
type,  the  PoPs  increase  or  decrease  in  magnitude  depending  on  the  degree  and 
type  of  change  (i.e.,  trough  becoming  ridge  or  vice  versa)  in  the  pattern. 


5.3.3  PoPs  Verification  Using  Brier  Scores 

In  order  to  evaluate  the  ability  of  the  above  four  methods  to 
forecast  the  probability  of  rain  occurrence,  the  Brier  Score  and  related 
statistics  were  computed  using  the  development  sample.  Brier  Scores  provide  a 
convenient  numerical  verification  of  probability  forecasts.  Included  is  a 
measurement  of  both  "reliability"  and  "resolution"  factors.  The  "reliability"  of 
the  forecasts  is  the  agreement  between  the  forecast  probability  and  observed 
frequency  of  occurrence  of  the  event.  The  "resolution"  measures  the  skill  in 
assigning  to  the  forecast  values  of  probability  which  deviate  from  the 
climatological  frequency. 

In  Appendix  11,  it  is  shown  that  the  average  Brier  Score 
for  M,  forecasts  in  the  probability  category  (i.e.,  0,  10,  20,  .  .  .  ,  100%) 
is: 


Bk  =  <Fk  -  V  +\v-V' 


(2) 


where  is  the  forecast  probability  and  0^  is  the  actual  relative  frequency  of 
occurrence.  The  first  term  on  the  right  hand  side  of  (2)  is  the  reliability  term, 
while  the  second  is  the  resolution.  The  reliability  term  becomes  a  minimum  when 
equals  0^.  On  the  other  hand,  the  resolution  term  becomes  a  minimum  when 
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0^  is  either  zero  or  one  hundred  percent  (i.e,,  equal  to  one)  and  a  maximum  at 
0^  equals  fifty  percent.  It  can  be  seen  that  if  the  reliability  is  fairly  good  (i.e., 
that  term  is  small),  the  poorest  separation  of  forecasts  (resolution  term)  results  in 
the  highest  Brier  Score. 


The  average  Brier  Score  for  all  probability  categories  is 
obtained  by  averaging  the  mean  scores  for  each  category,  weighted  by  the  number 
of  forecasts  in  that  category.  Thus, 


B  =  -jir  L  Mk  Bk  -  (3) 

where  M  equals  the  total  number  of  forecasts.  Possible  values  of  "B"  range 
from  zero  to  one,  with  zero  the  best  score. 

Two  additional  statistics  which  are  useful  in  comparing 
forecasts  are  the  root-mean-square  error  in  reliability  (RMSE)  and  the  skill  score 
based  on  a  comparison  with  a  simple  control  forecast.  The  first  of  these  is 
defined  by: 

RMSE  =121^  MkW  -  0k}2  .  (4) 

and  the  skill  score  by: 

s  =  100  { 1  -  7W } 
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where  B  is  the  average  Brier  Score  obtained  by  employing  the  forecast 
probabilities  and  cB,  the  Brier  Score  obtained  by  using  the  climatological 
frequency  or  persistence  as  the  forecast  probability. 

All  the  above  parameters  were  computed  for  each  of  the 
PoPs  obtained.  Eleven  subsets  of  forecast  probabilities  were  defined  and  are 
indicated  in  Table  7.  In  applying  the  derived  PoPs,  only  the  primary  map  types 
v/ ere  used  in  arriving  at  a  probability  forecast,  whereas  both  primary  and 
secondary  types  had  been  used  in  the  derivation.  Skill  scores  were  computed  using 
both  climatological  frequency  (based  on  period  May  to  September)  and  persistence 
as  the  control  forecast.  It  was  found  that  the  map  type  PoPs  were  better  than 


Table  7 

Probability  Subsets  and  Ranges  for  Brier  Score  Computations 

Probability  Range 

Probability  Subset  K  for  Subset 


0 

0  - 

4 

10 

5  - 

14 

20 

15  - 

24 

30 

25  - 

34 

40 

35  - 

44 

50 

45  - 

54 

60 

55  - 

64 

70 

65  - 

74 

80 

75  - 

84 

90 

85  - 

94 

100 

95  - 

100 

' 
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persistence  by  a  substantia!  amount,  however,  climatology  was  also  better  than 
persistence.  As  a  result,  persistence  was  disregarded  as  a  control  forecast  for 
summer  precipitation.  The  calculated  RMSEs  are  exceedingly  small  (less  than  0.5 
percent)  for  all  four  methods  and  for  all  stations  as  expected.  The  reliability 
term,  therefore,  contributes  little  to  the  Brier  Score  and  it  is  the  resolution  term 
which  is  the  dominant  factor.  Table  8  illustrates  the  computed  skill  scores  for 
each  station  and  each  map  type  PoPs  classification  using  the  climatological 
frequency  of  measureable  precipitation  as  the  control  forecast.  The  underlined 
values  for  each  station  refer  to  the  highest  computed  skill  score  while  the  dashed 
values  are  scores  within  one  percent  of  the  highest  value. 

The  best  improvement  over  climatology  is  obtained  for  most 
stations  by  using  the  combination  of  12Z  500  mb.  and  12Z  surface  types.  Only 
seven  stations  differed  from  the  above  and  five  of  these  were  using  the 
classification  method  of  12Z  500  mb.  and  following  00Z  500  mb.  map  type.  The 
tv/o  remaining  stations  (CT  and  RM)  obtained  the  best  results  from  using  only  the 
12Z  surface  type.  It  should  be  noted  however  that  for  these  seven  stations,  the 
500  mb.  -  surface  combination  skill  scores  were  within  one  percent  of  the  highest 
results.  The  500  mb.  -  500  mb.  classification  is  preferable  to  using  only  the 
12Z  500  mb.  type  in  all  cases.  This  result  v/as  anticipated  since  a  twenty-four 
hour  precipitation  period  v/as  used  but  this  may  not  be  the  case  if  a  shorter  time 
period  is  used . 
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TABLE  8 

Skill  Scores  (Percent)  Using  Primary  Map  Types  for  Development  Sample 


(May  - 

September,  1963-1971) 

Station 

500  mb. 

500  mb.  - 
500  mb. 

Surface 

500  mb. 
-  Surface 

XJ 

8.8 

7.5 

10.5 

12.6 

YE 

5.0 

5.2 

8.5 

8.0 

CT 

9.9 

14.4 

11.2 

13.4 

EG 

9.9 

11.7 

11.2 

13.3 

XD 

9.2 

11.2 

10.9 

12.8 

RM 

11.9 

17.7 

14.2 

16.8 

VG 

10.5 

9.3 

11.7 

16.5 

QF 

13.2 

16,0 

16.0 

18.5 

YC 

16.3 

19.5 

21.7 

20.8 

QL 

16.9 

16.3 

19,8 

21,8 

XH 

13.5 

11.1 

16.5 

16.8 

PC 

19.2 

17.6 

22.2 

24.5 

BA 

17.9 

12.4 

18.9 

21.8 

JA 

10.8 

8.1 

12.8 

14.9 

ET 

11.0 

13.4 

12.3 

16.3 

MM 

7.2 

6.8 

11.3 

10.6 

ZH 

10.5 

7.2 

13.8 

14.2 

ZU 

8.9 

12.2 

9.2 

14.5 

PY 

6.3 

8.3 

8.7 

9.8 

OJ 

6.5 

4.7 

8.4 

10.0 

QU 

9.3 

8.1 

10.6 

14.5 

PE 

8.4 

7.3 

11.9 

n.7 

OD 

10.3 

7.9 

13.7 

13.4 
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It  is  also  interesting  to  note  that  the  greatest  improvements 
over  climatology  is  obtained  for  stations  in  central  and  southern  Alberta. 

This  result  indicates  that  the  synoptic  flow  is  an  important  factor  for 
precipitation  in  central  and  southern  Alberta.  This  may  be  related  to  the  height 
and  proximity  of  the  Rocky  Mountains. 

A  comparison  of  the  observed  frequency  of  precipitation  to 
forecast  probability  for  each  probability  category  (not  shown)  indicates  that  for 
low  probabilities,  the  forecast  probabilities  are  generally  too  low,  while  in 
higher  probabilities,  the  forecast  values  are  too  high.  This  problem  has  been 
reported  in  previous  studies,  e.g.,  Saunders  (1963),  and  represents  "over-confidence" 

by  the  method,  that  is,  the  method  over-estimates  its  skill  in 

separating  the  forecasts  into  the  extreme  ends  of  the  probability  range.  This  may 

be  the  result  of  too  few  cases  in  the  derived  PcPs. 

The  skill  scores  for  any  method  show  only  moderate 
improvements  over  climatology  at  any  station.  This  is  because  too  many  forecasts 
are  made  at  the  climatological  frequency  of  precipitation.  Ideally,  the  forecast 
probability  distribution  should  have  maxima  above  and  below  the 
climatological  probabilities  (ideally  at  zero  and  one  hundred  percent).  This 
occurs  only  for  those  stations  having  a  skill  score  of  greater  than  twenty  percent 
and  even  in  these  cases,  the  maximum  number  of  forecasts  is  much  greater  at  the 
low  ranges  than  at  the  high  ranges.  An  example  for  Lethbridge  (QL)  is 
illustrated  in  Figure  8,  Appendix  11  . 
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These  results  indicate  that  for  the  development  sample 
(May  to  September  1963-1971),  improvement  over  climatological  precipitation 
probabilities  is  obtained  by  taking  the  synoptic  flow  pattern  into  consideration 
PoPs  verifications  using  map  types  on  an  independent  data  sample  (May  - 
September  1972-1973)  were  calculated.  They  are  discussed  in  the  following 
chapter. 


• '  l 


CHAPTER  VI 


RESULTS  FOR  INDEPENDENT  DATA  SAMPLES 

6.1  Surface  and  500  mb.  Map  Types  1972  and  1973 

The  independent  data  sample,  to  test  the  results  of  the  map 
catalogues  and  PoPs,  consisted  of  612  days  for  the  period  May  to  September  1972 
and  1973.  The  data  for  this  period  consisted  of  500  mb.  heights  and  mean  sea 
level  pressure  values  covering  the  area  of  the  map  catalogues.  Twenty-four  hour 
precipitation  amounts  for  the  periods  06Z  -  06Z  and  18Z  -  18Z  were  obtained 
from  Monthly  Summaries  and  from  the  Edmonton  Weather  Office.  The 
verification  results  for  the  PoPs  are  discussed  in  section  6.2. 

The  map  times  of  00Z  and  12Z  for  each  day  were  correlated 
against  the  key  map  pattern  for  the  surface  and  500  mb.  catalogues.  Each  map 
was  associated  with  that  map  type  for  which  the  correlation  coefficient  was  the 
highest  provided  it  exceeded  or  equalled  the  threshold  values  of  0,80  for  the 
surface  or  0.85  for  the  500  mb.  types.  Any  map  which  correlated  below  the 
threshold  value  v/ith  all  key  maps  was  assigned  to  an  uncorrelated  group.  The 
frequencies  of  occurrence  of  each  map  type  during  the  independent  data  period 
are  presented  in  Table  9. 
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TABLE  9 

Number  and  Frequency  of  Occurrence  (%)  of  Primary  500  mb.  and 
Surface  Map  Types  for  May  -  September  1972-1973 
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Map  Type 

500  mb. 

Surface 

1 

113(18) 

13(  2) 

2 

76(12) 

47 (  8) 

3 

47(  8) 

10(  2) 

4 

23(  4) 

73(12) 

5 

14(  2) 

12(  2) 

6 

8(  1) 

13(  2) 

7 

4(  1) 

14(  2) 

8 

55  (  9) 

9(  1) 

9 

10(  2) 

35  (  6) 

10 

K  0) 

8(  1) 

11 

36(  6) 

1(  0) 

12 

5(  1) 

12(  2) 

13 

3(  0) 

21  (  3) 

14 

17(  3) 

6(  1) 

15 

11(  2) 

2(  0) 

16 

7(  1) 

1(  0) 

17 

1(  0) 

7(  1) 

18 

6(  1) 

- 

19 

5(  1) 

3(  0) 

20 

6(  1) 

3(  0) 

21 

2(  0) 

- 

22 

2(  0) 

- 

23 

2(  0) 

4(  1) 

24 

9(  1) 

15(  2) 

25 

9(  1) 

15(  2) 

26 

8(  1) 

12(  2) 

27 

2(  0) 

12(  2) 

28 

5(  1) 

1(  0) 

29 

2(  0) 

2(  0) 

30 

3(  0) 

5(  1) 

31 

2(  0) 

2(  0) 

32 

5(  1) 

2(  0) 

33 

5(  1) 

- 

34  or  U 

1(  0) 

265(43) 

35 

7(  1) 

U 

109(18) 

* 

-  - 

- 

' 

* 

i.  £1  - 
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The  results  for  500  mb.  Indicate  very  little  change  in  map  type 
frequency  between  the  dependent  and  independent  data  samples.  The  uncorrelated 
group  frequency  shows  the  nighest  increase  of  six  percent,  making  it  the  second  most 

recurring  map  type.  Map  Type  1  remains  the  predominant  pattern  with  the  five  most 

frequent  types  remaining  the  same  in  both  samples  although  their  relative  order  has  changed . 
The  surface  map  types  indicate  very  little  difference  in 

frequencies  of  occurrence  between  dependent  and  independent  data  samples,  the 
mean  change  being  one  percent.  As  with  the  500  mb.,  the  uncorrelated  group 
frequency  shows  an  increase  of  nine  percent,  with  this  group  remaining  the 
predominant  type.  The  five  most  predominant  types  and  their  relative  order  have 
shown  no  change,  while  Types  18,  21,  22  and  33  have  not  occurred. 

6.2  PoPs  Verification  for  1972  and  1973 

The  precipitation  probabilities  v/ere  verified  for  each  of  the 
derived  methods  based  on  the  Brier  Score  and  using  the  same  procedure  as  in 
section  5.3.3.  The  occurrence  of  precipitation  was  defined  as  an  amount  of  0.01 
inches  or  more. 

The  computed  root-mean-square  errors  in  the  reliability  (RMSE)  and 
the  skill  scores  for  the  precipitation  period  06Z-06Z  are  contained  in  Tables  10 
and  1 1 ,  respectively . 
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TABLE  10 

Root-Mean-Square  Errors  in  Reliability  (%)  for  Independent  Data 

(Precipitation  Period  06Z-06Z) 


Station 

500  mb . 

500  mb .  - 
500  mb. 

Surface 

500  mb. 
-  Surface 

XJ 

0.3 

0.3 

0.3 

0.5 

YE 

1.5 

1.8 

1.3 

1.4 

CT 

0.7 

1.0 

0.8 

1.1 

EG 

1.6 

1.2 

1.9 

0.8 

XD 

1.8 

1.9 

1.4 

1.6 

RM 

1.4 

0.8 

•  0.6 

0.5 

VG 

1.1 

1.2 

0.8 

0.9 

QF 

1.6 

0.7 

0.5 

1.7 

YC 

1.5 

1.0 

0.8 

1.2 

QL 

0.4 

0.4 

0.6 

0,3 

XH 

1.0 

0.7 

0.7 

1.7 

PC 

3.0 

1.7 

0.8 

1.6 

BA 

1.0 

1.1 

1.1 

1.4 

JA 

0.2 

0.9 

1.1 

0.8 

ET 

1.2 

1.4 

1.0 

1.9 

MM 

1.8 

1.9 

1.1 

1.6 

ZH 

2,3 

1.4 

1.7 

1.4 

ZU 

2.5 

2.5 

3.1 

2.5 

PY 

0.9 

2.5 

0.3 

1.3 

OJ 

1.0 

1.1 

1.5 

1.2 

QU 

0.8 

1.2 

1.1 

1.3 

PE 

1.2 

1.7 

1.6 

2.1 

OD 

1.7 

2.3 

0.8 

1.6 

. 
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TABLE  11 

Skill  Scores  for  Independent  Data  (Precipitation  Period  06Z-06Z) 


Station 

500  mb . 

500  mb .-  500  mb . 

Surface 

500  mb.  - 
Surface 

XJ 

3.0 

3.5 

5.9 

4.4 

YE 

-  0.5 

-  2.4 

2.8 

0.7 

CT 

5.3 

5.9 

5.7 

4.1 

EG 

4.0 

4.2 

3.0 

6.0 

XD 

1.3 

0.3 

2.3 

2.2 

RM 

10.2 

13.6 

20.6 

12.5 

VG 

2.1 

1.7 

6.9 

2.0 

QF 

5.2 

9.4 

12.1 

5.2 

YC 

13.0 

17.8 

12.8 

17.3 

QL 

19.7 

19.1 

in 

00 

22.1 

XH 

7.6 

7.4 

7.0 

10.6 

PC 

5.0 

9.5 

6.6 

8.9 

BA 

14.4 

16.0 

6.5 

13.9 

JA 

6.4 

4.9 

1.5 

CO 

ET 

2.9 

3.1 

6.7 

3.1 

MM 

-  0.1 

1.0 

10.4 

3.6 

ZH 

-  0.7 

2.0 

5.6 

4.1 

ZU 

-  2.7 

-  1.9 

0.8 

-  2.2 

PY 

-  1.5 

-  2.6 

7.3 

1.5 

OJ 

2.7 

3.7 

3.4 

3.4 

QU 

2.3 

4.1 

1.8 

0.2 

PE 

2.3 

o  1 
00  1 

3.2 

-  0.6 

OD 

0.2 

-  3.5 

7.9 

2.3 

•  -  -  * 
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The  RMSEs  for  all  stations  have  shown  an  increase  over  these  of  the 
development  sample  but  the  maximum  value  is  less  than  3.5  percent.  This 
remains  wel I  below  the  average  value  of  11.3  percent  reported  by  Nikleva  (1971) 
for  subjective  probability  forecasts  for  six  fire  weather  forecast  locations  for  the 
period  May  to  September  1965  to  1969.  The  present  study's  results  appear  to  be 
highly  favourable  in  this  regard. 

The  low  RMSE  values  indicate  that  the  derived  PoPs  are  reliable. 
The  contribution  of  the  reliability  term  to  the  total  Brier  Score  is  less  than  ten 
percent,  the  maximum  being  thirteen  percent  for  Whitecourt  (ZU)  using  the 
surface  map  type  alone. 

The  computed  skill  scores  (Table  11)  exhibit  some  interesting 
results.  The  most  noticeable  aspect  is  the  occurrence  of  negative  skill  scores, 
implying  that  no  improvement  over  climatology  has  been  obtained  by  using  map 
types.  It  has  previously  been  mentioned  that  the  improvement  over  climatology 
for  stations  north  of  a  line  ET-XD-VG  was  generally  much  lower  than  those  for 
southern  stations.  The  negative  skill  scores  reaffirm  this  pattern.  Also,  the 
classification  method  producing  the  highest  skill  scores  are  greatly  different. 

The  underlined  values  in  Table  11  indicate  the  highest  skill  score  for  that 
station  and  the  dashed  line  indicates  differences  of  less  than  one  percent  from 
the  highest  score.  For  the  dependent  data,  the  best  combination  of  types 
appears  to  be  the  500  mb.  -  surface,  since  these  prod  iced  the  highest  scores 
most  frequently.  For  the  independent  data,  this  combination  produced  the 
highest  scores  for  only  three  stations  (EG,  QL,  and  XH)  compared  with  sixteen 
stations  previously.  Similarly,  the  surface  map  type  produces  the  best  results  for 
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Rocky  Mountain  House  (RM)  for  both  data  samples  while  for  Coronation  (CT), 
this  is  true  for  the  dependent  data  only.  The  500  mb.  -  500  mb.  and  surface 
map  type  skill  scores  for  CT  differ  by  only  0.2  percent  however.  The  motion  of 
the  upper  air  pattern  (500mb.  -  500  mb.  combination)  results  in  the  highest 
verification  for  Calgary  (YC)  in  both  data  samples.  All  stations  except  RM  and 
QL  experienced  decreases  in  skill  scores  compared  to  the  dependent  data,  with 
QL  continuing  to  show  the  highest  results  of  all. 

It  is  of  interest  to  point  out  the  results  for  XJ  and  QL.  The 
RMSEs  in  reliability  (Table  10)  are  nearly  identical  for  both  stations  for  all 
methods,  averaging  0.4  percent.  The  skill  scores  (Table  11)  for  the  same  two 
stations  show  substantial  differences,  average  4.2  percent  for  XJ  compared  v/ith 
17.4  percent  for  QL.  These  differences  in  the  skill  scores  are  the  result  of  the 
resolution  of  the  forecasts.  A  greater  proportion  of  the  QL  forecasts  are  above 
or  below  the  climatological  frequency  of  measureable  precipitation  compared  to 

XJ. 

Since  only  a  twenty-four  hour  precipitation  amount  was  available, 
one-half  of  the  612  maps  in  the  independent  data  sample  could  not  be  used  in 
the  previous  PoPs  verification  (namely,  00Z  map  types).  A  twenty-four  hour 
precipitation  amount  for  the  period  18Z-18Z  was  available  for  1972-1973  but 
not  for  the  development  sample.  The  derived  PoPs  were  used  for  this 
precipitation  period  with  the  00Z  map  type  replacing  the  12Z  type.  If  the 
verification  results  were  acceptable,  a  benefit  would  be  obtained  since  all  map 
types  could  be  used. 
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The  RMSEs  and  the  skill  scores  for  the  precipitation  period 
18Z-18Z  are  shown  in  Tables  12  and  13,  respectively.  The  RMSEs  in  reliability 
are  similar  to  those  in  Table  11,  implying  that  the  reliability  term  contributes 
little  to  the  total  Brier  Score. 

The  skill  scores  indicate  differences  between  the  development 
data  and  the  results  discussed  above.  The  skill  scores  are  lower  than  the 
dependent  data  values  for  all  stations  while  there  has  been  increases  and 
decreases  in  scores  from  the  1972-1973  data  using  the  12Z  types.  Medicine  Hat 
(XH)  is  the  only  station  which  has  consistently  experienced  the  highest  value  with 
one  map  type  combination,  namely  the  500  mb.  -  surface.  The  high  scores  are 
now  mostly  produced  by  the  classifications  using  either  the  500  mb.  only  or  the 
500  mb.  -  500  mb.  types  (thirteen  of  twenty-three);  the  greatest  increases 
occurring  for  the  500  mb.  type  only.  This  effect  may  indicate  that  convective 
precipitation  in  summer  is  controlled  by  preferred  500  mb.  flow  patterns. 

An  investigation  of  the  number  of  forecasts  in  each  probability 
category  for  all  stations  indicate  that  there  are  many  forecasts  at  or  near  the 
climatic  frequency  of  measureable  precipitation.  This  explains  the  lack  of  great 
improvement  over  climatology  and  indicates  that  more  variables  are  required. 
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TABLE  12 

Root-Mean -Square  Errors  in  Reliability  (%)  for  Independent  Data 

(Precipitation  Period  18Z-18Z) 


Station 

5C0  mb. 

500  mb. - 
500  mb. 

Surface 

500  mb. 
-  Surface 

XJ 

0.2 

0.5 

1.0 

1.0 

YE 

0.6 

1.2 

1.1 

1.0 

CT 

0.4 

1.0 

1.3 

1.2 

EG 

0.3 

0.9 

0.6 

0.5 

XD 

0.6 

0.5 

0.8 

0.7 

RM 

1.1 

1.0 

•  1.1 

0.8 

VG 

1.6 

1.0 

1.0 

1.5 

QF 

1.2 

0.7 

0.6 

1.3 

YC 

1.5 

1.1 

0.5 

1.6 

QL 

0.6 

0.8 

1.8 

1.1 

XH 

1.2 

0.4 

0.5 

1.3 

PC 

1.3 

0.3 

0.8 

0.8 

BA 

0.7 

0.7 

0.9 

0.4 

JA 

0.4 

1.2 

0.6 

0.9 

ET 

1.0 

1.2 

0.6 

0.7 

MM 

1.4 

1.6 

0.7 

1.3 

ZH 

1.3 

1.0 

0.8 

2.4 

ZU 

1.9 

1.6 

1.1 

2.1 

PY 

0.5 

0.7 

0.4 

1.1 

OJ 

0.6 

1.6 

0.5 

1.0 

QU 

0.4 

0.7 

1.5 

0.8 

PE 

1.0 

0.9 

0.9 

1.2 

OD 

1.6 

1.4 

0.7 

2.1 

* 
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TABLE  13 

Skill  Scores  (%)  for  Independent  Data 
(Precipitation  Period  18Z-18Z) 


Station 

500  mb. 

500  mb. - 
500  mb. 

Surface 

500  mb. 
-  Surface 

XJ 

4.0 

2.3 

1.6 

0.2 

YE 

6.5 

3.6 

5.3 

6.6 

CT 

7.2 

6.1 

-  0.1 

2.8 

EG 

11.3 

11.7 

6.2 

8.4 

XD 

6.6 

7.5 

6.7 

5.3 

RM 

13.4 

12.5 

11.3 

12.0 

VG 

0.6 

3.2 

3.0 

0.6 

QF 

9.1 

12.6 

8.0 

9.1 

YC 

17.9 

16.2 

9.5 

17.7 

QL 

20.3 

18.4 

-  2.0 

13.6 

XH 

8.7 

7.4 

2.7 

9.3 

PC 

13.0 

13.3 

5.2 

13.4 

BA 

13.8 

18.7 

10.7 

12.7 

JA 

10.3 

4.6 

3.3 

3.5 

ET 

4.6 

5.5 

11.5 

4.3 

MM 

0.4 

-  0.4 

9.5 

• 

CO 

ZH 

6.0 

6.4 

10.1 

6.8 

ZU 

-  0.3 

1.9 

5.4 

-  2.0 

PY 

1.3 

2.6 

4.6 

-  0.3 

OJ 

1.3 

-  2.8 

2.3 

0.1 

QU 

6.6 

4.8 

1.4 

4.1 

PE 

7.7 

5.3 

4.3 

5.1 

OD 

-  2.4 

-  2.9 

2.6 

-  2.0 

- 
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When  the  PoPs  statistics  were  initially  computed  for  the  500  mb.  - 
surface  and  500  mb.  -  500  mb.,  490  and  415  combinafions,respectiveIy,  had 
occurred  (section  5.3.2).  It  was  arbitrarily  decided  that  any  combination  having 
less  than  ten  occurrences  was  probably  not  significant  and  these  cases  were  not 
retained  as  a  separate  entity.  In  order  to  test  the  effect  of  this  decision, 
verification  statistics  were  computed  for  1972-1973  using  all  possible 
combinations  for  the  500  mb.  -  surface  classification.  The  results  indicate  a 
degradation  in  the  skill  scores  for  all  except  a  few  stations.  Negative  skill 
scores  are  obtained  for  a  greater  proportion  of  stations,  even  for  those  which  had 
positive  skill  scores  previously.  There  is,  however,  some  improvement  for  three 
stations  in  both  the  06Z-06Z  and  18Z-18Z  precipitation  periods.  For  the  period 
06Z-06Z,  RM  (20.0%),  XH  (10.7%),  and  ET  (4.0%)  show  an  increase,  while 
for  18Z-18Z,  RM  (14.5%),  XH  (8.7%)  and  MM  (6.2%).  It  is  especially 
noteworthy  that  RM  consistently  shows  an  improvement,  in  fact,  the  resulting 
skill  scores  would  be  near  or  the  highest  for  the  independent  data.  Although  the 
verification  was  not  performed  using  all  500  mb.  -  500  mb.  combinations,  there 
is  no  reason  to  believe  that  the  results  would  be  much  different. 

The  verification  generally  confirms  the  necessity  of  rejections  of 
combinations  which  have  few  occurrences.  It  is  probable,  however,  that  the 
arbitrary  value  of  ten  in  nine  summers  is  not  the  most  desirable. 

The  verifications  indicate  that  improvement  over  climatic 
frequency  of  precipitation  is  obtained,  notably  for  southern  stations,  by 
classifying  the  data  into  different  flow  patterns.  The  improvement  is  highly 
variable  over  the  area  verified  and  may  be  result  of  the  highly  irregular 
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topographic  features  or  not  distinguishing  between  convective  and  large-scale 
precipitation.  No  one  classification  method  produces  the  best  results  for  every 
station  and  this  may  also  be  a  result  of  topography  or  convective  activity.  The 
derived  probabilities  of  precipitation  for  06Z-06Z  can  also  be  used  for  the  period 
18Z-18Z  provided  the  00Z  map  types  are  used  instead  of  the  12Z  types. 

6.3  PoPs  and  Verification  Using  a  Higher  Precipitation  Threshold 

Initially  PoPs  were  derived  using  a  threshold  of  0.01  inches  or 
greater.  This  value  was  chosen  following  the  work  of  other  authors  but  may  not 
be  realistic.  The  synoptic  flow  pattern  may  not  be  as  dominant  an  influence  on 
low  precipitation  amounts  as  on  higher  values.  In  order  to  provide  an 
indication  of  this  possible  effect,  PoPs  were  calculated  for  the  12Z  500  mb.  and 
12Z  surface  map  types  for  the  development  sample  using  a  threshold  value  of  0.10 
inches  or  greater  and  the  precipitation  period  06Z-06Z.  This  higher  threshold 
value  should  tend  to  eliminate  the  precipitation  occurrences  due  to  weak  synoptic 
features  and  convective  activity. 

Both  primary  and  secondary  map  types  were  used  to  obtain  the  PoPs 
to  increase  .the  data  base  and  obtain  more  stable  probabilities.  The  climatic 
frequency  of  precipitation  occurrence  was  computed  for  this  threshold  value  and 
are  shown  in  Table  14.  The  values  are  relatively  constant  for  all  stations  and 
show  a  decreasing  tendency  towards  the  south. 
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TABLE  14 

Climatic  PoPs  (  >  0.10  in.)  and  Verification  Skill  Scores  For 
May  -  September,  1972-1973  (Precipitation  Periods  06Z-06Z  and  18Z-18Z) 


Station 

Climatic  PoPs 

]  2Z  500  mb . 

12Z  Sfc . 

00Z  500 

00Z  Sfc . 

XJ 

17 

-  3.0 

-  4.6 

5.0 

3.2 

YE 

17 

-  3.6 

-  0.6 

-  0.2 

2.5 

CT 

12 

2.0 

8.0 

3.8 

-  1.5 

EG 

17 

1.4 

6.2 

2.5 

5.1 

XD 

17 

1.1 

5.6 

4.6 

4.5 

RM 

22 

7.1 

15.3 

9.7 

4.6 

VG 

17 

1.8 

3.5 

-  0.8 

1.2 

QF 

18 

6.0 

11.3 

7.3 

3.8 

YC 

15 

7.2 

11.2 

6.9 

5.2 

QL 

12 

10.4 

8.3 

8.0 

0.1 

XH 

10 

1.2 

4.8 

4.1 

2.1 

PC 

13 

3.8 

9.7 

2.6 

10.6 

BA 

17 

3.2 

4.0 

-  0.2 

4.1 

JA 

15 

4.3 

-  0.6 

4.2 

2.0 

ET 

22 

2.2 

-  1.9 

4.5 

5.9 

MM 

21 

-  3.9 

0.4 

0.6 

3.1 

ZH 

20 

-  2.5 

0.6 

2.7 

3.7 

ZU 

19 

1.7 

K5 

0.5 

4.2 

PY 

17 

-  1.6 

2.2 

-  0.4 

2.9 

OJ 

15 

1.1 

-  2.8 

1.3 

-  0.9 

QU 

15 

0.6 

-  4.9 

5.7 

0.1 

PE 

14 

-  0.3 

-  3.4 

7.8 

-  1.2 

OD 

18 

-  2.1 

4.8 

-  0.2 

0.7 

. 


•  .  •  .  -  - 

• 

* 

.  -  • 

* 

' 


■  .  - 


*  r 

.  1 

. 

i  *  *  •  *  •  l 


79 

The  resulting  PoPs  for  both  500  mb.  and  surface  are  contained  in 
Appendices  9  and  10,  respectively.  The  probabilities  are  lower  than  those  at 
0.01  inches  but  the  departures  from  the  climatic  values  are  similar.  Areas  of  high 
and  low  probability  are  grouped  together  in  different  areas  depending  upon  the 
flow  pattern.  The  ranges  of  values  are  similar  for  both  types,  namely  from  zero 
to  eighty  percent,  with  the  low  values  predominating.  The  PoPs  again  exhibit  the 
tendency  of  being  too  close  to  the  climatic  PoPs.  No  attempt  was  made  to  repeat 
the  computations  for  combinations  of  map  patterns. 

Verification  results  were  computed  for  1972-1973,  using  both 
precipitation  periods  as  before.  The  RMSEs  (not  shov/n)  are  less  than  five  percent, 
indicating  that  the  forecast  probabilities  are  reliable.  The  skill  scores  are  shown 
in  Table  14  with  the  underlined  values  being  the  highest  value  and  the  dashed 
line,  within  one  percent  of  the  highest. 

The  skill  scores  for  the  12Z  types  are  highly  similar  to  previous 
results  for  the  majority  of  stations.  There  are  notable  exceptions  such  as  the 
improvement  at  ZU  and  the  degradation  at  Fort  St.  John  (XJ)  and  Fort  Nelson 
(YE).  The  best  results  are  evident  once  more  in  southern  stations  while  northern 
ones  produce  most  of  the  negative  scores.  The  greatest  proportion  of  stations 
experience  the  best  scores  with  the  surface  map  type,  notably  RM  (15.3%)  as  was 
the  case  previously.  The  skill  scores  for  the  precipitation  period  18Z-18Z  yield 
similar  results  as  for  the  0.01  inch  threshold  (Table  14).  The  highest  values  are 
nearly  equally  divided  between  the  500  mb.  and  surface  patterns.  Shifts  in  the 
highest  values  are  evident  along  with  increases  and  decreases  in  values.  The 
results  indicate  that  the  derived  PoPs  are  suitable  for  the  18Z-18Z  precipitation 
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period  provided  the  00Z  map  type  is  used. 

Although  the  evidence  is  not  conclusive,  it  appears  that  more 
experimentation  is  required  in  deciding  on  the  threshold  precipitation  value. 
Significant  precipitation,  including  convective  activity,  would  be  more 
dependent  on  the  synoptic  flow  pattern  and  therefore  a  higher  threshold  value 
may  be  preferrable.  The  question  of  what  is  significant  precipitation  would 
depend  on  the  type  of  climatological  aid  that  is  required. 

6.4  Comparing  Summer  500  mb.  Map  Types  to  Those  of  Winter 

Although  the  500  mb.  map  types  were  produced  from  the  May  to 
September  period,  weather  patterns  which  occur  in  other  seasons  should  bear  a 
resemblance  to  the  map  types  already  obtained,  probable  differences  being 
evident  in  frequencies  of  occurrence  and  other  map  type  characteristics.  Upper 
air  map  types  were  therefore  computed  using  the  predetermined  summer  key 
patterns  for  a  small  data  sample  of  281  maps  (OOZ  and  12Z)  for  the  period 
November  1,  1973  to  March  21,  1974  inclusive. 

The  number  of  occurrences  and  corresponding  frequency  for  each 
map  type  is  shown  in  Table  15.  A  total  of  thirteen  types  did  not  occur  during 
this  data  sample  and  is  probably  due  to  the  small  data  period  used.  Some  may  be 
missing  because  the  pattern  is  not  possible  at  this  time  of  year,  the  result  of  the 
shift  in  the  thermal  gradient  southwards.  Definite  conclusions  require  a  larger 
data  sample  to  be  processed. 
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TABLE  15 

500  mb.  Map  Type  Occurrences  and  Frequency  (%)  (Winter,  1973-1974) 


81 


500  mb.  Map  Type 

Number  of  Occurrences  (Frequency  %) 

1 

68(24) 

2 

38(14) 

3 

14(  5) 

4 

9(  3) 

5 

7(  2) 

6 

8(  3) 

7 

- 

8 

14(  5) 

9 

- 

10 

12(  4) 

11 

29(10) 

12 

- 

13 

- 

14 

24(  9) 

15 

H  0) 

16 

5(  2) 

17 

U  0) 

18 

- 

19 

5(  2) 

20 

— 

21 

- 

22 

- 

23 

H  0) 

24 

- 

25 

5(  2) 

26 

2(  1) 

27 

- 

28 

- 

29 

K  0) 

30 

K  o) 

31 

- 

32 

6(  2) 

33 

K  0) 

34 

- 

35 

2(  1) 

U 

27(10) 
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Changes  in  the  relative  frequencies  of  occurrence  between  summer 
and  winter  are  suprising!/  very  minor.  Map  Type  1  frequency  has  increased  by 
five  percent  and  this  pattern  continues  to  be  the  predominant  one.  The 
uncorrelated  group  frequency  has,  on  the  other  hand,  decreased  slightly.  The 
increased  frequency  of  Map  Types  11  and  14  has  resulted  in  their  being  included 
in  the  five  most  predominant  patterns  for  winter.  The  summer  frequencies  of 
these  map  patterns  were,  however,  much  lower. 

Although  greater  amounts  of  data  must  be  processed  before  making 
firm  conclusions,  present  indications  suggest  that  the  derived  map  types  as  a 
basic  catalogue  could  be  used  in  other  seasons  as  well.  Map  type 
characteristics,  i.e.,  following  map  types  and  persistence,  may  indicate 
significant  departures  from  those  for  the  summer,  as  would  precipitation 
probabilities. 

6.5  Map  Types  and  Numerical  Weather  Prognoses 

Tests  with  independent  data  in  this  chapter  indicate  that  map 
typing  and  derived  climatological  aids,  namely  PoPs,  show  skill  in 
forecasting.  In  order  for  the  results  of  this  study  to  be  used  in  a  forecast  office, 
numerical  weather  prognoses  (NWP)  of  500  mb.  and  surface  maps  would  be 
required  to  obtain  corresponding  map  types,  to  forecast  for  tomorrow  or  longer 
because  visual  map  typing  is  not  reliable  (section  3.6,  Hartranft  et  al,  1970a). 

With  the  cooperation  of  the  Edmonton  Weather  Office,  a  test  of 
the  numerical  weather  prognoses  was  undertaken.  Upper  air  and  surface  data 
were  abstracted  from  1 :40  million  scale  facsimile  charts  of  the  six-layer  NMC 
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primitive  equations  model.  These  prognoses  were  valid  thirty-six  hours  (at  12Z) 
after  the  initial  data  time.  The  period  used  was  from  April  1,  1974  to  July  31, 
1974,  a  total  of  112  surface  forecasts  and  116  500  mb.  forecasts  because  of  some 
missing  or  unreadable  charts. 

The  forecast  maps  were  correlated  with  the  map  catalogues  and 
verification  of  the  forecasts  was  carried  out  with  the  actual  map  types  which 
developed.  Although  the  primary  and  secondary  map  types  were  obtained,  only 
the  primary  map  types  were  verified.  The  results  are  indicated  in  Table  16. 

TABLE  16 

Forecast  Map  Type  Verification  Using  $ix~Layer 
Thirty-six  Hour  P.E.  Prognosis 


500  mb.  Surface 


Forecast  Primary 

Verifies  43  (37%)  27  (24%) 

Forecast  Uncorrelated 

Verifies  8  (  7%)  19  (16%) 

Forecast  Secondary 
Verifies  as  Actual 

Primary  22  (19%)  11  (10%) 

Forecast  Primary 

Incorrect  43  (37%)  55  (49%) 
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The  most  noticeable  result  is  that  the  NVVP  is  able  to  forecast  the 
500  mb.  map  type  much  better  than  the  surface  type.  This  result  could  be  expected 
since  there  is  greater  variation  and  complexity  in  the  surface  pressure  pattern  and 
the  degree  of  development  of  surface  NVVP.  The  forecast  primary  500  mb.  map 
types  verify  one  and  one-half  times  as  often  as  that  of  the  surface  but  the  variation 
decreases  when  the  correct  uncorrelared  map  types  are  included.  The 
uncorrelated  group  forecasts  for  the  surface  types  are  larger  since  this  is  the 
dominant  group.  There  is  some  correspondence  between  the  forecast  and  the  actual 
map  type  for  the  third  group  in  Table  16  since  the  forecast  secondary  verifies  as 
the  actual  primary  type.  By  combining  the  first  three  groups  as  being  correct 
forecasts,  the  numerical  prognoses  indicate  that  fifty  percent  of  the  maps  will 
verify  correctly  for  the  surface  with  better  results  for  500  mb. 

The  results  are  somewhat  discouraging  but  it  should  be  noted  that  the 
data  sample  was  small  and  the  thirty-six  hour  forecast  period  v/as  relatively  long. 
Further  testing  is  definitely  required  not  only  to  increase  the  data  sample  but  to 
test  the  various  forecast  periods.  It  should  be  emphasized  that  only  the  forecast 
map  types  were  verified.  It  may  be  that  the  forecast  probability  of  precipitation 
verification  will  not  be  changed  to  a  large  degree  even  with  the  errors  in  the 
forecast  type.  This  aspect  needs  to  be  investigated  when  a  iarger  data  sample  is 
obtained . 

All  numerical  weather  prognoses  produce  errors  under  a  variety  of 
synoptic  situations  because  of  simplifying  assumptions,  unrealistic  topographic 
modelling  or  data  initialization .  By  comparing  the  forecast  and  actual  map  types, 
it  may  be  possible  to  distinguish  consistent  errors  in  the  prognosis,  thereby 
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enabling  forecast  errors  to  be  corrected,  not  only  in  the  map  types  but  also  in  the 
weather  occurrence. 
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CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 

7.1  Summary  and  Conclusions 

This  study  has  given  a  means  of  classifying  synoptic  maps  for 
Alberta  using  the  objective  correlation  method  developed  by  Lund  (1963)  and 
later  refined  by  the  Map  Typing  Project  Office,  12th  Weather  Squadron, 
Colorado  Springs.  Classification  of  500  mb.  flows  for  the  summer  period.  May  to 
September  1963-1971,  has  delineated  thirty-five  distinct  map  types  (see  the  500 
mb.  summer  map  catalogue.  Appendix  4).  Statistical  information  on 
occurrence  and  frequency,  map  types  most  frequently  following  and  map  type 
persistence  are  contained  in  the  map  catalogue.  The  500  mb.  catalogue  can  be 
used  as  a  climatological  forecast  aid  for  predominant  summer  flow  patterns  and 
the  occurrence  of  precipitation  in  Western  Canada  and  the  Northwestern  United 
States. 

A  year-round  surface  map  catalogue  was  also  produced  by  the  Map 
Typing  Project  Office,  for  Alberta  and  most  of  British  Columbia.  It  used  data 
covering  a  twenty-six  year  period  (January  to  December,  1946-1971).  Appendix 
2  contains  this  catalogue  of  thirty-three  distinct  surface  map  patterns.  Since  the 
catalogue  contains  statistical  information  on  map  types  for  each  month  of  the 
year,  the  seasonal  change  of  weather  systems  is  indicated. 
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Probabilities  of  twenty-four  hour  summer  (06Z-06Z) 
precipitation  for  both  500  mb.  and  surface  catalogue  types  have  been  derived  and 
they  illustrate  the  catalogues'  application  to  precipitation  forecasting.  Four 
methods  of  analysis  were  employed  for  a  threshold  precipitation  value  of  0.01 
inches  (Appendices  5,  6,  7,  8).  The  results  indicate  that  significant  variation  in 
precipitation  probabilities  are  obtained  with  different  map  types  or 
combinations  of  types.  Skill  scores  indicate  the  best  classification  method  and 
improvement  over  climatological  forecasting.  The  best  results  were  obtained  for 
stations  in  southern  Alberta;  this  is  possibly  related  to  topographic  effects.  In 
general,  results  indicate  low  skill  scores  because  too  many  forecasts  are  made  at 
the  climatic  frequency  of  precipitation.  This  implies  that  more  parameters  are 
required . 

The  derived  probabilities  of  precipitation  were  verified  for  an 
independent  data  sample  (May-September,  1972-1973)  using  the  Brier  Score.  The 
forecast  probabilities  are  highly  reliable  but  low  skill  scores  were  obtained 
because  there  were  too  many  forecasts  near  the  climatological  frequency  of 
precipitation.  Negative  skill  scores  were  obtained  for  many  stations  indicating 
the  need  for  more  research. 

The  derived  precipitation  probabilities  were  also  verified  for  the 
same  period  (May  -  September,  1972-1973)  using  18Z-18Z  precipitation  amounts 
rather  than  06Z-06Z  and  the  00Z  map  type  rather  then  the  12Z.  Comparable 
results  were  obtained. 

A  higher  precipitation  threshold  (0.10  inches  instead  of  0.01 
inches)  was  considered  in  section  6.3  and  concluded  to  be  important  in  light  of 
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a  twenty-four  hour  precipitation  period.  Probabilities  of  twenty-four  hour 
precipitation  (0.10  inches  or  more)  were  therefore  computed  (Appendices  9,  10) 
and  verified  for  1972  and  1973.  The  results  are  comparable. 

The  results  on  precipitation  probabilities  are  good 
considering  that  only  two  predictors,  namely  the  500  mb.  and  surface  flows,  are 
used.  Although  convective  precipitation  predominates  during  summer,  no  attempt 
was  made  to  distinguish  between  the  two  types  in  these  summer  precipitation 
studies. 

Two  additional  studies  were  carried  out.  The  first  consisted  of 
obtaining  map  types  for  a  winter  period  using  the  summer  500  mb.  catalogue. 
Indications  are  that  the  500  mb.  catalogue  types  can  be  used  in  other  seasons. 

The  second  study  tested  the  ability  of  numerical  weather  prognoses  to  forecast  the 
occurrence  of  the  correct  map  type.  The  numerical  prognoses  forecast  the 
occurrence  of  the  correct  type  nearly  firry  percent  of  the  time.  The  500  mb.  flow 
is  forecast  correctly  more  often  than  the  surface  flow.  With  improved  numerical 
prognoses,  the  forecast  verifications  should  increase.  It  is  surmised  that 
consistent  errors  in  the  numerical  prognoses  during  certain  synoptic  situations, 
could  be  revealed  by  map  typing  studies.  This  study  has  shown  the  usefulness  of 
map  types  in  the  verification  of  NWP  prognoses  and  is  currently  being 
implemented  at  the  Edmonton  Weather  Office. 

The  use  of  the  product  moment  correlate  coefficient  to  classify 
weather  maps  was  discussed  in  section  2.2.1 .  The  major  advantage  of  its  use  is 
that  the  classifications  are  determined  solely  by  the  orientation  of  pressure 
systems  and, therefore,  the  directions  of  flow.  It  must  be  remembered  that  the 
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pressure  and  the  pressure  gradient  at  an/  location  is  not  taken  into  account. 

These  two  factors  will  have  a  bearing  on  the  weather  phenomena  being 
investigated.  It  would  be  advantageous  to  define  subsets  within  a  particular  map 
type  based  on  the  magnitude  and  gradients  of  the  weather  system. 

A  further  sophistication  to  map  typing  was  presented  by  Paegle  and 
Kierulff  (1974).  Besides  500  mb.  typing  they  considered  850-500  mb. 
thicknesses,  500-300  mb.  thickness,  500  mb.  absolute  geostrophic  vorticities, 
thermal  advections  and  vertical  velocities.  Their  results  indicate  that  a  similar 
extension  of  this  study  would  be  extremely  beneficial  to  our  understanding  of 
Western  Canada  weather  systems. 

The  classification  criterion  used  in  this  thesis  guarantees  that  each 
map  is  placed  within  the  most  analogous  type.  In  many  cases,  a  particular  map 
may  be  similar  to  two  or  more  predetermined  types.  In  such  situations,  future 
studies  may  find  it  advantageous  to  define  "hybrid"  types,  provided  that  the 
duplicated  classifications  are  numerous. 

The  classification  scheme  does  not  guarantee  that  any  one  map  will 
resemble  the  'key'  map  at  all  locations  or  grid  points.  Paegle  and  Kierulff  (1974) 
found  that  variation  of  maps  belonging  to  a  catalogue  type  resulted  in  a 
significant  categorizing  of  the  atmospheric  fields  over  most  of  the  "window." 

The  derived  precipitation  probabilities  indicated  that 
improvement  over  forecasting  by  climatological  probability  is  obtained  by 

t 

considering  synoptic  pressure  patterns,  a  result  well  known  by 
meteorologists.  The  relatively  low  skill  scores  obtained  may  be  related  to  the 
Twenty-four  hour  precipitation  amounts  which  were  used.  The  data  base  would 
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be  used  more  effectively  if  the  precipitation  period  were  reduced  to  twelve 
hours.  Precipitation  studies  could  also  be  extended  to  account  for  different 
amounts,  types,  and  seasons.  The  consideration  of  too  many  subgroups  of 
precipitation  events  would,  of  course,  develop  statistical  problems  due  to  the 
finite  nature  of  the  data  base.  Some  work  on  manually  recognizing  significant 
subgroups  has  been  done  in  this  study  in  section  5.3.2  and  6.2  and  has  given 
useful  results. 

Although  this  study  has  used  map  types  to  obtain  probabilities  of 
precipitation,  there  are  a  variety  of  other  climatological  aids  which  could  also  be 
produced.  Examples  include  cloud  amounts,  ceiling  and  visibility  tables  and  any 
other  phenomena  for  which  the  synoptic  flow  pattern  is  the  dominant  influence. 

The  500  mb.  and  surface  map  classifications  obtained  in  this  study 
provide  a  good  addition  to  our  knowledge  of  the  weather  systems  in  Western 
Canada.  These  map  catalogues  also  provide  a  good  foundation  for  the  development 
of  climatological  or  forecast  aids  of  weather  phenomena  for  which  the  synoptic 
flow  pattern  is  the  dominant  factor. 


.  1 


.1 

>  •  • 


-  t  •  • 


m  » 

. 


i 


.  \ 


#>  • 


«•  V 


«  » 


BIBLIOGRAPHY 


Augulis,  R.  P.,  1969:  Precipitation  Probabilities  in  the  Western  Region 

Associated  with  Winter  500  mb.  Map  Types.  ESSA  Technical  Memorandum, 
WBTM  WR45-1,  91  pp. 

,  1970:  Precipitation  Probabilities  in  the  Western  Region 
Associated  with  Spring  500  mb.  Map  Types.  ESSA  Technical  Memorandum, 
WBTM  WR  45-2,  75  pp. 

,  1970:  Precipitation  Probabilities  in  the  Western  Region 
Associated  with  Summer  500  mb.  Map  Types.  ESSA  Technical  Memorandum, 
WBTM  WR  45-3,  74  pp. 

,  1970:  Precipitation  Probabilities  in  the  Western  Region 
Associated  with  Fall  500  mb.  Map  Types.  ESSA  Technical  Memorandum, 
WBTM  WR  45-4,  65  pp. 

Barry,  R,  G.,  1960:  A  Note  on  the  Synoptic  Climatology  of  Labrador-Ungava. 
Quart.  J,  Roy.  Meteor.  Soc.,  86,  557-565. 

,  and  A.  H.  Perry,  1973:  Synoptic  Ciimatology  -  Methods  and 
Applications.  Methuen  and  Co.  Ltd.,  London,  555  pp. 

Bjorem,  D.  L.,  1966:  Application  of  Map  Types  to  Fire  Weather  Forecasts  of 
Temperature  and  Relative  Humidity,  M.Sc.  Thesis,  University  of 
Washington,  96  pp. 

Brier,  G.  W.,  1950:  Verification  of  Forecasts  Expressed  in  Terms  of 
Probability.  Mon.  Wea,  Rev.,  78,  1-3. 

,  and  R.  A.  Allen,  1951:  Verification  of  Weather  Forecasts. 
Compendium  of  Meteorology.  Amer.  Meteor.  Soc.,  841-848. 


Burns,  B.  M.,  1970:  Cold  Outbreaks  in  Alberta .  M.Sc„  Thesis,  University  of 
Alberta,  216  pp. 

California  Institute  of  Technology,  1943:  Synoptic  Weather  Types  of  North 
America.  Meteorological  Department,  Calif.  Inst.  Tech.,  Pasadena, 
Calif, ,  243  pp . 


9] 


*  - 


. 


# 


'  * 


- 

-  * 


. 


»  < 


t 


•  * 


■»  ■* 


... 


- 


#  :  i 


92 


Elliot,  R.  D.,  1949:  Forecasting  the  Weather.  Weatherwise,  2,  15-18,  40-43, 
64-67,  86-88,  110-113,  136-138. 

_ j  1951:  Ext  ended -Range  Forecasting  by  Weather  Types. 

Compendium  of  Meteorology,  Amer.  Meteor.  Soc.,  834-840. 

Epstein,  E.S.,  1966:  Quality  Control  for  Probability  Forecasts.  Mon.  Wea.  Rev., 
94,  487-494. 

,  1966:  Point  and  Area  Precipitation  Probabilities.  Mon.  Wea. 
Rev.,  94,  595-598. 


Hage,  K.  D.,  1957:  On  Summer  Cyciogenesis  in  Western  Canada  with  Upper  Cold 
Lows.  Scientific  Report  No.  1,  Contract  %  AF  19(604)  -  2179,  Dept.  of~ 
Meteorology,  University  of  Chicago,  87  pp. 

Hartranft,  F.  R.,  Restivo,  J.  S,  and  R.  C.  Sabin,  1970:  Computerized  Map 
Typing  and  Their  Application  in  t he  Development  of  Forecast  Aids. 

Technical  Paper,  Aerospace  Science  Division,  Headquarters  4th  Weather 
Wing,  No.  70^2,  57  pp. 


,  1970:  Predicting  Heavy 

Snowfall  for  Colorado  Springs  Based  on  Computer  Derived  Synoptic  Map 
Types.  Technical  Paper,  Aerospace  Science  Division,  Headquarters  4th 
Weather  Wing,  No.  70-5,  203  pp. 

Jacobs,  W.  J.,  1946:  Synoptic  Climatology.  Bull.  Amer.  Meteor.  Soc.,  27, 
306-31 1 . 

Knox,  J.  L.,  1969:  The  Use  of  Numerical  Probability  Factors  in  Public  Forecasts 
for  the  Prediction  of  Precipitation  Occurrence.  Technical  Memoranda, 
Department  of  Transport,  Meteorological  Branch,  No.  719,  32  pp. 

Lund,  I.  A.,  1963:  Map-Pattern  Classification  by  Statistical  Methods.  J.  Appl. 
Meteor.,  2,  56-65. 

Map  Typing  Project  Office,  1973:  Map  Type  Catalogue  for  the  Great  Lakes 
Window.  Technical  Paper,  Map  Typing  Project  Office,  12th  Weather 
Squadron,  3rd  Weather  Wing,  No.  73-3,  169  pp. 


,  1973:  Map  Type  Catalogue  for  the  Montana  Window. 
Technical  Paper,  Map  Typing  Project  Office,  12th  Weather  Squadron,  3rd 
Weather  Wing,  No.  73-4,  169  pp. 

,  1973:  Map  Type  Catalogue  for  the  Washington  Window . 
Technical  Paper,  Map  Typing  Project  Office,  12th  Weather  Squadron,  3rd 
WeatTier  Wing,  No.  73-7,  179  pp. 


...  * 


*■  '• 


•  V  r 


.  . 


: 


.  •  : 

'  * 


, 


-  *  .« 


.  /  ‘  *  .  *  e  >1 


.  » 


<v 


- 


*  v  ** 


•  - 

*  V 

.  -  ■  V 


- 


♦  - 


...  I 


*  »".* 

-  * 

•  -  -  ,* 

♦*  •• 


- 


- 


»  *■  •  * 

*  -  I 

K  * 


i*«  m  *• 


* 


- 


93 


Nieman,  R.  L.  and  R.  C.  Sabin,  1972:  The  Application  of  Computer  Map  Typing 
to  Terminal  Forecasting.  Technical  Paper,  Map  Typing  Project  Office,  12th 
Weather  Squadron,  3rd  Weather  Wing,  No.  72-3,  7  pp. 

Nikleva,  S.,  1971:  Verification  of  Probability  of  Precipitation  Occurrence 
Forecasts.  Technical  Memoranda,  Department  of  Transport, 

Meteorological  Branch,  No.  749,  19  pp. 

Paegle,  J.  N.  and  L.  R.  Kierulff,  1974:  Synoptic  Climatology  of  500-mb. 

Winter  Flow  Types.  J.  Appl.  Meteor.,  13,  205-212. 

,  1974:  Prediction  of  Precipitation  Probability 
Based  on  500-mb.  Flow  Types.  J.  Appl.  Meteor.,  13,  213-220. 

Root,  H.  E.,  1962:  Probability  Statements  in  Weather  Forecasting .  J.  Appl. 
Meteor.,  1,  163-168. 

Sabin,  R.  C.,  1972:  Applying  a  Window  Pane  Technique  to  the  Colorado  Springs 
Snow  Study.  Technical  Paper,  Aerospace  Science  Division,  Headquarters 
4th  Weather  Wing,  No.  72-1,  14  pp. 

,  and  R.  L.  Nieman,  1972:  A  Comparison  of  Two  Techniques  Used 
to  Develop  a  Computer  Generated  Map  Type  Catalogue.  Technical  Paper, 
Aerospace  Science  Division,  Headquarters  4th  Weather  Wing,  No.  72-2, 

17  pp. 

Sanders,  R.,  1963:  On  Subjective  Probability  Forecasting,  J.  Appl.  Meteor., 

2,  191-201. 

,  1967:  The  Verification  of  Probability  Forecasts.  J.  Appl. 
Meteor.,  6,  756-761. 

Scholefield,  P.  R.,  1972:  The  Persistence  of  Map  Types.  Technical  Paper,  Map 
Typing  Project  Office,  12th  Weather  Squadron,  3rd  Weather  Wing,  No. 
72-1,  24  pp. 


,  1973:  Comparing  Correlations  Between  Weather  Maps  With 
Similar  Isobcric  Configurations  but  Varying  Pressure  Intensities.  Technical 
Paper,  Map  Typing  Project  Office,  12th  Weather  Squadron,  3rd  Weather 
Wing,  No.  73-8,  11  pp. 

Storr,  D.,  1953:  An  Investigation  into  Means  of  Forecasting  Blizzards  on  the 
Western  Prairies.  Technical  Memoranda,  Department  of  Transport, 
Meteorological  Division,  No.  145,  25  pp. 


*  it 


f  -  -  *  *  .  * 

... 

■  ‘  - 

-  -  • 

*  V 


. 


94 


Thompson,  C.  E.,  1950:  Cold  Lows  Affecting  Alberta  in  Summer  (June  - 
September).  Technical  Memoranda,  Department  of  Transport, 
Meteorological  Division,  No.  76,  14  pp. 

Weatherburn,  C.  E.,  1949:  Mathematical  Statistics.  Cambridge  University  Press, 
London,  271  pp. 

Wold,  H.W.,  1973:  Map  Typing  Computer  Programs.  Technical  Paper,  Map 
Typing  Project  Office,  1 2th  Weather  Squadron,  3rd  Weather  Wing,  No. 
73-6,  170  pp. 


APPENDIX  1 


SURFACE  MAP  TYPE  AREA  AND  GRID  POINT 
PRESSURE  VALUES 


The  surface  map  typing  window  is  illustrated  in  Figure  6,  and  consists  of 
twenty-three  NMC  grid  points  (spacing  381  Km.).  The  window  is  outlined  with 
a  dashed  line.  The  date  and  grid  point  pressure  values  for  each  map  type  are 
given  in  Table  17. 
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TABLE  17  .  GRID  POINT  PRESSURE  VALUES  (CONTINUED) 
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APPENDIX  2 


SURFACE  MAP  TYPES  AND  CHARACTERISTICS 


The  surface  pressure  maps  in  this  Appendix  represent  the  pressure  patterns 
associated  with  the  map  types.  The  window  is  in  the  center  of  the  map,  and  is 
outlined  with  a  dashed  line.  The  characteristics  of  each  map  type  are  summarized 
by  month.  The  number  of  occurrences  indicates  the  number  of  times  the  pressure 
pattern  occurred  during  each  month  for  the  entire  period  of  record  (1946-1971). 
The  frequency  of  occurrences  gives  the  percentage  of  time  that  the  map  type 
occurred  during  the  month.  The  next  set  of  data  indicates  the  map  types  which 
most  frequently  followed  the  given  type  and  the  corresponding  percentage. 
Finally,  the  duration  of  the  map  type  is  indicated  by  the  number  of  twelve  hour 
periods  that  the  map  type  persisted,  and  the  corresponding  percentage. 
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APPENDIX  3 


SUMMER  500  MB.  MAP  TYPE  AREA  AND 
STATION  HEIGHT  VALUES 


The  500  mb.  map  typing  window  is  illustrated  in  Figure  7,  and  consists  of 
nineteen  irregularly  spaced  rawinsonde  stations.  The  date  and  station  height 
values  for  each  map  type  are  given  in  Table  18.  The  station  height  values  are 
expressed  in  decameters  (i.e.,  58  signifies  5,580  meters  . 
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FIGURE  7  • 

500  MB#  MAP  TYPE  AREA 
NUMBERED  STATIONS 
REFER  TO  THE  ORDER 
IN  WHICH  THE  HEIGHT 
VALUES  (DKM.)  ARE 
LISTED  IN  TABLE 
(SEE  BELOW)# 


TABLE  18  .  500  MB# 


MAP  TYPE 

DATE 

500  1 

MB  • 

1 

6 

JUL 

67 

1  2Z 

52 

55 

76 

63 

2 

27 

JUL 

67 

1  2Z 

54 

70 

73 

71 

3 

18 

JUL 

64 

12Z 

62 

60 

88 

70 

4 

23 

JUL 

68 

1  2Z 

74 

76 

76 

69 

5 

24 

JUN 

63 

1  2Z 

38 

44 

75 

66 

6 

26 

JUL 

71 

1  2Z 

77 

83 

73 

76 

7 

25 

JUL 

66 

1  2  Z 

75 

73 

87 

68 

8 

19 

JUL 

70 

1  2  Z 

49 

63 

84 

76 

9 

2 

JUL 

66 

12Z 

6  1 

70 

78 

71 

1  0 

12 

JUN 

69 

1  2Z 

79 

82 

59 

60 

1  1 

13 

JUN 

66 

1  2Z 

55 

57 

65 

60 

12 

13 

JUN 

68 

1  2Z 

65 

55 

65  54 


STATION  HEIGHT  VALUES 

STATION  HEIGHT  VALUES  (DKM.) 


52 

51 

58 

61 

67 

69 

74 

77 

53 

60 

77 

81 

81 

79 

88 

57 

63 

73 

77 

78 

79 

79 

77 

60 

59 

80 

85 

82 

75 

86 

62 

62 

60 

57 

66 

62 

75 

84 

65 

76 

71 

78 

82 

90 

9  1 

72 

69 

71 

75 

72 

76 

77 

77 

60 

62 

81 

85 

84 

80 

86 

37 

49 

61 

61 

68 

68 

65 

71 

63 

75 

69 

71 

65 

89 

76 

79 

77 

83 

87 

84 

82 

84 

80 

63 

59 

87 

88 

87 

66 

84 

71 

68 

66 

67 

62 

63 

66 

79 

67 

72 

73 

78 

8  1 

88 

89 

53 

71 

74 

70 

82 

79 

86 

85 

77 

76 

84 

88 

90 

78 

91 

65 

67 

67 

67 

64 

59 

65 

75 

68 

76 

5C 

59 

69 

82 

83 

76 

72 

74 

79 

73 

76 

72 

66 

61 

41 

76 

75 

71 

58 

69 

59 

58 

60 

63 

66 

70 

72 

69 

57 

58 

79 

86 

82 

62 

76 

63 

59 

56 

50 

55 

52 

57 

64 

48 

48 

59 

69 

68 

72 

75 

- 


. 

• 

' 

- 

.  - 

' 
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TABLE  18.  500  MB,  STATION  HEIGHT  VALUES  (CONTINUED) 

MAP  TYPE  DATE  500  MB,  STATION  HEIGHT  VALUES  ( DKM • ) 


13 

17 

JUN 

65 

1  2Z 

46 

64 

48 

68 

75 

77 

73 

73 

7  1 

69 

78 

77 

71 

83 

71 

67 

65 

36 

7  l 

14 

4 

SEP 

68 

0  z 

66 

78 

61 

57 

71 

82 

74 

83 

74 

65 

59 

59 

56 

53 

86 

89 

78 

64 

68 

1  5 

15 

SEP 

69 

oz 

56 

59 

53 

54 

55 

62 

57 

64 

59 

60 

64 

49 

55 

58 

71 

76 

72 

74 

77 

16 

22 

MAY 

71 

1  2Z 

43 

58 

43 

52 

63 

70 

68 

73 

65 

59 

63 

62 

53 

52 

72 

71 

59 

65 

6  1 

17 

4 

JUN 

67 

OZ 

73 

75 

68 

61 

66 

71 

63 

67 

63 

66 

63 

57 

50 

51 

68 

71 

70 

71 

77 

18 

24 

MAY 

69 

12Z 

62 

61 

62 

65 

66 

57 

7  1 

62 

74 

75 

72 

69 

56 

53 

67 

70 

79 

67 

8  1 

19 

2 

JUL 

64 

12Z 

57 

66 

57 

66 

71 

67 

73 

67 

74 

80 

80 

72 

65 

63 

67 

70 

79 

79 

87 

20 

15 

SEP 

67 

12Z 

52 

77 

56 

81 

88 

87 

90 

87 

87 

82 

73 

88 

76 

75 

85 

85 

84 

62 

78 

21 

21 

AUG 

68 

12Z 

50 

67 

53 

59 

62 

65 

59 

62 

60 

68 

75 

61 

61 

68 

61 

62 

65 

82 

76 

22 

29 

MAY 

67 

1  2  Z 

47 

44 

52 

57 

52 

40 

55 

45 

58 

63 

66 

63 

64 

66 

52 

58 

62 

69 

64 

23 

5 

JUN 

67 

1  2Z 

58 

70 

61 

60 

69 

75 

7  1 

71 

66 

70 

71 

64 

50 

56 

68 

68 

69 

74 

76 

24 

24 

AUG 

63 

1  2Z 

77 

75 

77 

68 

71 

69 

68 

64 

69 

79 

77 

66 

62 

71 

61 

67 

76 

8  1 

83 

25 

9 

MAY 

71 

OZ 

39 

45 

38 

39 

50 

59 

60 

64 

62 

60 

64 

54 

29 

63 

67 

64 

63 

68 

62 

26 

1  1 

JUN 

68 

OZ 

55 

51 

60 

64 

60 

52 

64 

61 

70 

69 

66 

66 

62 

64 

67 

71 

76 

64 

72 

27 

12 

JUL 

65 

oz 

72 

77 

75 

70 

71 

72 

66 

70 

68 

73 

78 

66 

60 

68 

70 

73 

74 

84 

83 

28 

9 

JUN 

65 

oz 

56 

63 

62 

65 

72 

67 

73 

70 

74 

76 

69 

67 

52 

4  1 

73 

73 

75 

65 

78 

29 

26 

MAY 

65 

12Z 

58 

68 

57 

65 

70 

73 

68 

70 

67 

55 

51 

62 

57 

50 

71 

72 

65 

45 

57 

30 

8 

JUL 

69 

OZ 

60 

74 

60 

64 

71 

76 

71 

76 

7  1 

68 

70 

65 

65 

72 

76 

78 

74 

77 

79 

31 

7 

MAY 

66 

1  2Z 

52 

61 

59 

70 

69 

62 

64 

66 

66 

70 

72 

67 

67 

57 

72 

75 

78 

72 

82 

32 

1  1 

SEP 

70 

12Z 

80 

70 

73 

56 

55 

68 

6  1 

73 

63 

59 

60 

50 

50 

46 

79 

85 

77 

66 

77 

33 

22 

JUL 

70 

12Z 

6  1 

70 

61 

53 

63 

72 

63 

71 

64 

71 

77 

65 

50 

71 

73 

77 

74 

8  1 

8  1 

34 

30 

JUN 

68 

OZ 

74 

74 

73 

73 

71 

74 

66 

7  1 

63 

57 

58 

63 

66 

56 

71 

76 

64 

69 

64 

35 

16 

AUG 

63 

1  2Z 

78 

72 

76 

77 

73 

71 

75 

74 

78 

80 

76 

70 

73 

70 

80 

85 

85 

78 

87 

APPENDIX  4 


SUMMER  500  MB.  MAP  TYPES  AND  CHARACTERISTICS 


The  500  mb.  contour  maps  in  this  Appendix  represent  the  contour  patterns 
associated  with  the  map  types.  The  window  is  contained  within  the  numbered 
stations.  The  characteristics  of  each  map  type  are  summarized  by  month.  The 
number  of  occurrences  indicates  the  number  of  times  the  contour  pattern  occurred 
during  each  month  for  the  entire  period  of  record  (May  -  September  1963-1971). 
The  frequency  of  occurrences  gives  the  percentage  of  time  that  the  map  type 
occurred  during  the  month.  The  next  set  of  data  indicates  the  map  types  which 
most  frequently  followed  the  given  type  and  the  corresponding  percentage. 
Finally,  the  duration  of  the  map  type  is  indicated  by  the  number  of  twelve-hour 
periods  that  the  map  type  persisted,  and  the  corresponding  percentage. 
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*************** 

♦  500  MB*  * 

*  MAP  TYPE  1  * 
*************** 


**** 

PERIOD  OF  RECORD 
1 963-197 l 
**** 

THIS  TYPE 

OCCURRED  525  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 


MAY  JUNE  JULY 


AUG  SEPT 


NUMBER  OF 


OCCURRENCES 

91 

81 

132 

1  19 

102 

FREQUENCIES  OF 

OCCURRENCES 

16 

15 

24 

21 

19 

MAP  TYPE  MOST 

1/ 

62 

1/ 

63 

1/ 

67 

1  / 

67 

1/ 

64 

FREQUENTLY 

3/ 

13 

3/ 

14 

3/ 

9 

2/ 

8 

3/ 

1  3 

FOLLOWING 

2/ 

6 

2/ 

5 

2/ 

4 

3/ 

5 

2/ 

10 

(MAP  TYPE/ 

l  5/ 

4 

5/ 

5 

4/ 

4 

25/ 

5 

4/ 

3 

FREQUENCY) 

25/ 

4 

4/ 

4 

2  5/ 

3 

4/ 

4 

8/ 

3 

MOST  FREQUENT 

1/ 

34 

2/ 

29 

1/ 

36 

1  / 

46 

1/ 

44 

DURATION  (NUMBER 

2/ 

31 

3/ 

29 

2/ 

22 

2/ 

1  3 

2/ 

26 

OF  1 2  HR  PERIODS/ 

3/ 

1  1 

1/ 

26 

3/ 

13 

5/ 

1  0 

3/ 

8 

FREQUENCY) 

4/ 

9 

6/ 

6 

4/ 

7 

4/ 

8 

4/ 

8 

MAP  TYPE  1 


. 


H 


- 

• 

. 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  2  * 
*************** 


**** 

PERIOD  OF  RECORD 
1 963— 1 97 1 
**** 

THIS  TYPE 

OCCURRED  385  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

72 

45 

93 

80 

95 

FREQUENCIES  OF 

OCCURRENCES 

13 

8 

17 

14 

18 

MAP  TYPE  MOST 

2/ 

64 

2/ 

60 

2/ 

70 

2/ 

59 

2/ 

61 

FREQUENTLY 

1/ 

1  5 

8/ 

18 

1/ 

1  3 

1  / 

1  3 

1/ 

1  3 

FOLLOWING 

8/ 

8 

1/ 

1  1 

8/ 

7 

8/ 

10 

8/ 

6 

(MAP  TYPE/ 

6/ 

3 

16/ 

4 

1  1/ 

3 

28/ 

4 

1  1/ 

5 

FREQUENCY) 

1  1/ 

3 

4/ 

2 

4/ 

2 

4/ 

3 

4/ 

3 

MOST  FREQUENT 

1/ 

3  1 

1/ 

56 

2/ 

31 

1/ 

44 

1/ 

50 

DURATION  (NUMBER 

2/ 

31 

2/ 

1  1 

1/ 

2; 

3/ 

21 

2/ 

13 

OF  1 2  HR  PERIODS/ 

5/ 

12 

3/ 

1  1 

5/ 

1  7 

2/ 

18 

3/ 

1  1 

FREQUENCY) 

3/ 

8 

8/ 

1  1 

3/ 

1  0 

4/ 

6 

4/ 

1  1 

i 


MAP  TYPE  2 


■  *3* 


. 


. 

r  •  ’  .  ' 

s 

-  ■  . 
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*************** 

♦  5C0  MB.  * 

♦  MAP  TYPE  3  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  249  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

4 

2 

49 

65 

59 

34 

FREQUENCIES  OF 

OCCURRENCES 

8 

9 

12 

11 

6 

MAP  TYPE  MOST 

3/ 

52 

3/ 

67 

3/ 

67 

3/ 

68 

3/ 

48 

FREQUENTLY 

7/ 

12 

U/ 

1  0 

1/ 

13 

15/ 

1  0 

7/ 

1  5 

FOLLOWING 

1/ 

10 

1/ 

4 

7/ 

8 

1  / 

8 

1  5/ 

9 

(MAP  TYPE/ 

1  5/ 

1  0 

7/ 

4 

15/ 

3 

U/ 

5 

21/ 

6 

FREQUENCY) 

U/ 

10 

5/ 

2 

21/ 

3 

7/ 

3 

3  3/ 

6 

MOST  FREQUENT 

1/ 

55 

2/ 

29 

1/ 

32 

3/ 

32 

2/ 

44 

DURATION  (NUMBER 

2/ 

20 

1/ 

24 

2/ 

32 

1  / 

21 

1/ 

39 

OF  1 2  HR  PERIODS/ 

3/ 

10 

3/ 

24 

3/ 

9 

2/ 

21 

3/ 

1  1 

FREQUENCY) 

6/ 

10 

5/ 

12 

8/ 

9 

4/ 

1  1 

5/ 

6 

MAP  TYPE  3 


•>.  ‘  *  *  <• 


- 


. 

. 

• 

*1  • 
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*************** 

*  500  MB •  * 

*  MAP  TYPE  4  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  126  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

14 

16 

45 

35 

16 

FREQUENCIES  OF 

OCCURRENCES 

3 

3 

8 

/ 

6 

3 

MAP  TYPE  MOST 

4/ 

71 

4/ 

40 

4/ 

60 

4/ 

49 

4/ 

38 

FREQUENTLY 

2/ 

7 

2/ 

20 

2/ 

12 

2/ 

23 

1  4/ 

1  9 

FOLLOWING 

l  2/ 

7 

1/ 

13 

1/ 

5 

6/ 

6 

1/ 

6 

(MAP  TYPE/ 

1  7/ 

7 

1  4/ 

7 

18/ 

5 

18/ 

6 

2/ 

6 

FREQUENCY ) 

U/ 

7 

15/ 

7 

2  8/ 

5 

23/ 

6 

1  1/ 

6 

MOST  FREQUENT 

1  / 

25 

1/ 

70 

1/ 

47 

1/ 

39 

1/ 

60 

DURATION  (NUMBER 

3/ 

25 

2/ 

20 

2/ 

21 

2/ 

33 

2/ 

30 

OF  1 2  HR  PERIODS/ 

4/ 

25 

5/ 

1  0 

3/ 

21 

3/ 

22 

4/ 

10 

FREQUENCY) 

6/ 

25 

5/ 

5 

4/ 

6 

MAP  TYPE  4 


* 


' 


4 

• 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  5  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  58  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 

MONTH  MAY  JUNE  JULY  AUG  SEPT 

NUMBER  OF 


OCCURRENCES 

9 

28 

4 

4 

13 

FREQUENCIES  OF 

OCCURRENCES 

2 

5 

1 

1 

2 

MAP  TYPE  MOST 

5/ 

44 

5/ 

43 

5/ 

50 

5/  50 

5/ 

62 

FREQUENTLY 

25/ 

22 

U/ 

21 

8/ 

25 

3/  25 

1/ 

15 

FOLLCWI NG 

1  / 

1  1 

1/ 

7 

30/ 

25 

8/  25 

2  5/ 

15 

(MAP  TYPE/ 

8/ 

1  1 

3/ 

7 

8/ 

8 

FREQUENCY) 

U/ 

1  1 

2  5/ 

7 

MOST  FREQUENT 

1/ 

40 

1/ 

63 

1/ 

50 

2/100 

2/ 

40 

DURATION  (NUMBER 

2/ 

40 

2/ 

19 

3/ 

50 

1/ 

20 

OF  1 2HR  PERIODS/ 

3/ 

20 

3/ 

6 

3/ 

20 

FREQUENCY) 

4/ 

6 

5/ 

20 

MAP  TYPE  5 


. 


- 

. 

• 

* 

- 

% 

.... 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  6  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  65  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 

MONTH  MAY  JUNE  JULY  AUG  SEPT 

NUMBER  OF 


OCCURRENCES 

2 

1 

10 

14 

10 

10 

FREQUENCIES  OF 

OCCURRENCES 

4 

2 

3 

2 

2 

MAP  TYPE  MOST 

6/ 

62 

6/ 

50 

6/ 

71 

6/ 

57 

6/ 

40 

FREQUENTLY 

U/ 

14 

28/ 

40 

2/ 

14 

28/ 

22 

2/ 

30 

FGLLOWI NG 

28/ 

10 

U/ 

10 

28/ 

7 

U/ 

1  1 

1  4/ 

10 

<  MAP  TYPE/ 

1  1/ 

5 

29/ 

7 

28/ 

10 

FREQUENCY ) 

1  6/ 

5 

U/ 

10 

MOST  FREQUENT 

1/ 

38 

1/ 

40 

1/ 

25 

1  / 

25 

1/ 

67 

DURATION  (NUMBER 

2/ 

25 

2/ 

40 

2/ 

2  L 

2/ 

25 

2/ 

17 

OF  12HR  PERIODS/ 

3/ 

25 

4/ 

20 

3/ 

25 

3/ 

25 

4/ 

17 

FREQUENCY) 

8/ 

13 

8/ 

25 

4/ 

25 

MAP  TYPE  6 


*  ' 

- 

• 

- 

• 

* 

- 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  7  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  40  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

5 

4 

21 

3 

7 

FREQUENCIES  OF 

OCCURRENCES 

1 

1 

4 

1 

1 

MAP  TYPE  MOST  27/ 

40 

U/ 

50 

7/ 

43 

7/ 

33 

21/ 

29 

FREQUENTLY  U/ 

40 

7/ 

25 

1  5/ 

14 

27/ 

33 

U/ 

29 

FOLLOWING  1/ 

20 

27/ 

25 

4/ 

1  0 

33/ 

33 

7/ 

14 

(MAP  TYPE/ 

2  7/ 

1  0 

9/ 

1  4 

FREQUENCY) 

9/ 

5 

27/ 

1  4 

MOST  FREQUENT  l/l 

00 

1/ 

67 

1/ 

58 

1  / 

50 

1/ 

83 

DURATION  (NUMBER 

2/ 

33 

2/ 

1  7 

2/ 

50 

2/ 

17 

OF  1 2  HR  PERIODS/ 

3/ 

17 

FREQUENCY) 

4/ 

8 

MAP  TYPE  7 


♦ 


-  -  .  ' 

- 

•  . 

* 

• 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  8  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-197 1 
**** 

THIS  TYPE 

OCCURRED  187  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 


MAY  JUNE  JULY 


AUG  SEPT 


NUMBER  OF 
OCCURRENCES 


38  29  36 


52 


32 


FREQUENCIES  OF 


OCCURRENCES 

7 

MAP  TYPE  MOST 

8/ 

53 

8/ 

FREQUENTLY 

1/ 

18 

1/ 

FOLLOWI NG 

22/ 

8 

5/ 

(MAP  TYPE/ 

25/ 

5 

2/ 

FREQUENCY ) 

26/ 

5 

MOST  FREQUENT 

2/ 

39 

1/ 

DURATION  (NUMBER 

1/ 

28 

2/ 

OF  1 2HR  PERIODS/ 

3/ 

28 

3/ 

FREQUENCY ) 

4/ 

6 

4/ 

5 

6 

9 

6 

59 

8/ 

67 

8/ 

67 

3/ 

53 

28 

1/ 

1  9 

1  / 

22 

1/ 

28 

1  0 

2/ 

3 

3/ 

8 

5/ 

9 

3 

3/ 

3 

26/ 

4 

3/ 

6 

1  9/ 

3 

25/ 

3 

42 

2/ 

33 

1  / 

28 

1/ 

47 

25 

3/ 

2; 

3/ 

22 

3/ 

20 

8 

1/ 

17 

2/ 

17 

4/ 

20 

8 

6/ 

1  7 

5/ 

1  7 

2/ 

13 

MAP  TYPE  8 


-  5 


*  M  .•> 

. 


’  ' 

-  •  •  ~ 


145 


*************** 

*  500  MB.  * 

*  MAP  TYPE  9  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-197 1 
**** 

THIS  TYPE 

OCCURRED  19  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

2 

5 

9 

1 

2 

FREQUENCIES  OF 

OCCURRENCES 

0 

1 

2 

0 

0 

MAP  TYPE  MOST 

U/l  00 

9/ 

60 

9/ 

56 

21/100 

21/  50 

FREQUENTLY 

3/ 

20 

2  2/ 

22 

U/  50 

FOLLOWING 

24/ 

20 

1  9/ 

1  1 

(MAP  TYPE/ 

21/ 

1  1 

FREQUENCY) 

MOST  FREQUENT 

1/100 

1/ 

50 

3/ 

50 

1/100 

1/100 

DURATION  (NUMBER 

4/ 

50 

1/ 

25 

OF  12HR  PERIODS/ 

2/ 

25 

FREQUENCY ) 


MAP  TYPE  9 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  10  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  27  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 

MONTH  MAY  JUNE  JULY  AUG  SEPT 

NUMBER  OF 


OCCURRENCES 

7 

1  4 

0 

0 

6 

FREQUENCIES  OF 

OCCURRENCES 

1 

3 

0 

0 

1 

MAP  TYPE  MOST 

l  0/ 

71 

10/ 

79 

6/ 

33 

FREQUENTLY 

6/ 

14 

U/ 

14 

1  0/ 

33 

FOLLOWING 

U/ 

1  4 

14/ 

7 

14/ 

17 

{MAP  TYPE/ 

17/ 

17 

FREQUENCY) 

MOST  FREQUENT 

1/ 

50 

2/ 

33 

1/ 

50 

DURATION  (NUMBER 

6/ 

50 

3/ 

33 

2/ 

50 

OF  12HR  PERIODS/ 

9/ 

33 

FREQUENCY ) 


MAP  TYPE  10 


*  ♦ 


J 


■ 


i  .  • 


:  -  ■ 


* 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  11  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  75  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 


MAY  JUNE  JULY 


AUG  SEPT 


NUMBER  OF 


OCCURRENCES 

1  6 

23 

5 

10 

21 

FREQUENCIES  OF 

OCCURRENCES 

3 

4 

1 

2 

4 

MAP  TYPE  MOST 

2/ 

20 

1  1/ 

73 

1/  40 

1 1  / 

50 

1  1/ 

43 

FREQUENTLY 

1  1/ 

20 

1/ 

9 

2/  40 

u/ 

20 

1/ 

24 

FOLLOWI  NG 

1/ 

13 

2/ 

9 

8/  20 

1  / 

l  0 

1  4/ 

10 

(MAP  TYPE/ 

14/ 

13 

8/ 

5 

8/ 

1  0 

2/ 

5 

FREQUENCY) 

4/ 

7 

26/ 

5 

18/ 

10 

4/ 

5 

MOST  FREQUENT 

1/ 

77 

1/ 

57 

1/100 

2/ 

60 

1/ 

58 

DURATION  (NUMBER 

2/ 

23 

4/ 

1  4 

1/ 

20 

2/ 

25 

OF  1 2  HR  PERIODS/ 

> 

7/ 

14 

3/ 

20 

3/ 

8 

FREQUENCY ) 

8/ 

14 

5/ 

8 

MAP  TYPE  11 


. 


• 
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*************** 

*  *  500  MB*  * 

*  MAP  TYPE  12  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  28  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERI ST  ICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 
OCCURRENCES 

6 

3 

1  1 

6 

2 

FREQUENCIES  OF 
OCCURRENCES 

1 

1 

2 

1 

0 

MAP  TYPE  MOST 

12/ 

50 

12/  33 

12/  36 

U/  50 

7/  50 

FREQUENTLY 

1  / 

17 

15/  33 

3/  18 

12/  33 

12/  50 

FOLLOWING 

27/ 

17 

35/  33 

7/  18 

4/  17 

(MAP  TYPE/ 

U/ 

17 

U/  18 

FREQUENCY) 

4/  9 

MOST  FREQUENT 

1/ 

33 

1/  50 

2/  57 

1/75 

2/100 

DURATION  (NUMBER 

2/ 

33 

2/  50 

1/  43 

3/  25 

OF  12HR  PERIODS/ 
FREQUENCY ) 

3/ 

33 

MAP  TYPE  12 


149 


*************** 

*  500  MB.  * 

*  MAP  TYPE  13  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  15  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 


MAY  JUNE  JULY  AUG  SEPT 


NUMBER  OF 
OCCURRENCES 


0  14 


0  0  1 


FREQUENCIES  OF 

OCCURRENCES  0  3 


0  0 


0 


MAP  TYPE  MOST 
FREQUENTLY 
FOLLOWING 
(MAP  TYPE/ 
FREQUENCY)- 

MOST  FREQUENT 
DURATION  (NUMBER 
OF  12HR  PERIODS/ 
FREQUENCY ) 


13/  50 
5/  29 
U/  21 


1/  43 
2/  43 
5/  14 


U/l  00 


1/100 


MAP  TYPE  13 


• 
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******  ********* 

*  500  MB*  * 

*  MAP  TYPE  14  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  59  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

1  2 

6 

1 

9 

31 

FREQUENCIES  OF 

OCCURRENCES 

2 

1 

0 

2 

* 

6 

MAP  TYPE  MOST 

1  4/ 

36 

6/ 

33 

6/100 

14/ 

67 

l  4/ 

70 

FREQUENTLY 

2/ 

1  8 

14/ 

33 

2/ 

1  1 

2/ 

10 

FOLLOWING 

1  1/ 

18 

2/ 

1  7 

6/ 

1  1 

6/ 

10 

(MAP  TYPE/ 

6/ 

9 

U/ 

17 

U/ 

1  1 

10/ 

3 

FREQUENCY) 

16/ 

9 

29/ 

3 

MOST  FREQUENT 

1/ 

63 

1/ 

75 

1/100 

2/ 

67 

2/ 

30 

DURATION  (NUMBER 

2/ 

25 

3/ 

25 

5/ 

33 

1/ 

20 

OF  12HR  PERIODS/ 

3/ 

13 

3/ 

20 

FREQUENCY)  4/  10 


MAP  TYPE  14 


- 


, 


•  . 

. 

- 

-  •  • 

- 

I 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  15  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  66  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 


MAY  JUNE  JULY 


AUG  SEPT 


NUMBER  OF 


OCCURRENCES 

12 

1  0 

17 

16 

1  1 

FREQUENCIES  OF 

OCCURRENCES 

2 

2 

3 

3 

2 

MAP  TYPE  MOST 

1  5/ 

25 

15/ 

40 

15/ 

41 

15/ 

25 

1  5/ 

45 

FREQUENTLY 

1/ 

17 

35/ 

20 

4/ 

24 

2/ 

19 

1/ 

9 

FOLLOWING 

1  1/ 

17 

1/ 

10 

1/ 

1  8 

1  / 

1  3 

4/ 

9 

(MAP  TYPE/ 

U/ 

1  7 

4/ 

1  0 

2/ 

12 

4/ 

13 

1  1/ 

9 

FREQUENCY) 

2/ 

8 

14/ 

10 

35/ 

6 

1  1  / 

6 

1  4/ 

9 

MOST  FREQUENT 

1/ 

73 

1/ 

50 

1  / 

50 

1  / 

67 

1/ 

50 

DURATION  (NUMBER 

2/ 

1  1 

2/ 

33 

2/ 

30 

2/ 

33 

3/ 

33 

OF  1 2  HR  PERIODS/ 

3/ 

11 

3/ 

1  7 

3/ 

20 

2/ 

17 

FREQUENCY) 


MAP  TYPE  15 


.... 

•- 


-  • 


r' 


* 

* 

' 

.  . 
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*************** 

♦  500  MB,  * 

*  MAP  TYPE  16  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  51  TIMES 
DURING  THE  PERIDD 
OF  RECORD 


CHARACTERI STICS 


MONTH 


MAY  JUNE  JULY  AUG  SEPT 


NUMBER  OF 

OCCURRENCES  26  8  0  2  15 

FREQUENCIES  OF 


OCCURRENCES 

5 

1 

0  0 

3 

MAP  TYPE  MOST 

l  6/ 

64 

16/ 

38 

2/  50 

16/ 

40 

FREQUENTLY 

2/ 

16 

2/ 

25 

11/50 

2/ 

27 

FOLLOWING 

8/ 

12 

8/ 

25 

1  1/ 

1  3 

(MAP  TYPE/ 

2  5/ 

8 

2  5/ 

13 

25/ 

13 

FREQUENCY) 

8/ 

7 

MOST  FREQUENT 

2/ 

50 

2/ 

60 

1/100 

1/ 

56 

DURATION  (NUMBER 

1/ 

20 

1/ 

40 

2/ 

22 

OF  1 2  HR  PERIODS/ 

3/ 

10 

3/ 

22 

FREQUENCY) 

4/ 

10 

M*P  TYPE  16 


*  . 


- 


. 


: 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  17  ♦ 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  24  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 


MAY  JUNE  JULY  AUG  SEPT 


NUMBER  OF 


OCCURRENCES 

4 

4 

2 

0 

14 

FREQUENCIES  OF 

OCCURRENCES 

1 

4 

1 

0 

0 

3 

MAP  TYPE  MOST 

1  7/ 

50 

1  7/ 

50 

4/  50 

l  7/ 

43 

FREQUENTLY 

28/ 

25 

23/ 

25 

17/  50 

32/ 

29 

FOLLOWING 

U/ 

25 

U/ 

25 

4/ 

14 

(MAP  TYPE/ 

6/ 

7 

FREQUENCY) 

U/ 

7 

MOST  FREQUENT 

1/ 

50 

1/ 

50 

2/100 

1  / 

63 

DURATION  (NUMBER 

3/ 

50 

3/ 

50 

2/ 

1  3 

OF  1 2  HR  PERIODS/ 

3/ 

13 

FREQUENCY) 

4/ 

13 

MAP  TYPE  17 


y  - 


► 

' 


* 

•  * 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  18  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  35  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 

MONTH  MAY  JUNE  JULY  AUG  SEPT 

NUMBER  OF 


OCCURRENCES 

4 

6 

5 

17 

3 

FREQUENCIES  OF 

OCCURRENCES 

1 

1 

1 

3 

1 

MAP  TYPE  MOST 

3/  50 

18/ 

33 

1/ 

20 

18/ 

59 

2/ 

33 

FREQUENTLY 

18/  50 

26/ 

33 

2/ 

20 

1  / 

1  2 

18/ 

33 

FOLLOWING 

U/ 

33 

3/ 

20 

3/ 

12 

19/ 

33 

(MAP  TYPE/ 

19/ 

20 

4/ 

6 

FREQUENCY) 

35/ 

20 

12/ 

6 

MOST  FREQUENT 

2/100 

1/ 

50 

1/ 

1  00 

1/ 

43 

1/ 

50 

DURATION  (NUMBER 

2/ 

50 

2/ 

29 

2/ 

50 

OF  1 2  HR  PERIODS/ 

• 

5/ 

29 

FREQUENCY) 

MAP  TYPE  18 


' 

• 

*  .  \ ' 

* 

• 

- 

- 
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*  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

*  500  MB.  * 

*  MAP  TYPE  19  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  55  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

8 

4 

15 

19 

9 

FREQUENCIES  OF 

OCCURRENCES 

1 

1 

3 

3 

2 

MAP  TYPE  MOST 

3/ 

25 

9/ 

25 

19/ 

47 

19/ 

61 

2/ 

33 

FREQUENTLY 

19/ 

25 

19/ 

25 

1/ 

20 

1/ 

6 

3/ 

22 

FOLLOWING 

21/ 

25 

21/ 

25 

3/ 

20 

2/ 

6 

19/ 

22 

( MAP  TYPE/ 

31/ 

13 

31/ 

25 

5/ 

7 

3/ 

6 

5/ 

1  1 

FREQUENCY) 

U/ 

13 

1  5/ 

7 

5/ 

6 

18/ 

1  1 

MOST  FREQUENT 

1/ 

83 

1/ 

67 

1/ 

50 

1/ 

50 

1/ 

71 

DURATION  (NUMBER 

3/ 

17 

2/ 

33 

3/ 

38 

3/ 

25 

2/ 

29 

OF  1 2  HR  PERIODS/ 

2/ 

13 

4/ 

13 

FREQUENCY) 

5/ 

13 

MAP  TYPE  19 


f 


--T  I 

/. 


•  V 
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*************** 

*  500  MB*  * 

*  MAP  TYPE  20  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  18  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 


MAY  JUNE  JULY  AUG  SEPT 


NUMBER  OF 

OCCURRENCES  3  2  0  3  10 


FREQUENCIES  OF 

OCCURRENCES  1  0  0 


MAP  TYPE  MOST 
FREQUENTLY 
FOLLOWI NG 
(MAP  TYPE/ 
FREQUENCY) 

MOST  FREQUENT 
DURATION  (NUMBER 
OF  1 2  HR  PERIODS/ 
FREQUENCY) 


11/  33  2/  50 

20/  33  20/  50 

U/  33 


1/  50  2/100 

2/  50 


U/  67  20/  80 

20/  33  16/  10 

U/  10 


1/  50  1/  50 
2/  50  9/  50 


MAP  TYPE  20 


. 


- 


. 

* 


157 


*************** 

*  500  MB.  * 

*  MAP  TYPE  21  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  27  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

4 

4 

7 

5 

7 

FREQUENCIES  OF 

OCCURRENCES 

1 

1 

1 

1 

1 

MAP  TYPE  MOST 

5/ 

50 

U/ 

50 

21/ 

57 

21  / 

60 

25/ 

29 

FREQUENTLY 

9/ 

25 

5/ 

25 

30/ 

29 

5/ 

20 

1/ 

14 

FOLLOWING 

21/ 

25 

13/ 

25 

15/ 

14 

25/ 

20 

2/ 

14 

(MAP  TYPE/ 

16/ 

14 

FREQUENCY) 

2  1/ 

1  4 

MOST  FREQUENT 

1/ 

67 

1/1 

00 

2/ 

67 

2/ 

50 

1/ 

83 

DURATION  (NUMBER 

2/ 

33 

3/ 

33 

3/ 

50 

2/ 

17 

OF  1 2HR  PERIODS/ 
FREQUENCY) 


MAP  TYPE  21 


.  •  .  .  . 

* 


* 


- 


' 

•  • 

•  . 
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**♦♦** **♦* ***** 

♦  500  MB.  * 

♦  MAP  TYPE  22  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  26  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

15 

2 

7 

2 

0 

FREQUENCIES  OF 

OCCURRENCES 

3 

0 

1 

0 

0 

MAP  TYPE  MOST 

22/ 

53 

5/  50 

22/ 

43 

22/100 

FREQUENTLY 

U/ 

27 

U/  50 

3/ 

29 

FOLLOWI NG 

3/ 

7 

9/ 

14 

(MAP  TYPE/ 

5/ 

7 

U/ 

14 

FREQUENCY) 

9/ 

7 

MOST  FREQUENT 

1/ 

43 

1/100 

1/ 

50 

2/100 

DURATION  (NUMBER 

3/ 

29 

2/ 

25 

OF  1 2HR  PERIOOS/ 

2/ 

14 

3/ 

25 

FREQUENCY) 

4/ 

14 

MAP  TYPE  22 


<  : 


2!  -■'u 


.* 


- 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  23  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  17  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

3 

3 

2 

6 

3 

FREQUENCIES  OF 

OCCURRENCES 

1 

1 

0 

1 

1 

MAP  TYPE  MOST 

2/  33 

23/  67 

2/100 

2/ 

33 

23/  67 

FREQUENTLY 

1 6/  33 

2/  33 

4/ 

17 

17/  33 

FOLLOWING 

28/  33 

23/ 

17 

(MAP  TYPE/ 

28/ 

17 

FREQUENCY) 

U/ 

17 

MOST  FREQUENT 

1/  100 

3/100 

1/100 

1/ 

80 

3/100 

DURATION  (NUMBER 

2/ 

20 

OF  1 2 HR  PERIODS/ 
FREQUENCY! 


MAP  TYPE  23 


•  r  »  • 


■ 


' 


•  •  i  * 
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*************** 

♦  500  MB.  * 

*  MAP  TYPE  24  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  15  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

1 

2 

4 

7 

1 

FREQUENCIES  OF 

OCCURRENCES 

0 

0 

1 

1 

0 

MAP  TYPE  MOST 

U/  100 

24/  50 

7/ 

25 

24/  83 

U/100 

FREQUENTLY 

U/  50 

9/ 

25 

27/  17 

FOLLOWING 

24/ 

25 

(MAP  TYPE/ 

27/ 

25 

FREQUENCY) 

MOST  FREQUENT 

1/100 

2/100 

1/ 

67 

1/  50 

1/100 

DURATION  (NUMBER 

2/ 

33 

6/  50 

OF  12HR  PERIODS/ 
FREQUENCY) 


MAP  TYPE  24 


,  -  V 

- 


t  • 


* 
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*************** 

*  500  MB.  * 

♦  MAP  TYPE  25  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  64  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERI STICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

14 

8 

9 

1 1 

22 

FREQUENCIES  OF 

OCCURRENCES 

3 

1 

2 

2 

4 

MAP  TYPE  MOST 

2/ 

36 

U/ 

50 

1/ 

44 

1/ 

30 

2  5/ 

45 

FREQUENTLY 

25/ 

29 

2  5/ 

25 

2  5/ 

33 

8/ 

20 

2/ 

25 

FOLLOWING 

1/ 

14 

1/ 

13 

2/ 

1  1 

25/ 

20 

1/ 

20 

(MAP  TYPE/ 

8/ 

7 

4/ 

13 

19/ 

1  1 

4/ 

1  0 

8/ 

5 

FREQUENCY) 

33/ 

7 

19/ 

10 

19/ 

5 

MOST  FREQUENT 

1/ 

70 

1/ 

83 

1/ 

67 

1  / 

78 

1/ 

46 

DURATION  (NUMBER 

2/ 

20 

3/ 

17 

2/ 

17 

2/ 

22 

2/ 

38 

OF  1 2  HR  PERIODS/ 

3/ 

10 

3/ 

17 

3/ 

15 

FREQUENCY ) 


MAP  TYPE  25 


• 

. .  • 

* 

162 


***♦**♦♦***♦♦** 

*  500  MB.  * 

♦  MAP  TYPE  26  * 

*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  20  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

4 

10 

1 

5 

0 

FREQUENCIES  OF 

OCCURRENCES 

1 

2 

0 

1 

0 

MAP  TYPE  MOST 

U/ 

50 

26/ 

40 

1/100 

3/  40 

FREQUENTLY 

1/ 

25 

1/ 

20 

26/  40 

FOLLOWING 

26/ 

25 

3/ 

20 

8/  20 

(MAP  TYPE/ 

18/ 

10 

FREQUENCY) 

U/ 

10 

MOST  FREQUENT 

1/ 

67 

1/ 

67 

1/100 

2/  67 

DURATION  (NUMBER 

2/ 

33 

2/ 

17 

1/  33 

OF  1 2HR  PERIODS/ 
FREQUENCY) 

4/ 

17 

MAP  TYPE  26 


- 


♦  * 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  27  * 
*************** 


**** 

PER  100  OF  RECORD 
1963-1971 
♦  *** 

THIS  TYPE 

OCCURRED  17  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

3 

4 

5 

4 

1 

FREQUENCIES  OF 

OCCURRENCES 

1 

1 

1 

1 

0 

MAP  TYPE  MOST 

23/ 

67 

27/  75 

2  7/ 

40 

4/ 

25 

17/100 

FREQUENTLY 

4/ 

33 

33/  25 

4/ 

20 

15/ 

25 

FOLLOWI NG 

24/ 

20 

27/ 

25 

IMAP  TYPE/ 

FREQUENCY) 

U/ 

20 

U/ 

25 

MOST  FREQUENT 

1/100 

4/100 

2/ 

67 

1  / 

67 

1/100 

DURATION  (NUMBER 

1/ 

33 

2/ 

33 

OF  12HR  PERIODS/ 
FREQUENCY ) 


MAP  TYPE  27 


- 

-  ’ 


* 

«• 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  28  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  36  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CH ARACTERI STICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

6 

7 

5 

14 

4 

FREQUENCIES  OF 

OCCURRENCES 

1 

1 

1 

3 

1 

MAP  TYPE  MOST 

28/ 

33 

28/ 

43 

2/  50 

28/ 

36 

1/ 

25 

FREQUENTLY 

2/ 

17 

4/ 

29 

1/  25 

2/ 

21 

2/ 

25 

FOLLOWING 

6/ 

17 

2/ 

1  4 

18/  25 

U/ 

14 

8/ 

25 

{MAP  TYPE/ 

18/ 

17 

U/ 

14 

4/ 

7 

28/ 

25 

FREQUENCY 1 

23/ 

17 

8/ 

7 

MOST  FREQUENT 

1/ 

75 

1/ 

50 

1/100 

1/ 

56 

1/ 

67 

DURATION  (NUMBER 

3/ 

25 

2/ 

25 

2/ 

33 

2/ 

33 

OF  1 2 HR  PERIODS/ 

3/ 

25 

3/ 

1  1 

FREQUENCY ) 


MAP  TYPE  28 


* 


■ 


• 

V  • 

i  .  .  :  . 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  29  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  13  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

NUMBER  OF 

OCCURRENCES 

6 

FREQUENCIES  OF 

OCCURRENCES 

1 

MAP  TYPE  MOST 

29/ 

50 

FREQUENTLY 

6/ 

17 

FOLLOWING 

20/ 

17 

(MAP  TYPE/ 

U/ 

17 

FREQUENCY) 

• 

MOST  FREQUENT 

1/ 

67 

DURATION  (NUMBER 

4/ 

33 

OF  12HR  PERIODS/ 

FREQUENCY) 

* 

JUNE  JULY  AUG  SEPT 

0  3  0  4 

0  „  1  0  1 

6/  67  29/  75 

29/  33  U/  25 

1/  50  4/100 

2/  50 


MAP  TYPE  29 


' 

* 


•• 

*- 

- 

‘ 


- 

- 

•  ■  . 

• 

*  • 
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**  ************* 

*  500  MB.  * 

*  MAP  TYPE  30  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  18  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

3 

5 

8 

1 

1 

FREQUENCIES  OF 

OCCURRENCES 

1 

1 

1 

0 

0 

MAP  TYPE  MOST 

2/  67 

3  0/ 

40 

2/ 

50 

2/100 

16/100 

FREQUENTLY 

U/  33 

16/ 

20 

30/ 

38 

FOLLOWING 

25/ 

20 

5/ 

13 

(MAP  TYPE/ 

U/ 

20 

FREQUENCY) 

MOST  FREQUENT 

1/100 

2/ 

67 

1/ 

80 

1/100 

1/1  00 

DURATION  (NUMBER 

1/ 

33 

4/ 

20 

OF  1 2 HR  PERIODS/ 
FREQUENCY ) 


MAP  TYPE  30 


‘ 


, 

• 

* 

'  -  . 
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*************** 

*  500  MB.  * 

*  MAP  TYPE  31  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 

♦  *** 

THIS  TYPE 

OCCURRED  13  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

4 

4 

3 

2 

0 

FREQUENCIES  OF 

OCCURRENCES 

1 

1 

1 

0 

0 

MAP  TYPE  MOST 

31/  50 

19/ 

50 

31/  67 

8/100 

FREQUENTLY 

21/  25 

1/ 

25 

11/  33 

FOLLOWING 

(MAP  TYPE/ 
FREQUENCY) 

U/  25 

31/ 

25 

MOST  FREQUENT 

2/100 

1/ 

67 

3/100 

1/100 

DURATION  (NUMBER 
OF  1 2HR  PERIODS/ 

2/ 

33 

FREQUENCY ) 


MAP  TYPE  31 


* 


168 


*************** 

*  500  MB.  * 

*  MAP  TYPE  32  * 
*************** 


*♦♦* 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  15  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 


MAY  JUNE  JULY 


NUMBER  OF 

OCCURRENCES  000 


FREQUENCIES  OF 

OCCURRENCES  000 


MAP  TYPE  MOST 
FREQUENTLY 
FOLLOWING 
(MAP  TYPE/ 
FREQUENCY) 

MOST  FREQUENT 
DURATION  C  NUMBER 
OF  1 2  HR  PERIODS/ 
FREQUENCY) 


AUG  SEPT 


0  15 


0 


3 


32/ 

33 

17/ 

27 

U/ 

20 

10/ 

13 

34/ 

7 

1/ 

60 

2/ 

30 

3/ 

10 

MAP  TYPE  32 


- 


r5 


•  i  T 

•  - 


169 


*************** 

♦  500  MB.  * 

*  MAP  TYPE  33  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 
**  ** 

THIS  TYPE 

OCCURRED  19  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 

MONTH  MAY  JUNE  JULY  AUG  SEPT 

NUMBER  OF 


OCCURRENCES 

3 

4 

5 

4 

3 

FREQUENCIES  OF 

OCCURRENCES 

1 

1 

1 

1 

1 

MAP  TYPE  MOST 

1/ 

33 

1/ 

25 

33/ 

40 

1  / 

50 

2/ 

33 

FREQUENTLY 

4/ 

33 

17/ 

25 

2/ 

20 

33/ 

25 

2  5/ 

33 

FOLLOWI NG 

30/ 

33 

25/ 

25 

4/ 

20 

U/ 

25 

33/ 

33 

(MAP  TYPE/ 

30/ 

25 

25/ 

20 

FREQUENCY) 

MOST  FREQUENT 

1/100 

1/1 

00 

2/ 

67 

1  / 

67 

1/ 

50 

DURATION  (NUMBER 

1/ 

33 

2/ 

33 

2/ 

50 

OF  12HR  PERIOOS/ 
FREQUENCY) 


MAP  TYPE  33 


• 

• 
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** ****** ******* 

*  500  MB*  * 

*  MAP  TYPE  34  * 
♦****♦♦******♦♦ 


**** 

PERIOD  OF  RECORD 
1963-1971 
**** 

THIS  TYPE 

OCCURRED  7  TIMES 

DURING  THE  PERIOD 
OF  RECORD 


CHARACTERI STICS 


MONTH 


MAY  JUNE  JULY  AUG  SEPT 


NUMBER  OF 

OCCURRENCES  0200 


FREQUENCIES  OF 

OCCURRENCES  0000 


MAP  TYPE  MOST  34/100 

FREQUENTLY 

FOLLOWING 

(MAP  TYPE/ 

FREQUENCY) 

MOST  FREQUENT  2/100 

DURATION  (NUMBER 
OF  1 2 HR  PERIODS/ 

FREQUENCY) 


34/  60 
10/  20 
U/  20 


2/  50 
3/  50 


MAP  TYPE  34 


.  .  .  •  . 

* 


. 


.  , 

' 


171 


*************** 

*  500  MB.  * 

*  MAP  TYPE  35  * 
*************** 


**** 

PERIOD  OF  RECORD 
1963-1971 

**** 

THIS  TYPE 

OCCURRED  18  TIMES 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

■ 

OCCURRENCES 

3 

5 

6 

4 

0 

FREQUENCIES  OF 

OCCURRENCES 

1 

1 

1 

1 

0 

MAP  TYPE  MOST 

12/ 

67 

11/ 

40 

1  2/ 

50 

35/  75 

FREQUENTLY 

35/ 

33 

12/ 

20 

7/ 

17 

12/  25 

FOLLOWING 

35/ 

20 

35/ 

1  7 

(MAP  TYPE/ 

U/ 

20 

U/ 

17 

FREQUENCY) 

MOST  FREQUENT 

1/ 

50 

1/ 

75 

1/ 

80 

4/100 

DURATION  (NUMBER 

2/ 

50 

2/ 

25 

2/ 

20 

OF  1 2  HR  PERIODS/ 
FREQUENCY) 


MAP  TYPE  35 


• 

* 
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********************* 

♦  500  MB.  * 

♦  UNCORRELATED  MAPS  * 
********************* 


*♦** 

PERIOD  OF  RECORD 
1963-1971 
**** 

322  MAPS  WERE 
UNCORRELATED 
DURING  THE  PERIOD 
OF  RECORD 


CHARACTERISTICS 


MONTH 

MAY 

JUNE 

JULY 

AUG 

SEPT 

NUMBER  OF 

OCCURRENCES 

97 

119 

18 

48 

40 

FREQUENCIES  OF 

OCCURRENCES 

1  7 

22 

3 

9 

7 

MAP  TYPE  MOST 

U/ 

67 

U/ 

66 

U/ 

44 

U/ 

55 

U/ 

54 

FREQUENTLY 

1  4/ 

5 

1/ 

3 

2/ 

6 

14/ 

6 

1  6/ 

8 

FOLLOWING 

16/ 

4 

6/ 

3 

4/ 

6 

2/ 

4 

32/ 

8 

(MAP  TYPE/ 

6/ 

2 

1  0/ 

3 

7/ 

6 

3/ 

4 

14/ 

5 

FREQUENCY) 

8/ 

2 

13/ 

3 

8/ 

6 

4/ 

4 

20/ 

5 

MOST  FREQUENT 

1/ 

37 

1/ 

44 

1/ 

60 

1  / 

55 

1/ 

58 

DURATION  (NUMBER 

2/ 

20 

2/ 

23 

2/ 

20 

3/ 

18 

2/ 

16 

OF  1 2  HR  PERIODS/ 

3/ 

17 

4/ 

9 

3/ 

10 

2/ 

1  4 

3/ 

1  1 

FREQUENCY) 

4/ 

1  1 

3/ 

7 

5/ 

1  0 

4/ 

5 

4/ 

1  1 

UNCORRELATED  MAPS 


APPENDIX  5 


SUMMER  500  MB.  MAP  TYPE  PoPs 

The  twenty-four  hour  percent  probabilities  of  precipitation  (  >0.01  inches) 
are  listed  for  each  500  mb.  map  type.  The  number  of  cases  used  to  calculate  the 
probabilities  are  indicated  for  each  map  type. 

The  station  locations  are  shown  in  Figure  5  (page  52). 
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APPENDIX  6 


SUMMER  SURFACE  MAP  TYPE  PoPs 

The  twenty-four  hour  percent  probabilities  of  precipitation  (  >  0.01  inches) 
are  listed  for  each  surface  map  type.  The  number  of  cases  used  to  calculate  the 
probabilities  are  indicated  for  each  map  type. 

The  station  locations  are  shown  in  Figure  5  (page  52  ). 
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APPENDIX  7 


SUMMER  500  MB.  -  SURFACE  MAP  TYPE  PoPs 

The  twenty-four  hour  percent  probabilities  of  precipitation  (  >  0.01  inches) 
are  listed  for  each  500  mb.  -  surface  combination.  The  number  of  cases  used  to 
calculate  the  probabilities  are  indicated  for  each  map  type  combination  (only  > 
ten  cases  are  retained).  The  surface  type  labelled  OTHER  includes  all  surface 
map  types  which  had  fewer  than  ten  occurrences  initially.  This  classification  is  to 
be  used  when  the  surface  type  does  not  correspond  to  the  significant  combinations. 

The  station  locations  are  shown  in  Figure  5  (page  52  ). 
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APPENDIX  8 


SUMMER  500  MB.  -  500  MB.  MAP  TYPE  PoPs 


The  twenty-four  hour  percent  probabilities  of  precipitation  (  >  0.01  inches) 
are  listed  for  each  500  mb.  -  500  mb.  combination.  The  combinations  use  the 
initial  500  mb.  map  type  and  the  following  twelve-hour  500  mb.  map  type,  to 
distinguish  weather  system  motions.  The  number  of  cases  used  to  calculate  the 
probabilities  are  indicated  for  each  map  type  combination  (only  ^  ten  cases  are 
retained).  The  following  twelve-hour  500  mb.  map  type  labelled  OTHER  includes 
all  500  mb.  map  types  which  had  fewer  than  ten  occurrences  initially.  This  last 
classification  is  to  be  used  when  the  following  500  mb.  map  does  not  correspond 
to  any  significant  500  mb.  map  type. 

The  station  locations  are  shown  in  Figure  5  (page  52  ). 
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APPENDIX  9 


SUMMER  500  MB.  MAP  TYPE  PoPs  (  >0.10  INCHES) 

Twenty-four  hour  percent  probabilities  of  precipitation  (  >  0.10  inches)  are 
listed  for  each  500  mb.  map  type.  The  number  of  cases  used  to  calculate  the 
probabilities  are  indicated  for  each  map  type. 

The  station  locations  are  shown  in  Figure  5  (page  52  ). 
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APPENDIX  10 


SUMMER  SURFACE  MAP  TYPE  PoPs  (  >  0.10  INCHES) 

Twenty-four  hour  percent  probabilities  of  precipitation  (  >0.10  inches) 
are  listed  for  each  surface  map  type.  The  number  of  cases  used  to  calculate  the 
probabilities  are  indicated  for  each  map  type. 

The  station  locations  are  shown  in  Figure  5  (page  52  ). 
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APPENDIX  11 


THE  BRIER  SCORE 


One  essential  criterion  for  a  successful  verification  method  is  that 
improving  the  average  score  is  consistent  with  a  better  forecast  of  the  specific 
situation.  If  this  is  not  done,  on  occasion,  it  will  turn  out  that  the  "best" 
forecast  on  the  basis  of  score  is  not  the  most  useful  forecast.  The  Brier  Score 
meets  the  above  criterion  and  is  explained  in  an  article  by  Brier  and  Allen  (1951). 

For  a  single  forecast,  the  Brier  Score  can  be  written: 

B  =  (F-E)2  (11-1) 


where  F  equals  the  forecast  probability  of  occurrence  which  ranges  from  zero  to 
one  and  E  equals  zero  if  the  event  does  not  occur  or,  one  if  the  event  occurs. 

The  value  of  B  range  from  zero  to  one  with  zero,  the  best  score. 

The  Brier  Score  becomes  more  meaningful  when  averaged  over  a  large  number 
of  forecasts . 


Let 


N  =  number  of  forecasts 


th 

F.  =  forecast  probability  on  the  i  forecast 

th 

E.  =  1  if  event  occurs  on  the  i  occasion 

i 

th 

=  0  if  non-occurrence  on  the  i  occasion 


The  Brier  Score  for  the  N  forecasts  (B  )  is: 

n 


N 


Br 


lE  [F,-E,]: 


i :  1 


(11-2) 
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* 


•  , 


213 


Suppose  now  that  all  probabilities  F.  are  expressed  in  tenths  and  that  the 
set  of  N  forecasts  are  classified  into  subsets  according  to  the  value  of  F .  ,  the 
forecast  probabilities.  There  will  then  be  eleven  subsets  as  follows: 


F. 


0.0,  0.1,  0.2,  .  .  .0.9,  1.0 


to  be  designated  by  k  where 


k  =  1,  2,  3,  .  .  .  11. 

Consider  a  subset  "k11  which  is  made  up  of  forecasts  (e.g.,  if  0.6 
chance  of  rain  was  forecast  on  fifteen  occasions  then  =  15)  and  a  total  of 
event  occurrences  (e.g.,  if  rain  occurred  in  the  fifteen  forecasts  ten  times, 
=  10).  The  a  posteriori,  or  empirical  probability  will  be 

R. 


0, 


(11-3) 


M, 


for  the  example. 


0. 


10/15  =  0.67. 


From  equation  11-2,  the  contribution  to  the  Brier  Score  for  the  event 

2 

occurrences  is  1/M  ^  (F^  -  1)  R^  and  the  contribution  for  the  days  for  the 

2  th 

non-occurrences  is  1/M ^  (F^  -  0)  (M^  -  R^).  The  Brier  Score  for  the  k 


subset  is  then  given  by 
B, 


(F k  -  l)20k  +  F2k(l  -  0L) 


k  v'  k  ,/vk  '  k  v  ’  ~k 

The  mean  Brier  Score,  considering  all  subsets,  will  be 


(11-4) 


B 


=  ttE  MkB 


k<i 


w 


ii 

here  N  =  M 


k  :  1 


(11-5) 
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Equation  11-4  may  be  simplified  by  adding  and  subtracting  0^  and 


rearranging  terms  to  obtain 


B, 


(Fk  -  V  *  0k  (1  -  V 


01-6) 


with  the  mean  Brier  Score  becoming 


B  N  ^  MktFk  '  V  +  0k(1  '  V} 


(11-7) 


The  quantity  (F^  -  0^)  measures  how  close  to  the  actual  frequency 
0^  ,  the  group  of  forecasts  with  probability  predicted  the  event  to  occur. 
Sanders  (1963)  calls  this  term  a  measure  of  the  "validity"  of  the  forecasts,  while 
other  authors,  e.g.,  Root  (1962),  refer  to  it  as  the  "reliability."  To  minimize 
the  reliability  term,  the  forecast  probability  F^  should  equal  the  relative 
frequency  0^  . 

This  term  usually  turns  out  to  be  rather  small.  It  is  the  second  term  of 
equation  11-6  which  is  the  main  contributor  to  the  Brier  Score. 

The  quantity  0^(1  -  0^)  is  a  measure  of  the  "sharpness"  (Sanders,  1963) 
or  "resolution"  (Root,  1962)  of  the  forecasts.  To  demonstrate  this  term,  suppose 
that  the  forecaster  has  been  able  to  assess  his  probabilities  so  that  F^  always 
equals  0^  ,  i.e.,  perfect  reliability.  The  Brier  Score  then  becomes 
B, 


V1  -  V 


and  this  term  reaches  its  minimum,  zero,  only  when 


0, 


0 


or 


0, 


1  , 


k  w  k 

that  is,  when  all  instances  have  been  sorted  into  two  categories,  in  one 

which  the  event  occurs  and  in  the  other  of  which  the  event  does  not  occur. 


•  ■  4  fll  m  i  *  -  - 


' 

:W  ■  -  •  '  «•  ]' 
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The  term  reaches  its  maximum  value,  0.25,  when  0^  =  0.5. 

Remembering  the  assumption  that  the  forecaster  is  able  to  make  =  0^  , 
it  is  thus  seen  that  to  minimize  the  resolution  term,  the  forecast  probability  should 
be  as  close  as  possible  to  the  extreme  ends  of  the  probability  scale. 

The  reliability  of  forecasts  is  illustrated  in  Figure  8,  by  plotting,  for  each 
probability  interval,  the  forecast  probability  against  the  observed  rain 
occurrence.  Included  in  the  graph  is  the  number  of  forecasts  issued  in  each 
probability  category. 

Since  the  diagonal  line  represents  perfect  reliability,  (observed  frequency 

=  forecast  probability),  the  vertical  distance  between  the  diagonal  and  the 

plotted  line  (F  ^  -  0^)  indicates  the  amount  of  error  in  the  respective 

2 

probability  forecast.  Since  (F^  -  0^)  is  the  squared  deviation,  it  is, 
therefore,  a  measure  of  the  reliability  of  the  forecasts  in  the  subset. 

The  root-mean-square  error  (RMSE)  in  reliability  for  all  probability 
categories  can  therefore  be  defined  as 

RMSE  =  {i00/n}  £  Mk(Fk“0k)2  (,1’8) 


If  there  is  a  major  difference  in  the  climatic  frequency  of  particular 
weather  phenomena  between  one  locality  and  another,  the  respective  Brier  Scores 
will  not  be  comparable.  In  an  attempt  to  overcome  this,  a  skill  score  can  be 
defined  as 


100  (1  -  B/cB) 


01-9) 


where  cB  is  the  Brier  Score  obtained  by  employing  the  climatological  frequency 

as  the  forecast  probability,  and  B  equals  the  Brier  Score  obtained  by  using  the 
forecast  probabilities. 


' 


/  • 
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FIGURE  8 

Verification  of  Precipitation  Probability  Forecasts  Using  Primary  500  mb. 
and  Surface  Map  Types  (May  -  September,  1963-1971,  06Z-06Z) 


FORECAST  PRECIPITATION  PROBABILITY  (%) 


• 


*> 


(11-10) 


The  climatological  Brier  Score  is  defined  by 
cB  c  -  2clV'N  +  R/N 

where  c  is  the  climatic  frequency,  R  is  the  total  number  of  event 
occurrences  and  N  is  the  total  number  of  forecasts. 

The  skill  score,  thus  defined,  provides  a  measure  of  the  percentage 
improvement  of  the  forecast  probabilities  over  the  use  of  the  climatological 
frequency  as  the  probability  forecast. 


